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ABSTRACT 

i 

A  major  goal  in  the  study  of  plate  tectonics  is  the  acqui¬ 
sition  of  a  knowledge  of  the  history  of  relative  motion  among 
the  rigid  plates  of  the  earth's  lithosphere L_^The  three  papers 
of  this  thesis <contr4bute  to  this  effort  and>demonstrate  that 
studies  of  the  stability  and  evolution  of  triple  junctions  and 
of  the  finite  rotations  of  systems  of  three  plates  can  yield 
significantly  more  accurate  tectonic  histories  than  can  studies 
of  -the^f relative  motions  between  two  plates  alone.  Topographic 
and  magnetic  investigation  of  the  Southwest  Indian  Ridge  and 
reconstruction  of  the  plate  system  of  the  Indian  Ocean  shows  - 

that  both  Africa  and  Antarctica  are  rigid  plates  and  their  pole  '•Jcq 
of  relative  rotation  has  remained  fixed  near  8^N,  42J|W  since  ^ 

the  Eocene.  A  detailed  survey  of  the  Indian  Ocean  triple  ’ 

junction  reveals  that  the  Indian  Ocean  plate  motions  have 
remained  constant  since  10  Ma .  stability  conditions  of  the 

junction  show  that  the  general  morphology  of  the  Southwest 
Indian  Ridge  results  from  the  evolution  of  the  Indian  Ocean 
triple  junction.  A  method  is  presented  for  determining  the 
finite  rotations  best  reconstructing  the  past  relative  positions 
of  three  plates  around  a  triple  junction.  The  method  is 
illustrated  by  reconstructions  of  the  plates  around  the  Labrador 
Sea  triple  junction  at  the  times  of  anomalies  24  (56  Ma)  and 
21  (50  Ma).y^The  region  of  uncertainty  of  the  Greenland-North 
America  finite  pole  is  mapped  for  each  reconstruction,  and  it 
demonstrates  bhat  consideration  of  the  three  plate  system 
yields  more  well-constrained  results  than  does  a  treatment  of 
the  two  plates  alone. 


Thesis  Supervisor:  J.  G.  Sclater 
Title:  Professor 
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One  of  the  earliest  and  most  important  goals  in  the  study 
of  the  new  global  tectonics  has  been  to  determine  the  relative 
positions  of  the  earth’s  lithospheric  plates  at  various  times 
in  the  past.  Initially,  such  plate  reconstructions  were  for  the 
purpose  of  demonstrating  the  plausibility  of  the  continental 
drift  hypothesis  (e.g.,  Wegener,  1912;  Du  Toit,  1937;  Bullard 
et  al . ,  1965),  but  they  have  since  come  to  form  much  of  the 
basis  for  research  in  geology  and  in  the  physics  of  the  earth. 

An  accurate  and  precise  knowledge  of  the  history  of  relative 
motions  among  the  lithospheric  plates  is  imperative,  since  for 
several  hundred  million  years,  these  motions  have  been  the 
principal  mechanism  shaping  the  face  of  the  planet. 

The  history  of  sea-floor  spreading,  for  example,  has 
almost  entirely  determined  the  thermal  evolution  of  the  oceanic 
lithosphere  and  thus  the  general  topography  of  the  ocean  basins 
(see  discussion  in  Parsons  and  Sclater,  1977) .  From  plate 
reconstructions  one  may  then  discover  the  three-dimensional 
shape  of  the  ocean  basins  as  a  function  of  time  (Sclater  et  al . , 
1977) .  This,  in  turn,  has  heavily  influenced  paleo-oceanographic 
conditions  and  the  history  of  sedimentation  (van  Andel  et  al . , 
1977) .  Plate  reconstructions  thus  provide  both  a  background 
against  which  to  view  paleo-oceanographic  and  sedimentological 
studies  and  a  powerful  tool  with  which  to  pursue  them.  Further¬ 
more,  the  plate  motions  are  the  most  concrete  evidence  against 
which  theories  of  mantle  convection  and  of  the  driving  forces 
of  plates  must  be  tested  (Forsyth  and  Uyeda,  1975;  Solomon  et  al . , 
1975;  Solomon  et  al.,  1977). 
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The  bold  outlines  of  sea-floor  spreading  are  now 
established  over  most  of  the  earth,  but  the  scientific  consi¬ 
derations  above  and  the  related  search  for  economic  resources 
in  the  ocean  basins  call  for  ever  more  detailed  analyses  of 
plate  motions.  This  thesis  is  a  collection  of  three  papers 
which  contribute  to  this  effort.  In  particular,  they  focus 
on  the  value  of  studying  triple  junctions  and  systems  of  three 
or  more  plates  in  order  to  determine  both  instantaneous  and 
finite  relative  plate  rotations  with  the  greatest  possible 
accuracy  and  precision.  The  geometry  and  evolution  of  a  triple 
junction  are  highly  dependent  on  the  local  relative  motions 
among  the  three  plates,  and  thus  the  sea-floor  formed  near  a 
triple  junction  is  a  particularly  sensitive  recorder  of  any 
changes  in  the  relative  velocities.  Moreover,  the  relative 
rotations  within  a  system  of  three  rigid  plates  are  interdepen¬ 
dent.  This  is  true  for  both  the  instantaneous  and  finite  cases 
and  can  provide  important  additional  constraints  not  available 
within  a  two  plate  system.  Therefore,  one  may  more  precisely 
define  the  relative  plate  motions. 

*  *  * 

The  tectonic  evolution  of  the  southern  oceans  (Figure  1) 
is  the  least  well-known  of  all  the  earth's  oceanic  areas.  The 
region  is  vast,  remote,  and  relatively  poorly  surveyed.  It  is 
also  complex,  involving  a  large  number  of  major  plates  and 
having  experienced  severe  changes  of  spreading  poles  and  even 
of  plate  boundaries.  An  understanding  of  the  southern  oceans 
is  vital,  however,  to  a  general  understanding  of  the  global 
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FIGURE  1 


Chart  showing  the  major  tectonic  features  of  the  Indian  and 
South  Atlantic  Oceans.  From  Norton  and  Sclater  (in  press). 
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tectonic  history.  Knowledge  of  the  relative  motions  among 
Africa,  India,  and  Antarctica  is  necessary  in  order  to  directly 
connect  the  history  of  the  Atlantic  plate  system  with  that  of 
the  Pacific,  all  other  paths  being  broken  by  subduction  zones. 
Moreover,  the  regional  tectonic  evolution  has  had  a  dominant 
effect  on  paleo-oceanographic  conditions  as  the  surrounding 
continents  separated  from  Antarctica,  fostering  the  formation 
of  a  circumpolar  oceanic  current  system.  Even  at  present, 
regional  tectonic  conditions  influence  the  movement  of  water 
masses,  since  deep  fracture  zones  across  the  Southwest  Indian 
Ridge  provide  pathways  by  which  Antarctic  bottom  water  reaches 
the  western  Indian  Basin  (Warren,  1978) . 

Chapter  2  of  this  thesis  presents  the  results  of  a  topo¬ 
graphic  and  magnetic  study  of  the  eastern  portion  of  the 
Southwest  Indian  Ridge,  the  boundary  between  the  African  and 
Antarctic  plates,  an  area  not  heretofore  well  surveyed.  Because 
of  the  slow  spreading  rate  of  less  than  10  mm/yr,  magnetic 
anomalies  are  poorly  resolved,  yet  our  careful  survey  allows 
us  to  determine  the  pole  and  rate  of  relative  motion  across 
the  ridge.  This  pole  is  consistent  with  data  from  farther 
westward  (Sclater  et  a_l. ,  1976;  Sclater  et  al . ,  1978;  Bergh 
and  Norton,  1976),  indicating  that  both  the  African  and  Antarctic 
plates  are  behaving  as  rigid  bodies.  Consideration  of  the 
three  plate  system  of  the  Indian  Ocean  using  published  sea-floor 
spreading  data  and  finite  rotations  for  the  other  plate  pairs 
then  shows  that  the  Africa-Antarctica  relative  rotation  has 
been  constant  since  the  Eocene  to  within  a  temporal  resolution 
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of  about  10  m.y. .  Moreover,  the  dominant  morphologic  trends 
along  the  eastern  portion  of  the  Southwest  Indian  Ridge  are 
a  direct  result  of  the  evolution  of  the  Indian  Ocean  triple 
junction . 

*  *  * 

Since  McKenzie  and  Morgan  (1969)  pointed  out  the  sensitive 
dependence  of  the  stability  and  evolution  of  triple  junctions 
on  the  motions  of  the  plates  around  them,  several  investigators 
have  studied  triple  junctions  in  order  to  better  determine  the 
instantaneous  relative  motions  of  lithospheric  plates  (Falconer, 
1972;  Sclater  et  al^. ,  1976;  Hey,  1977).  In  addition  to  the 
sensitivity  mentioned  above,  these  efforts  are  aided  by  an 
important  possible  redundancy  of  sea-floor  spreading  data.  If 
one  assumes  the  plates  are  locally  rigid,  then  any  four  of  the 
six  quantities  describing  plate  motion  about  a  triple  junction 
(three  rates  and  three  directions  of  relative  motion)  are 
sufficient  to  determine  the  other  two.  If  one  is  able  to 
measure  all  six  quantities,  then  one  can  test  the  local 
rigidity  of  the  plates. 

Chapter  3  presents  results  of  a  study  of  the  evolution 
of  the  Indian  Ocean  triple  junction  near  25°S,  70°E.  We  are 
able  to  determine  rates  and  directions  on  all  three  plate 
boundaries  and  show  that  the  velocity  triangle  is  closed  and 
the  plates  are  locally  rigid.  The  plate  velocities  in  the 
vicinity  of  the  triple  junction  are  consistent  with  those 
elsewhere  along  the  Indo-Atlantic  plate  boundaries,  and  we 
determine  the  poles  of  instantaneous  relative  rotation.  A 
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careful  mapping  of  the  evolution  of  the  triple  junction  shows 
that  these  poles  have  remained  constant  since  10  Ma. 

*  *  * 

Since  Bullard  et  al.  (1965)  first  used  a  statistical 
criterion  to  fit  together  the  continental  shelves  around  the 
Atlantic  Ocean,  several  methods  have  been  developed  to  find  by 
similar  means  the  best  fit  of  crustal  isochron  patterns  as 
defined  by  sea-floor  spreading  data  (e.g.,  McKenzie  and  Sclater, 
1971;  Pilger,  1978) .  Recently,  Hellinger  (1979)  has  presented 
a  new  technique  for  determining  the  finite  rotation  defining  the 
best  reconstruction  of  the  past  relative  positions  of  two  plates. 
The  assumptions  about  the  data  implicit  in  the  method  are  well- 
suited  to  a  wide  class  of  sea-floor  spreading  data,  and,  signi¬ 
ficantly,  the  errors  in  the  data  contribute  rigorously  to  a 
statistically  defined  region  of  confidence  about  the  optimal 
rotation  tensor. 

In  Chapter  4,  Hellinger' s  technique  is  extended  to  treat 
the  relative  finite  rotations  of  a  system  of  three  plates  about 
a  triple  junction.  Dealing  with  a  three  plate  system  has  a 
number  of  advantages  over  dealing  with  pairs  of  plates  alone. 

For  example,  the  additional  constraints  imposed  by  simultaneously 
seeking  reconstructions  of  three  plate  pairs  might  overcome 
difficulties  arising  from  a  paucity  of  data  along  one  or  more 
of  the  boundaries.  Moreover,  it  is  important  that  any  proposed 
set  of  finite  rotations  for  a  reconstruction  be  internally 
consistent.  This  method  insures  that  this  condition  is  satis¬ 
fied.  If  no  such  satisfactory  set  can  be  found,  one  must 


-s~~s3araBgaB 


21 


re-examine  the  assumption  of  rigid  plate  behavior. 

The  method  is  illustrated  by  application  to  the  plates 
around  the  now  extinct  Labrador  Sea  triple  junction.  Using 
available  magnetic  anomaly  identifications  (Srivastava,  1978; 
Phillips  and  Tapscott,  in  preparation) ,  we  reconstruct  the 
relative  positions  of  Greenland,  North  America,  and  Rockall 
Bank  at  the  times  of  anomalies  21  (50  Ma)  and  24  (56  Ma)  and 
determine  the  relative  finite  rotations  of  the  plates.  These 
finite  rotations  yield  a  better  reconstruction  than  do  those 
of  other  authors.  They  are  also  more  well-constrained,  since 
their  error  ellipses  are  smaller  than  those  derived  by  consi¬ 
dering  only  two  plates  at  a  time.  The  three  plate  method  for 
calculating  finite  rotations  should  prove  a  valuable  tool  for 

refining  our  ideas  of  the  tectonic  history  of  the  earth. 

*  *  * 

These  three  studies  demonstrate  the  value  of  considering 
the  evolution  of  three  plate  systems  and  the  stability  of 
triple  junctions  in  determining  and  refining  relative  rotations, 
both  instantaneous  and  finite.  These  methods  will  continue  to 
contribute  to  a  more  complete  and  detailed  understanding  of  the 
history  of  sea-floor  spreading.  Perhaps  the  most  interesting 
immediate  result  of  these  papers  is  the  surprising  constancy 
of  the  rotations  between  the  Indian  Ocean  plates  since  40  Ma, 
and  especially  since  10  Ma.  More  detailed  reconstructions  of 
the  region  should  now  be  possible  with  relatively  little  new 
data.  It  was  originally  intended  that  Chapter  4  would  illustrate 
the  three  plate  finite  rotation  method  with  an  application  to 
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the  Indian  Ocean  system.  However,  preliminary  investigations 
revealed  that,  except  for  the  area  near  the  triple  junction, 
well-located  crossings  of  magnetic  anomaly  lineations  and 
fracture  zones  are  still  too  sparse.  A  meaningful  treatment 
would  first  require  extensive  reinterpretation  of  published 
data  along  the  Central  and  Southeast  Indian  Ridges,  and  this 
is  a  subject  for  future  work. 
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ABSTRACT 

The  Southwestern  Indian  Ridge,  the  contact  between  the 
African  and  Antarctic  plates,  lies  between  the  Bouvet  Triple 
Junction  in  the  South  Atlantic  and  the  Indian  Ocean  Triple 
Junction  about  2100  kins  east  of  Madagascar.  From  the 
vicinity  of  Prince  Edward  Island  at  40°E  it  trends 
northeasterly  and  it  is  segmented  by  a  suite  of  deep  north- 
south  gashes  terminating  on  the  northeast  with  two  spectacular 
meridional  fracture  zones,  the  "Atlantis  II"  and  the  "Melville", 
at  57°30'E  and  60°30'E  respectively.  From  there  northeast  to 
the  Indian  Ocean  Triple  Junction  at  25°30'S,  70°00'E  the 
ridge  trends  N75°E;  it  is  characterized  by  a  triangle  of 
rough  topography  with  the  triple  junction  at  the  eastern  apex. 
From  all  available  data  an  instantaneous  pole  of  relative 
motion  for  Africa/Antarctica  was  computed;  it  lies  at  8.4°N, 

42 . 4°W,  with  a  rate  of  0.15  degrees/my. 

Since  the  marked  change  in  the  direction  and  rate  of 
spreading  in  the  Madagascar,  Crozet,  and  Central  Indian 
Basins  that  occurred  in  the  Eocene  (44  Ma,  Anomaly  19) ,  the 
poles  of  relative  motion  for  the  African,  Indian,  and 
Antarctic  plates  have  changed  very  little.  We  fixed  the 
Africa/Antarctica  and  Africa/India  poles  and  computed  that 
for  India/Antarctica.  We  justified  this  pole  by  comparisons 
of  predicted  isochrons  with  observed  magnetic  lineations  and 
determined  the  tectonic  history  of  the  triple  junction. 
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Since  the  Eocene  (44  Ma,  Anomaly  19),  this  junction  has 
moved  as  rapidly  eastwards  with  respect  to  Africa  as 
Antarctica  has  moved  south.  The  resultant  geometry  and 
slow  spreading  account  for  the  triangle  of  rough  topography 
produced  by  the  Southwest  Indian  Ridge  east  of  the  Melville 
Fracture  Zone.  The  triple  junction  evolved  as  a  stable 
Ridge-Ridge-Ridge  type  with  the  Southeast  Indian  Ridge 
remaining  approximately  constant  in  length.  It  was  not 
resolved  whether  this  constancy  in  length  is  maintained  by 
frequent  ridge  jumps  or  by  oblique  spreading  on  the 
Southwest  and  Central  Indian  Ridges  near  the  triple  junction. 


30 


INTRODUCTION 

The  topography  of  the  Indian  Ocean  is  dominated  by  three 
major  active  mid-ocean  ridge  systems.  The  most  northerly  of 
these  features  starts  in  the  Gulf  of  Aden,  becomes  the 
southeast-trending  Carlsberg  Ridge  southeast  of  the  gulf, 
and  at  1°N  and  60°E  continues  as  the  en  echelon  offset, 
northerly-trending  Central  Indian  Ridge.  Near  25°S,  at 
the  so-called  Indian  Ocean  Triple  Junction,  this  ridge 
bifurcates  into  the  Southeast  Indian  Ridge  and  Southwest 
Indian  Ridge.  The  Southeast  Indian  Ridge,  with  fast 
spreading  (^3  cm/yr)  and  relatively  smooth  topography, 
is  offset  strikingly  to  the  south  by  two  fracture  zones 
close  to  Amsterdam  Island  and  St.  Paul  Island  (Schlich  and 
Patriat,  1971)  before  trending  southeastward  between  Australia 
and  Antarctica.  The  Southwest  Indian  Ridge,  by  contrast, 
spreads  more  slowly  (<1  cm/yr)  and  it  has  much  rougher 
topography  and  many  more  extremely  deep  fracture  zones.  This 
ridge  system  is  offset  markedly  to  the  south  in  the  vicinity 
of  Marion  and  Prince  Edward  Islands  and  it  terminates  at  the 
Bouvet  Triple  Junction  in  the  South  Atlantic. 

The  three  major  active  mid-ocean  ridges  mark  the  present 
boundaries  of  the  Indian,  African,  and  Antarctic  plates.  Each 
of  these  active  plate  boundaries  has  been  studied  in  detail 
in  areas  of  considerable  extent  along  its  length.  For 
example,  McKenzie  and  Sclater  (1971)  and  Fisher  and  others 
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(1971)  have  elucidated  the  present  tectonic  structure  of  the 
Carlsberg  and  Central  Indian  Ridges,  and  Weissel  and  Hayes 
(1972)  have  sorted  out  the  history  of  separation  of 
Australia  and  Antarctica  from  initial  break-up  to  the  present. 
The  southwest  branch  has  been  studied  south  of  Africa  by 
Norton  (1976)  and  Sclater  and  others  (1978)  and  Bergh  and 
Norton  (1976)  have  examined  the  northward  extension  of  the 
Prince  Edward  Fracture  Zone.  Apart  from  the  study  of  Fisher 
and  others  (1971)  on  the  Central  Indian  Ridge,  however, 
none  of  the  detailed  work  has  been  concentrated  close  to  the 
Indian  Ocean  Triple  Junction,  and  hence  little  is  known  about 
the  tectonic  history  of  this  feature.  Furthermore,  until 
recently  the  geometry  of  the  Southwest  Indian  Ridge  was 
poorly  known.  Poles  determined  for  the  Africa/Antarctica 
plate  boundary  near  Bouvet  Island  (Sclater  and  others,  1978) 
and  by  closure  around  the  Bouvet  Triple  Junction  (Forsyth, 
in  Sclater  and  others,  1976)  differ  strikingly  from  those 
determined  by  closure  around  the  Indian  Ocean  Triple  Junction 
(McKenzie  and  Sclater,  1971)  and  from  global  closure  (Minster 
and  others,  1974) .  This  has  led  to  the  tentative  suggestion 
that  the  African  plate  might  consist  of  two  subplates  (Forsyth, 
1976;  Sclater  and  others,  1976). 

This  lack  of  detailed  knowledge  of  the  late  Tertiary 
history  of  the  Indian  Ocean  contrasts  markedly  with  a 
surprisingly  complete  understanding  of  the  late  Cretaceous 
and  early  Tertiary  history.  For  example,  well  identified 
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and  clearly  lineated  anomalies  have  been  identified  in  the 
Central  Indian  Basin  by  McKenzie  and  Sclater  (1971) ,  in  the 
Wharton  Basin  by  Sclater  and  Fisher  (1974),  and  in  the 
Mascarene,  Madagascar,  and  Crozet  basins  by  Schlich  (1975) . 
These  authors  have  demonstrated  that  from  80  million  years 
ago  (Anomaly  34)  to  50  million  years  ago  (Anomaly  21)  the 
spreading  centers  trended  almost  east-west.  In  the  early 
Eocene,  India  collided  with  Eurasia  (Molnar  and  Tapponier, 
1976),  there  was  an  abrupt  slowdown  in  spreading,  and 
between  anomalies  21  and  16  (39  million  years  ago)  the 
direction  of  spreading  changed  to  northeast-southwest. 
Finally,  the  Chagos  Fracture  Zone  terminated  as  a  north- 
south  lineation,  the  Chagos-Maldives  region  split  off 
from  the  Mascarene  Plateau  and  the  en  echelon  Central  Indian 
Ridge  was  formed  (Fisher  and  others,  1971) .  However,  aside 
from  a  study  at  the  southwestern  end  of  the  Ninetyeast  Ridge 
by  Sclater  and  others  (1976),  almost  nothing  is  known  about 
the  history  of  the  Central  Indian  Ridge  and  Southeast 
Indian  Ridge  from  then  until  9  million  years  ago  (Anomaly  5) . 

Norton  and  Sclater  (in  preparation)  have  used  the  fit 
of  the  Chagos  Bank  into  Nazareth  Bank  in  the  Mascarene 
Plateau,  DSDP  Site  238  (Fisher  and  others,  1974)  and  a  recent 
reanalysis  of  the  magnetic  data  south  of  Australia  by 
Weissel  (1977)  to  reconstruct  the  positions  of  the  African, 
Antarctic,  and  Indian  Plates  39  Ma  (Anomaly  16) .  The 
striking  feature  of  this  reconstruction  is  that  the 
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positions  of  the  finite  poles  describing  the  motion  of  the 
three  plates  are  little  different  from  the  instantaneous 
poles.  These  authors  have  speculated  that  all  of  the 
changes  in  direction  of  the  three  plates  took  place  prior 
to  39  Ma  (Anomaly  16)  and  that  since  then  their  directions 
of  relative  motion  have  changed  very  little. 

In  this  paper  a  topographic  and  a  magnetic  anomaly 
chart  of  the  Southwest  Indian  Ridge  east  of  53°E  and  of 
the  Indian  Ocean  Triple  Junction  are  presented.  In  the 
analysis  of  these  data  the  following  problems  are 
addressed: 

(a)  What  is  the  geometry  of  the  Southwest  Indian 
Ridge? 

(b)  Is  the  direction  of  motion  observed  on  the  Africa/ 
Antarctic  plate  boundary  close  to  the  Indian  Ocean  Triple 
Junction  compatible  with  that  close  to  Bouvet  Island? 

(c)  What  has  been  the  tectonic  history  of  the  Southwest 
Indian  Ridge  and  the  triple  junction  east  of  53°? 

TOPOGRAPHY 

Only  since  the  International  Geophysical  Year  1957-58 
and  the  International  Indian  Ocean  Expedition  1960-65  has 
there  been  any  attempt  to  gather  substantial  quantities  of 
well-controlled,  precise  bathymetric  data  in  the  Indian 
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Ocean.  Most  of  the  topographic  information  on  the  Southwest 
Indian  Ridge  east  of  40°E  has  come  from  recent  cruises  of 
the  Scripps  Institution  of  Oceanography,  the  Woods  Hole 
Oceanographic  Institution,  and  the  Institut  de  Physique  du 
Globe  de  Paris.  The  first  detailed  attempts  to  understand 
the  tectonics  and  crustal  composition  were  carried  out  from 
R/V  Argo  by  Scripps  Institution  of  Oceanography  on  CIRCE 
Expedition  (1968) .  This  cruise  was  followed  by  more 
detailed  work  in  this  area  and  on  the  Central  Indian  Ridge 
from  R/V  Melville  during  S.I.O.'s  expedition  ANTIPODE 
1970-71.  The  final  work  reported  here  was  two  months  of 
R/V  Atlantis  II  Cruise  93,  legs  5  and  6,  to  the  Southwest 
Indian  Ridge  and  the  Indian  Ocean  Triple  Junction  in  1976 
and  two  lines  of  Mar ion  Duf resne  Cruise  11  run  in  the 
same  year.  The  precisely  controlled  topographic 
information  gathered  on  the  R/V  Melville  and  R/V  Atlantis  II 
cruises  marks  the  main  data  source;  it  is  complemented  by 
•  the  long,  generally  northwest-southeast  lines  run  by  the 

Institut  de  Physique  du  Globe  between  Reunion  and  Crozet, 

Reunion  and  Kerguelen,  and  Kerguelen  and  Amsterdam  (Gallieni 
prior  to  1971  and  Marion  Duf resne  since  then) .  A  few  other 
lines,  such  as  from  S.I.O.  expeditions  MONSOON  (1960), 

LUSIAD  (1962-63)  and  DODO  (1964),  and  some  isolated  lines 
from  reconnaissance  of  the  Lamont-Doherty  Geological 
Observatory  or  Glomar  Challenger  (in  1972)  also  have  been 
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used . 

The  most  spectacular  features  on  this  topographic 
chart  of  the  Southwest  Indian  Ridge  (Figure  1)  are  the 
three  deep  north-south  trending  gashes,  with  shallow  lips, 
between  57°  and  61°E  and  the  rough  v-shaped  axial 
topography  between  the  easterly  fracture  zone  and  the 
Indian  Ocean  Triple  Junction  near  25°30'S,  70°00'E.  The 
three  fracture  zones  were  recognized  and  delineated  aboard 
R/V  Melville  (Engel  and  Fisher,  1975)  and  R/V  Atlantis  II 
(this  paper) ;  the  two  most  prominent,  intersecting  the 
ridge  crest  at  60°30'E  and  57°00'E,  have  been  named  after 
these  two  research  vessels,  respectively.  Two,  at  least, 
reach  depths  of  more  than  6000  meters;  all  three  provide 
avenues  for  the  passage  of  deep  Antarctic  bottom  water  from 
the  Crozet  Basin  north  into  the  Mascarene  Basin  (Warren, 
1978) .  The  gashes  extend  well  to  the  north  and  south  of 
the  deepest  sections  which  characteristically  occur  at 
or  near  their  intersections  with  the  ridge  crest.  They 
all  offset  the  ridge  crest  as  it  is  identified  by 
earthquake  epicenters,  and  they  all  appear  to  have  some 
7  to  9  degrees  of  north-south  extent.  The  north-south 


cross-grain  of  the  Southwest  Indian  Ridge  terminates 
sharply  at  the  Melville  Fracture  Zone  at  61°E  (Figure  1). 
To  the  east  the  rough  topography  appears  to  form  a  right 
triangle  with  the  most  acute  apex  at  the  triple  junction 
and  the  shortest  side  being  the  Melville  Fracture  Zone. 
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FIGURE  1 

Topographic  chart  of  the  Southwest  Indian  Ridge 
between  53°E  and  the  Indian  Ocean  triple  junction. 


Within  this  triangle  there  is  a  spectacular  deep  that  extends 
overall  N75°E  in  a  sinuous  trend  from  that  fracture  zone  to 
the  triple  junction,  which  is  marked  by  an  isolated  pocket  or 
median  valley  greater  than  5000  m  deep  (Tapscott  and  others, 
in  preparation) . 

The  rough  topography  on  the  Southwest  Indian  Ridge 
contrasts  markedly  with  the  smoother,  much  blockier  relief 
of  the  Southeast  Indian  Ridge.  On  the  latter  the  contours 
are  dominated  by  N45°E  trending  fracture  zones  and  a  200  to 
300  m  deep  that  marks  the  axis  at  the  crest  of  the  ridge. 

The  Central  Indian  Ridge  is  also  characterized  by  a  deep  at 
the  axis  but  the  topography  is  more  rugged  and  the  ridge  axis 
is  more  segmented  by  huge  fracture  zones  that  expose  rocks  of 
the  lower  crust  (Engel  and  Fisher,  1975) .  This  is  especially 
true  to  the  north  of  the  area  shown  in  Figure  1,  notably 
for  the  Marie  Celeste  Fracture  Zone  at  18 °S. 

Away  from  the  active  ridge  sections  there  are  other 
topographic  provinces  of  importance .  East  and  southeast  of 
Mauritius  for  several  hundred  kilometers  the  topographic 
relief  is  extreme,  with  elongated  ridges  and  faulted  slivers 
separated  by  flat  floored  basins  of  turbidites  and  volcanic 
ash.  This  ridge  and  pond  topography  trends  nearly  north- 
south  in  a  large  region  west  of  64 °E,  except  for  the  sharply 
defined  Rodriguez  Ridge  at  about  20 °S  and  a  discordant  knot 
near  23°S,  59°30'E.  This  area  contrasts  markedly  with  the 
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smooth  volcanogenic  apron  surrounding  and  south  of  Reunion, 
and  also  with  the  sedimented  and  smoothed  tract  from 
Reunion  south  to  the  northeast-trending  northern  foothills 
of  the  Southwest  Indian  Ridge  at  52°-59°E.  South  and 
southwest  of  the  triple  junction  for  a  considerable  distance 
are  northeast-southwest  trending  ridges  and  swales  whose 
origin  is  obscure,  but  which  also  occur  northeast  of  the 
Southeast  Indian  Ridge  as  pictured  on  Figure  1. 

MAGNETIC  ANOMALIES  AND  SEISMICITY 

Figure  2  portrays  the  residual  magnetic  anomalies  (total 
field  with  IGRF  removed  for  data  prior  to  1971;  total  field 
with  IGRF  +  400  gammas  removed  for  data  since  then)  along 
trace  for  all  the  available  Scripps  Institution  of  Oceanography, 
Woods  Hole  Oceanographic  Institution,  and  Institut  de  Physique 
du  Globe  cruises.  A  few  reconnaissance  lines  from  CONRAD  16 
and  VEMA  20  of  the  Lamont-Doherty  Geological  Observatory  are 
also  shown.  Anomalies  in  the  area  around  the  triple  junction 
are  not  drawn.  These  data  are  shown  at  a  larger  scale 
elsewhere  (Tapscott  and  others,  in  preparation).  Also 
superimposed  upon  the  lineations  are  the  relocated  earthquake 
epicenters  from  the  ESSA  earthquake  tape  between  1971  and 
1976.  These  earthquakes  were  relocated  by  Sean  Solomon  of 
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FIGURE  2 

Profiles  of  residual  magnetic  anomalies  at  right 
angles  to  track  in  the  southwest  Indian  Ocean.  All 
cruises  after  1971  have  400  gammas  added  to  the  IGRF 
(IAGA,  1965) .  The  magnetic  anomalies  in  the  area  of 
the  rectangle  have  been  interpreted  by  Tapscott  and 
others  (in  preparation) .  The  circled  crosses  represent 
relocated  earthquake  epicenters,  the  relative  size 
indicating  magnitude.  The  black  dots  with  lines  through 
them  represent  magnetic  anomalies.  The  heavy  lines 
mark  the  active  spreading  centers  and  transform  faults 
and  the  light  dashed  lines  mark  fracture  zone  traces. 
The  hatched  line  marks  the  rough-smooth  boundary  around 
the  Southwest  Indian  Ridge. 
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On  superimposing  the  magnetic  anomaly  and  seismicity 
chart  upon  the  topographic  contours ,  in  general  the  earthquake 
epicenters  and  the  pronounced  positive  anomalies  lie  close  to 
the  elevated  topography.  When  examined  in  detail  both  the 
center  of  the  Central  Indian  and  of  the  Southeast  Indian 
Ridges  are  associated  with  a  500  m  depression.  The  axis  of 
the  Southwest  Indian  Ridge  is  marked  by  a  much  more  pronounced 
cleft  that  generally  is  deeper  than  1000  m  and  in  some 
localities  extends  to  2000  m.  In  the  case  of  the  Central 
and  Southeast  Indian  Ridges  these  depressions  correlate  well 
with  the  center  of  the  central  magnetic  anomaly  and  hence  they 
are  interpreted  to  mark  the  axis  of  spreading.  In  general, 
the  offsets  on  the  central  anomalies  in  these  two  regions  are 
small  and  few  very  deep  transform  faults  are  observed.  However 
on  the  Central  Indian  Ridge  at  20°S  and  23°S  and  on  the 
Southeast  Indian  Ridge  at  26°S,  27°30'S  and  28°30'S  small 
but  observable  transform  faults  are  revealed  both  by 
topography  and  by  offsets  of  the  magnetic  lineations. 

In  contrast  to  the  other  ridge  axis  sections  in  this 
region,  the  Atlantis  II  and  Melville  Fracture  Zones  are 
spectacular  features.  They  show  a  pronounced  north-south 
trend  in  both  the  topography  and  the  relocated  epicenters. 

There  probably  are  two  north-south  transform  faults  between 
them  and  together  these  four  features  offset  the  active 


43 


spreading  centers  5°  northward  in  an  east-west  distance  of 
4°.  All  the  epicenters  lie  either  on  transform  faults  or 
spreading  centers . 

East  of  61°E  (Melville  Fracture  Zone)  the  morphology 
and  activity  of  the  ridge  trend  plainly  N75°E  overall. 

Though  the  general  trend  of  the  topography  is  northeast,  the 
ridge  is  characteristically  blocky,  with  generally  east-west 
troughs  joined  and  offset  by  broader  north-south  depressions. 
Between  66°E  and  67°E  there  is  a  clear  central  magnetic 
anomaly  associated  with  a  deep  depression  in  the  topography. 
It  is  assumed  that  except  for  the  area  between  61°S  and  63°E 
this  depression  marks  the  active  spreading  region,  here  a 
zone  consisting  of  east-west  trending  spreading  centers 
offset  by  short  north-south  transform  faults.  However,  to 
be  consistent  one  is  constrained  by  the  data  control  to 
assume  that  between  61°  and  63°E  the  ridge  axis  likewise 
follows  the  central  depression  and  hence  runs  about  30° 
oblique  to  the  normal  or  at  60°  to  the  direction  of  the 
transform  faults. 

Older  Anomalies 

The  magnetic  anomalies  on  the  Southwest  Indian  Ridge 
are  difficult  to  interpret  owing  to  the  slow  spreading 
rate  (< 1  cm/yr)  and  the  associated  rough  topography.  Further 
complications  arise  because  the  axis  consists  of  short 
sections  of  spreading  ridge  offset  by  fracture  zones. 
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Thus  it  is  difficult  to  obtain  a  magnetic  profile  that  has 
not  crossed  a  fracture  zone.  The  ridge  crest  profiles  that 
were  run  north-south  on  ATLANTIS  II  (All)  93-5  and  one 
carefully  controlled  line,  MARION  DUFRESNE  (MD)  11-2,  are 
presented  as  Figure  3.  The  latter  line  was  run  between  two 
clearly  identified  central  anomalies  recorded  by  Atlantis  II. 
An  attempt  was  made  to  match  the  lineations  on  the  ridge  axis 
profiles  with  the  standard  block  model  and  a  'contaminated' 
model.  In  the  'contaminated'  model  the  box-like  feature  of 
blocks  has  been  smoothed  by  a  Gaussian  filter  of  variable 
length  (Tisseau,  1978).  The  Gaussian  filter  models  the 
random  nature  of  the  intrusion  process  at  the  spreading 
center.  The  exact  nature  of  the  computation  process  can  be 
found  in  Tapscott  and  others  (in  preparation) .  With  the  aid 
of  the  contamination  model  one  is  able  to  observe  Anomaly  5 
on  both  ends  of  the  MD-11  profile,  and  on  one  or  perhaps 
two  of  the  All  93-5  profiles.  We  were  not  able  to  identify 
with  certainty  any  anomalies  older  than  Anomaly  5. 

The  magnetic  anomalies  on  the  Central  Indian  Ridge  are 
much  easier  to  identify  than  those  on  the  Southwest  Indian 
Ridge  because  its  spreading  rate  is  faster  and  the  topography 
less  rough.  A  previous  compilation  by  Fisher  and  others 
(1971)  identified  many  anomalies,  from  the  central  anomaly 
to  Anomaly  5  either  side  of  the  ridge  axis.  Here  we  present 
a  more  extensive  compilation  of  the  data  but,  except  for 


FIGURE  3 


Residual  magnetic  anomaly  profiles  across  the  South¬ 
west  Indian  Ridge  (for  position  see  numbers  on  profiles 
on  Figure  2)  projected  onto  000°  and  compared  with 
synthetics,  skewed  60°,  unaltered  and  contaminated  by 
a  4  km  Gaussian  filter.  Spreading  rate  0.8  cm/yr. 
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the  detailed  survey  near  the  triple  junction,  few  new 
identifications.  However  one  is  able  to  recognize  with 
more  certainty  than  Fisher  and  others  (1971)  the  transform 
faults  and  offsets  on  the  lineations.  In  the  present  paper 
the  major  advance  has  been  on  anomalies  older  than  Anomaly  5. 
All  93-5,  All  93-6,  and  MD-11  profiles  were  run  closely 
parallel  to  what  was  thought  to  be  the  direction  of 
spreading.  These  lines  were  most  successful,  and  on  three 
profiles  we  identified  with  certainty  anomalies  6,  8,  13 
and  possibly  Anomaly  18  on  two  lines  and  on  one,  MD-11, 
anomalies  19  and  20  that  obviously  represent  spreading  at 
a  much  faster  rate  than  the  younger  anomalies  (Figure  4). 
Anomaly  6  also  is  recognizable  on  a  VEMA-20  profile  to  the 
east  of  the  ridge  axis.  From  the  presumed  lineation  of  these 
anomalies  and  the  apparent  rate  it  is  clear  that  the  present 
phase  of  spreading  extends  back  to  at  least  the  time  of 
Anomaly  18  and  probably  until  about  Anomaly  19.  At  that 
time  there  must  have  been  a  major  change  in  direction, 
following  the  formation  of  the  northwest-southeast  trending 
anomalies  in  the  Madagascar  Basin. 

Of  the  three  spreading  centers,  the  Southeast  Indian 
Ridge  has  by  far  the  easiest  anomalies  to  identify.  Though 
there  are  considerably  fewer  tracks  on  this  section  of  the 
ridge  the  trends  and  rates  of  spreading  are  very  easily 
determined  from  the  magnetic  lineations.  As  all  the 
profiles  over  the  Southeast  Indian  Ridge  shown  here  have 


FIGURE  4 


Residual  magnetic  anomaly  profiles  across  the  Central 
Indian  Ridge,  projected  onto  N060°E.  These  are 
compared  with  synthetics,  skewed  30°,  unaltered  and 
contaminated  by  a  4  km  Gaussian  filter.  Spreading  rate 
2.0  (anomalies  8  to  5)  and  1.8  (anomalies  21  to  8)  cm/yr. 
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been  discussed  elsewhere  we  have  not  presented  a  figure  of 
the  projected  anomalies.  North  of  the  ridge  axis  and  just 
west  of  the  Ninetyeast  Ridge,  Sclater  and  others  (1976) 
reported  a  suite  of  anomalies  6  through  17  with  one  clear 
Anomaly  19  at  25°30'S  and  85°E.  To  the  south  of  the 
Southeast  Indian  Ridge  anomalies  6,  8,  13  and  19  are  all 
clearly  identified  on  profiles  presented  by  Schlich  (1975). 
Also  one  may  recognize  in  the  topography  possible  fracture 
zones  trending  northeast-southwest  at  right  angles  to  the 
proposed  magnetic  lineations.  The  magnetic  lineations  on 
the  southeast  branch  reveal  three  main  features.  First, 
though  the  ridge  has  a  relatively  complicated  blocky 
morphology,  the  ridge  axis  sections  are  much  longer  than 
the  transform  faults,  and  in  general  there  is  a  relatively 
simple  evolutionary  history.  Second,  on  the  profiles  both 
north  and  south  of  the  ridge  the  direction  of  spreading  and 
the  rate  of  spreading  changed  at  Anomaly  19  time.  Finally, 
the  identified  anomalies  7  through  13  north  of  the  ridge 
axis  are  about  1  1/2  degrees  of  latitude  closer  to  the 
ridge  axis  than  are  those  to  the  south.  Preliminary 
comparison  of  profiles  suggests  that  this  is  due  to  a 
general  asymmetry  by  a  jump  or  continuous  asymmetric 
spreading  between  anomalies  5  and  6. 


51 


TECTONIC  CHART 

Figure  5  is  a  tectonic  chart  of  the  area  between  15 °S 
and  35°S,  and  between  54 °E  and  85 °E,  based  upon  the  topographic 
and  magnetic  data  presented  in  Figures  1  and  2  and  that 
presented  by  Tapscott  and  others  (in  preparation)  at  the 
triple  junction.  For  the  areas  outside  these  figures  we 
complemented  our  data  by  accepting  the  interpretation  of 
Schlich  (1975)  for  the  Mascarene  and  Crozet  Basins,  Fisher 
and  others  (1971)  for  the  Central  Indian  Ridge,  and  Sclater 
and  others  (1976)  and  Sclater  and  Fisher  (1974)  for  the 
Central  Indian  Ridge. 

West  of  the  Melville  Fracture  Zone  at  61°E  the 
Southwest  Indian  Ridge  consists  of  a  series  of  diffuse  and 
not  clearly  identified  spreading  centers  offset  in  significant 
degree  by  an  en  echelon  suite  of  north-south  trending  transform 
faults.  East  of  the  Melville  Fracture  Zone  the  ridge  becomes 
a  tenuous  set  of  east-west  trending  spreading  centers  offset 
by  short  and  difficult  to  identify  north-south  transform 
faults.  The  center  of  spreading  is  remarkably  well  delineated 
by  the  4000  m  contour  (Figure  5) . 

The  north-south  trending  fracture  zones  are  spectacular 
in  both  depth  and  extent  and  their  topographic  expression 
can  be  extended  to  just  south  of  the  well-identified  magnetic 
lineations  in  the  Madagascar  Basin.  To  the  east  of  the 
Melville  Fracture  Zone  this  pronounced  north-south  grain 


FIGURE  5 


A  tectonic  chart  of  the  area  surrounding  the  Indian 
Ocean  triple  junction  superimposed  upon  the  4000  m 
contour.  The  continuous  or  dashed  lines  represent 
the  central  anomalies.  The  single  dashed  lines 
represent  active  transform  faults  and  fracture  zones. 
The  hatched  line  marks  the  boundary  of  the  Southwest 
Indian  Ridge.  The  black  dots  with  lines  through  them 
are  magnetic  anomalies.  The  circled  crosses  are  recent 
epicenters,  relative  size  indicating  magnitude. 
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disappears.  The  trends  well  south  of  the  Southwest  Indian 
Ridge  spreading  center  appear  northeast-southwest  in  direction. 

The  rough  to  smooth  topographic  transitions  north  and 
south  of  the  axes,  that  mark  the  presumed  boundaries  between 
the  Southwest  Indian  Ridge  and  other  ridges,  have  a  complex 
character.  Between  the  triple  junction  and  the  Melville 
Fracture  Zone  at  61°E  the  northern  boundary  runs  east-west 
but  is  not  very  clear,  whilst  the  southerly  boundary  runs 
southwest-northeast  and  is  well-defined.  The  area  enclosed 
by  those  boundaries  is  a  right  triangle,  with  the  most  acute 
apex  at  the  triple  junction.  West  of  the  Melville  Fracture 
Zone  the  rough/smooth  boundary  between  the  Southwest  Indian 
Ridge  and  the  other  two  ridges  trends  at  a  much  steeper  angle 
both  north  and  south  of  the  ridge  than  to  the  east  of  61°E. 

The  Central  Indian  Ridge,  north  of  the  Marie  Celeste 
Fracture  Zone  at  18°S,  is  marked  by  very  deep  en  echelon 
fracture  zones  with  clearly  identified  and  massively  offset 
ridge  segments  between  them.  South  of  this  fracture  zone 
the  ridge  segments  appear  longer  and  the  transform  faults 
are  less  obvious  on  the  topography.  The  magnetic  anomalies 
can  be  traced  with  confidence  to  Anomaly  8  and  with  less 
confidence  to  Anomaly  19.  However,  between  Anomaly  19 
and  Anomaly  21  time  there  is  a  striking  change,  since  at 
50  Ma  (Anomaly  21)  the  spreading  direction  apparently  was 


N35°E . 


The  Southeast  Indian  Ridge,  in  contrast,  has  the 
characteristic  blocky  morphology  of  relatively  fast  (^3  cm/yr) 
spreading  centers  though  it  still  possesses  a  central  valley. 
Most  of  the  recent  epicenters  lie  on  transform  faults  and 
few  lie  on  the  active  spreading  center  sections.  A  surprisingly 
large  grouping  of  earthquakes  is  observed  at  32°S  and  79°E. 

These  earthquakes  apparently  lie  on  an  inactive  section  of  a 
transform  fault.  Their  origin  and  the  explanation  of  the 
large  number  is  unknown  at  this  time.  From  the  present  back 
until  approximately  44  Ma  (Anomaly  19)  the  direction  of 
spreading  has  been  roughly  constant  at  N45°E  and  the  rate 
has  remained  between  2.5  and  3  cm/yr  with  some  asymmetry. 
However,  earlier  than  this,  between  50  Ma  (Anomaly  21)  and 
80  Ma  (Anomaly  34)  the  spreading  direction  was  very  different, 
N05°E  in  the  Central  Indian  Basin  (McKenzie  and  Sclater,  1971) , 
and  N30°E  in  the  Crozet  Basin  (Schlich,  1975) ,  the  apparent 
discrepancy  in  both  trends  being  due  to  the  curvature  of  the 
earth. 

PRESENT  PLATE  MOTION:  AFRICA/ANTARCTICA 

One  of  the  main  objectives  of  this  study  was  to  determine 
an  instantaneous  pole  of  relative  motion  for  Africa/Antarctica. 
To  obtain  this  pole,  we  determined  the  positions  along  the 
plate  boundary  where  we  believed  the  directions  of  motion 
were  well-known  (Table  1;  Fig.  6) .  These  positions  and 


TABLE  1 


Directions  of  Spreading  on  the  Southwest  Indian  Ridge 


Feature 

Latitude 

Longitude 

Triple  Junction 
(1) 

25.60°S 

70.00°E 

Spreading  Center 
(2) 

27. 30°S 

66. 63°E 

Melville  F.Z. 

(3) 

30.00°S 

60. 75°E 

"588E  F.Z." 

(4) 

31.70*8 

58. 35°E 

Atlantis  II  F.Z. 

(5) 

33.00°S 

57.00°E 

Prince  Edward  F.Z. 
(6) 

46.00°S 

35.15°E 

Shaka  F.Z. 

(7) 

53.50°S 

09.00°E 

I si as  Orcadas  F.Z. 
(8) 

54.25°S 

06.00°E 

Bouvet  F.Z. 

(9) 

54.258S 

02.008E 

Direction 

Observed  Predicted 

Date  Source 

01°±09® 

358° 

Tapscott  and  others 
(in  prep.) 

358°±05° 

359° 

This  paper 

3548±048 

O 

rH 

O 

This  paper 

028±05° 

o 

CM 
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FIGURE  6 


A  comparison  between  predicted  isochrons  and  observed 
magnetic  lineations  in  the  vicinity  of  the  Indian 
Ocean  Triple  Junction.  The  isochrons  were  computed 
by  rotating  the  present  ridge  axis  by  half  the  angle 
about  the  poles  given  in  Table  2b.  The  dotted  line 
represents  the  Melville  Fracture  Zone. 


directions  were  combined  two  at  a  time  to  determine  the  pole 
of  a  rotation  which  fit  both  points.  Only  positions  separated 
by  more  than  4°  in  latitude  were  selected,  and  these  were 
weighted  according  to  the  accuracy  with  which  it  was  believed 
the  direction  was  known.  These  statistics  have  a  weighted 
mean  of  the  individual  poles  of  8.46N,  42.2°W,  with  an  a95 
of  4.0°  and  a  rate  of  0.15  degrees/my.  A  comparison  of  the 
directions  predicted  by  this  pole  and  those  observed  at  the 
nine  locations  selected  shows  that  none  disagree  by  more  than 
02°.  This  fit  is  considered  excellent  and  strikingly 
demonstrates  that,  at  the  level  of  accuracy  at  which  we  are 
operating,  Africa  appears  to  be  a  single  plate  and  the 
Bouvet  Triple  Junction  data  is  compatible  with  that  from  the 
Indian  Ocean  Triple  Junction. 

It  is  of  interest  to  note  that  this  mean  pole  lies  more 
than  30°  from  that  predicted  by  Minster  and  others  (1974)  for 
their  global  closure  model.  It  also  lies  15°  away  from  the 
pole  calculated  by  Forsyth  (reported  in  Sclater  and  others, 
1976)  by  closure  around  the  Bouvet  Triple  Junction.  However, 
it  is  remarkably  close  to  that  predicted  by  Chase  (1972) . 
Clearly  some  reanalysis  of  the  closure  models  of  Minster  and 
others  (1974)  is  needed  before  they  are  preferred  over  those 
of  Chase  (1972) . 

In  an  imminent  paper  Tapscott  and  others  (in  preparation) 
compiled  60  observations  of  rate  and  instantaneous  direction 
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of  relative  motion  for  the  South  American,  African,  Indian, 
and  Antarctic  plates.  They  added  to  this  set  some  observations 
from  the  Bouvet  and  Indian  Ocean  Triple  Junctions  and  used  a 
program  written  by  Forsyth  (1976)  to  find  the  best  fit  set  of 
relative  rotation  vectors.  They  find  an  Africa/Antarctica 
pole  at  5.6°N,  39.4°W  with  a  rate  of  0.15  degrees/my.  It  is 
gratifying  to  see  how  well  this  more  general  set  gave  almost 
the  same  pole  for  Africa/Antarctica  as  given  by  our  independent 
calculation.  Our  pole,  which  agrees  with  the  data  and  is 
consistent  with  those  on  four  other  plate  boundaries,  is  to  be 
preferred  to  that  of  Minster  and  others  (1974) . 


TECTONIC  HISTORY 

The  fracture  zones  on  the  Southwest  Indian  Ridge  trend 
north-south  and  extend  a  considerable  distance  in  this  direction 
to  the  north  of  the  active  section.  But  no  anomalies  older 
than  Anomaly  5  were  identified  on  this  ridge  axis.  Thus  there 
is  no  possibility  with  the  available  data  set  of  directly 
determining  whether  or  not  the  present  spreading  regime  has 
persisted  since  44  Ma  (Anomaly  19)  as  is  observed  on  the 
Central  Indian  and  Southeast  Indian  Ridges.  However,  we  can 
indirectly  resolve  the  spreading  history  of  the  Southwest 
Indian  ridge  by  reconstructing  the  past  position  of  the  other 
two  ridges.  This  is  possible  to  the  east  of  the  Melville 
Fracture  Zone  because  the  southernmost  Anomaly  19  on  the 
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Central  Indian  Ridge  lies  directly  north  of  the  most 
northern  Anomaly  19  in  the  Crozet  Basin  (Figure  5) .  Further, 
both  anomalies  lie  just  east  of  the  Melville  Fracture  Zone 
and  the  gap  between  them  encloses  almost  the  entire  region 
of  rough  topography  between  the  Melville  and  the  triple 
junction. 

Norton  and  Sclater  (in  press)  have  used  the  closure  of 
the  Chagos  Bank  onto  the  Mascarene  Plateau  and  alignment  of 
the  Vishnu  Fracture  Zone  and  Mauritius  Trench  to  compute  an 
Anomaly  16  (39  Ma)  pole  for  the  motion  of  Africa  with  respect 
to  Australia.  Weissel  and  others  (1977)  have  used  the  magnetic 
anomalies  between  Australia  and  Antarctica  to  compute  an 
Anomaly  18  pole  for  the  motion  of  India  with  respect  to 
Antarctica.  It  is  apparent  that  since  Anomaly  19  time  the 
finite  poles  must  have  lain  fairly  close  to  the  present  poles 
for  the  motions  to  have  remained  so  constant  both  in  rate  and 
in  direction  on  the  three  ridge  axes.  As  a  consequence,  one 
can  assume  these  Africa/India  and  India/Antarctica  poles  to 
hold  back  to  Anomaly  19  time  (44  Ma)  (Table  2a) .  One  then 
computes  a  pole  for  the  motion  of  Antarctica  with  respect 
to  Africa.  This  lies  close  to  the  present  pole.  Rotating 
the  present  Central  Indian  Ridge  axis  by  the  Anomaly  19 
pole  for  Africa/India  gives  a  good  fit  to  the  observed  Anomaly  19 
in  the  Madagascar  Basin.  But  rotating  the  Southeast  Indian 
Ridge  axis  by  the  Weissel  and  others  (1977)  pole  gives  a  very 
poor  match,  even  after  allowing  some  asymmetry. 


TABLE  2a 


Poles  Computed  From  Tapscott  and  others  (in  preparation), 
Weissel  and  others  (1977)  and  Norton  and  Se.later  (in  preparation) 


Antarctica/Africa 

S.W.I.R. 


Anomaly 

Age 

Lat. 

Long. 

Angle 

5 

9 

7.8 

-41.8 

1.38* 

19 

44 

3.1 

-38.3 

8.0* 

Africa/India  Antarctica/India 


C, 

Lat . 

.I.R. 

Long. 

Angle 

S 

Lat . 

■  E  •  I  •  R  • 
Long. 

Angle 

15.5 

48.6 

5.71 

16.2 

34.6 

5.91 

13.3 

54.1 

23. 52 

9.9 

34.8 

24. 63 

*Resultant  pole 

Pole  from  Tapscott  and  others  (in  preparation) 

Pole  computed  from  Norton  and  Sclater  (in  press) 

3Pole  computed  from  Weissel  and  others  (1977) 

TABLE  2b 

Poles  Used  for  Reconstructing  Position  of  Ridge  Axis 

Antarctica/Africa  Africa/India  Antarctica/India 


Anomaly 

Age 

Lat. 

S.W.I. 

Long. 

,R. 

Angle 

Lat. 

C.I.R. 

Long. 

Angle 

S 

Lat . 

.E.I.R. 

Long. 

Angle 

5 

9 

8.4 

-42.2 

1.41 

16.0 

48.3 

5.82 

16.8 

34.3 

6.0* 

8 

27 

8.4 

-42.2 

4 . 91 » 3 

16.0 

48.3 

15. I2’3 

15.4 

29.8 

16.0* 

19 

44 

8.4 

-42.2 

8.01, 3 

16.0 

48.3 

23. 22, 3 

14.0 

28.8 

24.8* 

*Resultant  pole 

3This  paper 
2 

Pole  from  Fisher  and  others  (1971) 

'l 

Angle  adjusted  to  superpose  ridge  axis  on  proper  anomaly 

In  both  tables  for  the  poles  of  rotation,  positive  number  indicate  northern 
latitudes  and  eastern  longitudes  and  a  right-handed  screw  about  an  axis  passing 
through  the  pole  and  the  center  of  the  earth. 
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The  present  data  set  are  not  good  enough  to  obtain 
improved  finite  poles  but  they  are  good  enough  to  demonstrate 
the  general  development  of  the  triple  junction.  Thus  in 
order  to  illustrate  how  the  Southwest  Indian  Ridge  has  been 
created  by  the  rapid  eastward  motion  of  the  triple  junction 
we  opted  for  the  simplest  approach  compatible  with  the  data. 
First,  we  assumed  that  the  present  poles  of  relative  motion 
between  Africa  and  Antarctica  and  between  Africa  and  India 
have  remained  fixed  since  44  Ma  (Anomaly  19) .  We  justify 
this  by  the  long,  straight  fracture  zone  traces  observed  on 
the  Southwest  Indian  Ridge  and  by  the  major  fracture  zones 
on  the  Central  Indian  Ridge  and  the  close  agreement  between 
the  poles  of  Tapscott  and  others  (in  preparation)  and  Norton 
and  Sclater  (in  press) .  For  the  27  Ma  (Anomaly  8)  and  44  Ma 
(Anomaly  19)  rotations  for  the  Central  Indian  Ridge,  we 
computed  the  angles  of  rotation  by  fitting  the  ridge  axis 
to  the  observed  anomalies  8  and  19.  The  angle  of  rotation 
for  9  Ma  (Anomaly  5)  on  the  Southwest  Indian  Ridge  was  obtained 
by  assuming  the  present  rate  of  spreading  to  have  been  constant 
from  that  time.  The  angle  for  the  44  Ma  (Anomaly  14) 
rotation  we  determined  by  fitting  the  ridge  axis  to  the 
rough-smooth  topographic  boundary  just  south  of  Anomaly  19 
on  the  Central  Indian  Ridge,  the  assumption  being  that 
relative  to  Africa,  44  million  years  ago  the  triple  junction 
lay  directly  north  of  the  Melville  fracture  zone  and  just 
south  of  the  observed  Anomaly  19  on  the  MD-11  profile. 
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For  the  27  Ma  (Anomaly  8)  rotation  angle,  we  assumed  a 
constant  rate  of  spreading  on  the  Southwest  Indian  Ridge 
between  44  Ma  and  9  Ma.  We  believe  this  process  does 
not  introduce  significant  errors  because  the  average  rate 
for  the  44  Ma  rotation  is  only  slightly  greater  than  that 
for  9  Ma. 

Using  these  six  rotations  we  computed  the  poles  and 
angles  of  relative  rotation  between  India  and  Antarctica 
for  9  Ma,  27  Ma,  and  44  Ma  (Table  2b) .  These  poles  were 
checked  by  rotating  the  present  ridge  axis  by  half  the  total 
angle  to  either  side  to  obtain  predicted  isochrons.  The 
match  to  Anomaly  5  is  excellent.  For  Anomaly  8  the  overall 
angle  is  too  small,  but  satisfactory,  but  spreading  has  been 
strongly  asymmetric  and  there  is  considerably  more  ocean  south 
of  the  ridge  than  north  of  it  between  anomalies  5  and  8.  For 
Anomaly  19  the  asymmetry  is  still  in  existence  but  the  overall 
angle,  given  the  simplicity  of  the  method  used  to  obtain  it, 
is  surprisingly  good. 

With  these  poles,  one  rotates  the  Southwest  Indian  and 
Central  Indian  Ridges,  as  well  as  the  Southeast  Indian  Ridge, 
about  the  present  ridges  by  half  the  rotation  angles  of  the 
respective  9  Ma,  27  Ma,  and  44  Ma  rotations.  Some  interesting 
points  for  the  development  of  the  Southwest  Indian  Ridge 
become  apparent  after  these  rotations.  For  example,  the 
most  northerly  Anomaly  8  in  the  Crozet  Basin  and  the  shape 
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of  the  rough-smooth  boundary  between  the  Southwest  Indian 
and  Southeast  Indian  Ridges  demonstrates  that  the  Southeast 
Indian  Ridge  has  shortened  only  very  little  since  44  Ma 
(Anomaly  19) .  Because  one  knows  with  some  precision  the 
amount  by  which  this  ridge  has  shortened  with  time,  one 
can  determine  the  amount  that  the  Central  Indian  and 
Southwest  Indian  Ridges  have  changed  in  length.  Both  the 
Central  Indian  and  Southwest  Indian  Ridges  have  increased 
in  length,  the  Southwest  Indian  Ridge  by  1000  kms  and  the 
Central  Indian  Ridge  by  550  kms.  The  increase  in  length  of 
the  Southwest  Indian  Ridge  since  44  Ma  (Anomaly  19)  is 
greater  than  the  800  kms  it  has  opened.  It  is  this  relatively 
rapid  eastward  motion  of  the  Indian  Ocean  Triple  Junction  and 
the  relatively  slow  north-south  spreading  on  the  Southwest 
Indian  Ridge  that  creates  the  triangle  of  rough  topography 
with  the  apex  at  the  triple  junction.  The  right-sided  western 
boundary  of  the  triangle  is  the  Melville  Fracture  Zone.  This 
suggests  most  strongly  that  this  phase  of  eastward  drift  of 
the  triple  junction  started  immediately  after  the  formation 
of  the  Melville  Fracture  Zone  and  the  abrupt  change  of  spreading 
on  the  Central  and  Southeast  Indian  Ridges  at  44  Ma. 

The  general  trend  of  the  Southwest  Indian  Ridge  is 
oblique  to  the  Melville  and  Atlantis  II  fracture  zones. 

The  trace  of  the  triple  junction  is  east-west  on  the  African 
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plate  and  northeast  on  the  Antarctic  plate.  West  of  the 
Melville  Fracture  Zone  the  trend  of  the  Southwest  Indian 
Ridge  is  even  more  oblique,  and  if  the  same  model  holds  it  is 
tempting  to  suggest  that  this  rough/smooth  boundary  observed 
between  the  Southwest  Indian  Ridge  and  the  Madagascar  Basin 
is  the  trace  of  the  triple  junction  on  the  African  plate 
between  anomalies  25  and  19. 

DEVELOPMENT  OF  THE  TRIPLE  JUNCTION 

There  are  considerable  uncertainties  in  both  the  lengths 
of  the  individual  ridge  axis  sections  chosen  and  the  poles 
used  for  the  rotations.  Thus  one  cannot  give  an  exact 
representation  of  the  evolution  of  the  triple  junction. 

On  the  other  hand,  we  believe  that  for  the  most  part  it  has 
had  a  stable  RRR  configuration.  It  is  easy  to  see  on  the 
tectonic  history  diagram  (Fig.  6)  that  at  9  Ma,  (Anomaly  5) , 

27  Ma  (Anomaly  8) ,  and  44  Ma  (Anomaly  19)  the  ridge  axes  all 
meet  at  a  point  without  the  necessity  of  transform  offsets. 

The  tectonic  diagram  has  been  used  to  construct  a  mean 
velocity  triangle  for  this  time  span  (Figure  7a) .  This 
triangle  enables  one  to  discuss  the  general  evolution  of  the 
triple  junction.  Assuming  these  mean  directions  and  symmetric 
spreading  at  right  angles  to  the  spreading  centers,  the  triple 
junction  (J,  Figure  7a)  on  the  velocity  diagram  should  move 
southeastward  down  the  Southeast  Indian  Ridge  as  time  passes. 
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FIGURE  7 

a.  A  velocity  diagram  of  the  Indian  Ocean  triple  junction 
during  phases  of  normal  RRR  evolution.  For  symmetric 
and  non-oblique  spreading,  the  velocity  of  the  junction 
is  given  by  J.  To  maintain  an  approximately  constant 
length  for  the  Southeast  Indian  Ridge,  the  triple 
junction  must  periodically  jump  to  the  northwest  along 
the  direction  of  that  ridge,  creating  transform  faults 
on  the  Southwest  and  Central  Indian  Ridges  (light 
dashed  lines) . 

b.  A  schematic  representation  of  the  trace  of  migration 
of  the  triple  junction  on  the  three  plates,  assuming 
frequent  ridge  jumps. 

c.  A  velocity  diagram  of  the  Indian  Ocean  triple  junction 
assuming  oblique  spreading  on  the  Central  and  South¬ 
west  Indian  Ridges.  The  motion  of  the  junction  J 
causes  the  length  of  the  Southeast  Indian  Ridge  to 
decrease  only  very  slowly  compared  to  the  situation 

in  Figure  7a. 

d.  A  schematic  representation  of  the  trace  of  the  triple 
junction  on  the  three  plates. 


Thus  as  one  goes  back  in  time  and  reverses  the  effect,  the 
Southeast  Indian  Ridge  should  increase  in  length.  This  is  not 
observed;  on  the  contrary,  the  Southeast  Indian  Ridge  has 
remained  almost  exactly  the  same  length  between  44  Ma  and 
present. 

There  are  two  simple  ways  in  which  the  Indian  Ocean 
triple  junction  can  evolve  and  maintain  the  Southeast  Indian 
Ridge  at  an  almost  constant  length.  In  the  first  the  RRR 
junction  develops  smoothly  for  a  short  time  as  J  in  Figure  7a. 
Then  the  triple  junction  jumps  northwest  up  the  southeast  branch 
and  fracture  zones  are  formed  on  the  Southwest  and  Central 
Indian  Ridges.  The  frequency  and  magnitude  of  these  jumps 
is  such  that  the  length  of  the  Southeast  Indian  Ridge  never 
changes  by  a  significant  amount.  Because  the  jump  causes  some 
of  the  African  plate  to  be  added  to  the  Antarctic  plate,  the 
traces  of  this  jump  are  more  pronounced  on  the  Antarctic  plate 
than  on  the  African  plate.  Also  since  some  of  the  crust 
created  on  the  African  plate  by  the  Central  Indian  Ridge  is 
added  to  the  Indian  plate,  the  jumps  leave  a  distinctive  zig-zag 
shape  of  the  triple  junction  on  the  Indian  plate.  A  schematic 
development  of  the  triple  junction  is  shown  as  Figure  7b. 

Note  the  zig-zag  nature  of  the  east-west  boundary  of  the  trace 
of  the  triple  junction  on  the  African  plate  and  the  extra  crust 
observed  on  the  Antarctic  plate. 

An  alternate  and  mechanically  more  simple  scheme  for 
the  development  of  the  triple  junction  is  for  spreading  on 
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the  Central  Indian  and  Southwest  Indian  Ridges  to  be  oblique 
rather  than  at  right  angles  to  the  transform  faults.  In  this 
case  the  extension  of  the  three  ridge  axes  can  be  made  to  fit 
close  to  the  mid-point  of  the  motion  between  India  and 
Antarctica  and  the  velocity  triangle  resembles  Figure  7c. 

This  triple  junction  is  stable  and  the  Southeast  Indian  Ridge 
becomes  shorter  only  very  slowly  as  time  passes.  The  trace 
of  the  triple  junction  on  the  three  plates  is  simple,  with  no 
necessity  for  jumps  and  offsets  (Figure  7d) .  In  the  actual 
Southwest  Indian  Ridge,  oblique  spreading  does  not  occur  well 
away  from  the  triple  junction.  It  is  possible  that  the 
obliquity  could  decrease  with  increasing  distance  from  the 
triple  junction.  One  can  obtain  the  same  result  by  asymmetric 
rather  than  oblique  spreading  on  the  Southwest  Indian  and 
Central  Indian  Ridges.  Such  strong  asymmetry  can  be  ruled 
out  by  the  symmetric  magnetic  anomalies  observed  on  the 
Central  Indian  and  Southwest  Indian  Ridges. 

On  the  present  Southwest  and  Indian  Ridge  the  topographic 
data  (Figure  1)  suggest  strongly  that  there  is  oblique 
spreading  at  61 °E  to  63°E  and  normal  spreading  at  most  other 
places,  except  possibly  close  to  the  triple  junction.  Further, 
there  appears  to  have  been  both  oblique  and  normal  spreading 
on  the  Southwest  Indian  Ridge  in  the  past.  Though  it  will 
be  difficult  to  differentiate  between  these  two  extreme 
cases,  and  in  the  limit  of  small  ridge  jumps  they  are  the  same, 
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it  is  worthwhile  to  enumerate  some  of  the  differences 
expected  if  the  ridge  jumps  are  large. 

In  the  ridge- jumping/normal  spreading  model,  one 
would  expect  the  following  features: 

(1)  as  many  fracture  zones  on  the  Central  Indian  Ridge 
as  on  the  Southwest  Indian  Ridge  during  the  eastward  movement 
of  the  triple  junction. 

(2)  rough  but  continuously  lineated  topography,  associated 
with  the  jump,  on  the  trace  of  the  junction  on  the  Antarctic 
plate. 

(3)  rough  and  discontinuously  lineated  topography  along 
the  trace  on  the  African  plate. 

(4)  less  rough  but  still  observable  topographic  expression 
on  the  Indian  plate. 

In  the  case  of  oblique  spreading,  one  would  anticipate 

(a)  smooth  traces  on  all  three  plates  of  the  junction.  Further 

(b)  the  anomalies  on  the  Southwest  Indian  Ridge  would  be 
subdued  and  oblique  near  the  traces  of  the  triple  junction  and 
finally,  (c)  the  anomalies  on  the  Central  Indian  Ridge  would 
be  oblique  near  the  traces  of  the  triple  junction.  This 
obliquity  would  be  small  and  may  not  be  detectable  without 

a  detailed  survey. 


72 


CONCLUSIONS 

Between  54°E  and  61°E  the  Southwest  Indian  ridge  is 
dominated  by  a  suite  of  en  echelon  north-south  topographic 
gashes.  East  of  61°E  to  the  Indian  Ocean  triple  junction 
the  major  feature  is  a  triangular  region  of  rough  topography 
with  a  sinuous  northeast  trending  depression  running  from 
the  north-south  trending  gashes  to  the  eastern  apex.  We 
show  that  this  depression  probably  marks  the  axis  of  a 
spreading  center  and  that  spreading  occurs  approximately 
north-south  at  a  rate  of  0.8  cm/yr.  Further,  we  have 
found  that  the  directions  of  motion  predicted  for  the  Africa/ 
Antarctica  boundary  at  the  Bouvet  and  Indian  Ocean  triple 
junctions  are  compatible.  Thus  the  Southwest  Indian  Ridge 
is  a  continuous  boundary  between  only  two  plates.  The  directions 
and  rates  of  motion  at  both  triple  junctions  and  those  inferred 
from  transform  faults  and  magnetic  anomalies  between  them  give 
a  pole  for  the  present  motion  of  Africa/Antarctica  at  8.4°N, 

42 . 4°W.  This  pole  is  close  to  that  of  Chase  (1972)  but  more 
than  30°  from  that  predicted  by  Minster  and  others  (1974)  in 
their  respective  global  closure  models. 

A  major  change  of  spreading  direction  from  north-south 
is  northeast-southwest  occurs  in  the  Indian  Ocean  at  the 
time  of  Anomaly  19  (44  Ma) .  We  have  shown  that  from  this 
time  until  the  present  the  position  of  the  poles  of  relative 
motion  between  the  African,  Indian,  and  Antarctic  plates 
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change  remarkably  little.  We  use  this  information  and  the 
magnetic  anomalies  to  construct  a  tectonic  history  of  the 
three  major  plate  boundaries,  the  Central,  Southwest,  and 
Southeast  Indian  Ridges,  at  the  time  of  anomalies  19,  8,  and 
5  (44  Ma,  27  Ma,  and  9  Ma;  Figure  8) .  In  our  initial 
reconstruction  the  Central  and  Southeast  Indian  Ridges  meet 
the  Southwest  Indian  Ridge  at  the  Melville  Fracture  Zone 
at  61 °E.  The  right  triangular  region  of  rough  topography 
between  this  fracture  zone  and  the  triple  junction  is  created 
by  the  rapid  eastwards  motion  of  this  junction  and  slow 
spreading  on  the  Southwest  Indian  Ridge.  In  the  development 
of  the  triple  junction  the  Southeast  Indian  Ridge  remains 
almost  constant  in  length.  With  present  information  it  is 
not  possible  to  determine  whether  this  results  predominantly 
from  spreading  normal  to  the  transform  faults  on  the 
Southwest  Indian  Ridge  and  continuous  jumps  to  keep  the 
Southeast  Indian  Ridge  constant  or  to  oblique  spreading  on 
both  the  Southwest  Indian  and  Central  Indian  Ridges.  We 
believe  that  the  triple  junction  probably  evolves  in  both 
modes . 
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FIGURE  8 

The  position  of  the  Central,  Southeast,  and  Southwest 
Indian  Ridges  with  respect  to  a  fixed  Africa  at  the 
time  of  anomalies  19(A),  8(B),  5(C),  and  at  present  (D) . 
Also  shown  are  the  rotated  and  present  positions  of  the 
observed  anomalies  (symbols  as  for  Figure  6) .  The 
continuous  lines  represent  the  inferred  position  of 
anomalies  19, 8, and  5  on  the  reconstructions.  They 
were  determined  by  rotating  the  individual  ridge  axes 
through  the  appropriate  rotation  angle.  On  the  Anomaly 
19  reconstruction  (A)  the  magnetic  anomalies  from  the 
Indian  plate  (black)  overlap  those  from  the  Antarctic 
plate  but  do  not  match  the  inferred  position  of  the 
ridge  axis.  This  implies  that  the  pole  we  used  for 
India/Antarctica  for  this  reconstruction  is  reasonable 
and  that  more  crust  has  been  created  south  of  the  axis 


than  to  the  north. 
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CHAPTER  3 
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ABSTRACT 

The  boundaries  of  three  plates,  Africa,  India,  and 
Antarctica,  meet  in  a  triple  junction  in  the  Indian  Ocean 
near  25°S,  70°E.  Using  observed  bathymetry  and  magnetic 
anomalies,  we  locate  the  junction  to  within  5  km  and  show 
that  it  is  a  ridge-ridge-ridge  type.  Two  of  the  ridges  are 
normal  spreading  centers  with  well-defined  median  valleys, 
but  one,  the  Southwest  Indian  Ridge  (SWIR) ,  has  a  peculiar 
morphology  in  the  vicinity  of  the  triple  junction,  that  of 
an  elongate  triangular  deep  with  the  junction  at  its  apex. 
The  floor  of  this  deep  represents  crust  formed  at  the  SWIR, 
and  the  morphology  results  from  the  evolution  of  the  triple 
junction.  We  determine  directions  and  rates  of  spreading  on 
all  three  ridges,  and  poles  of  relative  motion  for  all  plate 
pairs.  Though  we  cannot  determine  the  precise  stability 
conditions  at  the  triple  junction,  we  are  able  to  map  the 
evolution  of  the  junction  over  the  last  10  million  years. 

We  detect  and  discuss  topographic  expressions  of  the  triple 
junction  traces  on  the  three  plates. 
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INTRODUCTION 

The  central  Indian  Ocean  triple  junction,  the  meeting 
point  of  the  African,  Indian,  and  Antarctic  plates,  is 
located  near  25°S,  70°E.  Previously,  little  has  been  known 
about  the  area  in  which  the  Central  (CIR) ,  Southeast  (SEIR) 
and  Southwest  (SWIR)  Indian  ridges  join,  as  there  have  been 
few  properly  oriented  ship  tracks  through  the  region.  Until 
recently,  there  has  not  even  been  any  reliable  indirect 
evidence  (through  global  closure  calculations)  of  the  local 
spreading  direction  of  the  SWIR  (Minster  et  al.,  1974;  Chase, 
1972) .  Since  the  SWIR  is  spreading  slowly,  at  less  than 
10  mm/yr  (Bergh  and  Norton,  1976) ,  its  data  are  often  difficult 
to  interpret,  and  the  more  well-constrained  observations  possible 
at  the  triple  junctions  at  either  end  can  be  important  additions 
to  our  knowledge.  It  was  largely  for  this  reason  that  Sclater 
et  al ,  (1976)  studied  the  Bouvet  Triple  Junction  at  the  South 

Atlantic  terminus  of  the  SWIR,  and  that  we  began  the  present 
study  of  the  Indian  Ocean  Triple  Junction. 

The  data  presented  here  come  principally  from  legs  5  and  6 
of  R/V  Atlantis  II  cruise  93  from  February  through  April  of 
1976.  The  survey  was  designed  to  delineate  the  tectonic 
history  within  250  km  of  the  triple  junction  by  running  five 
or  more  lines  spaced  at  a  few  tens  of  kilometers  across  each 
of  the  three  ridges  approximately  parallel  to  their  spreading 


directions.  By  utilizing  already  existing  data  as  well 
(Figure  1) ,  we  have  compiled  a  detailed  geophysical  survey 
within  a  roughly  500  km  square  around  the  triple  junction. 

Using  topographic  and  magnetic  data  from  these  cruises,  we 
locate  the  tr:  'le  junction  to  within  5  km  and  show  that  the 
instantaneous  relative  velocity  triangle  closes  and  that  the 
junction  is  of  the  ridge-ridge-ridge  (RRR)  type.  We  then 
examine  the  tectonic  history  of  the  region  over  the  past  10 
million  years,  showing  that  it  is  consistent  with  observations 
elsewhere  along  the  same  plate  boundaries  and  in  particular 
with  those  at  the  Bouvet  Triple  Junction.  Finally,  we  discuss 
the  implications  of  the  observed  relative  motions  of  the  plates 
and  the  junction  over  the  past  10  million  years. 

TOPOGRAPHY 

Topographic  information  from  recent  satellite  navigated 
cruises  of  the  Scripps  Institute  of  Oceanography  and  the 
Institute  de  Physique  du  Globe  de  Paris  have  been  combined  with 
that  of  the  Atlantis  II  survey.  To  these  have  been  added  older 
tracks  from  the  SIO  expeditions  DODO  (1964)  and  CIRCE  (1968) 
and  the  Institute  de  Physique  du  Globe.  We  have  contoured 
the  data  in  intervals  of  500  m  corrected  for  the  local  sound 
velocity  profile,  and  the  resulting  topographic  chart 
(Figure  2)  shows  several  features  of  interest  to  good  advantage. 
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FIGURE  1 

The  ship  tracks  on  which  we  base  our  study  of  the 
Indian  Ocean  triple  junction.  The  solid  lines  represent 
the  AII-93  track,  and  the  dotted  lines  represent  other 
tracks.  The  open  circles  with  crosses  denote  teleseis- 
mically  located  earthquake  epicenters  (1964-1976)  and 
indicate  the  general  location  of  the  ridge  axes.  The 
letters  designate  AII-93  profiles  discussed  in  the  text. 


FIGURE  2 


Bathymetry  in  the  area  of  the  Indian  Ocean  triple 
junction.  Depths  are  in  corrected  meters,  and  the 
contour  interval  is  500  m.  Ship  tracks  are  shown  as 
dotted  lines.  Contours  by  R.  L.  Fisher. 
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One  can  clearly  see  the  three  plate  boundaries  meeting  in  a 
triple  junction  near  25.7°S,  70°E,  and  the  junction  appears  at 
first  glance  to  be  of  the  RRR  type.  The  morphologies  of  the 
three  spreading  centers  are  of  markedly  different  character, 
however . 

The  SEIR,  which  separates  the  Indian  and  Antarctic  plates, 
has  the  morphology  of  a  "classic"  spreading  ridge.  It  is 
characterized  by  smooth,  blocky  relief  trending  sub-parallel  to 
a  well-defined  median  valley  offset  at  intervals  by  deep  linear 
fracture  zones.  These  fracture  zones,  which  dominate  the 
tectonic  grain  of  the  ridge,  trend  about  N45°E  and  occur  at 
70.5°E,  71°E ,  and  73°E.  They  offset  the  ridge  axis  in  a 
left- lateral  sense  and  are  associated  with  about  2000  m  of 
relief.  The  transform  at  71<>E  has  a  60  km  offset  and  constrains 
the  local  direction  of  relative  motion  between  the  Indian  and 
Antarctic  plate  to  be  within  2°  of  047°.  Between  the  fracture 
zones,  and  oriented  perpendicular  to  their  trend,  lie  linear 
sections  of  spreading  ridge,  their  axes  marked  by  a  median 
valley  about  500  m  deep  (Figure  3a)  and  well-delineated  by 
the  3500  m  isobath  (Figure  2) .  Additional  minor  transform 
faults  may  offset  the  median  valley  by  small  amounts.  Beyond 
the  rift  valley  mountains,  the  sea  floor  deepens  gradually 
in  a  subdued,  blocky  fashion  from  about  2500  m  depth  to  about 
3500  m  depth  300  km  from  the  axis.  We  have  compared  the  profile 


FIGURE  3a 


Bathymetric  and  magnetic  profiles  across  the  Southeast 
Indian  Ridge  projected  onto  a  heading  of  N45°E.  See 
Figure  1  for  profile  locations.  The  expected  range 
of  bathymetric  values  for  a  ridge  spreading  at  30  mm/yr 
according  to  the  relation  of  Parsons  and  Sclater  (1977) 
is  shown.  Also  shown  are  the  magnetization  sequence 
according  to  the  time  scale  of  LaBrecque  et  al .  (1977) 
and  the  phase  shifted  anomaly  sequence  calculated  for 
a  flat  layer  source  extending  from  3500  to  4000  m  depth, 
both  before  and  after  convolution  with  a  Gaussian  filter 
of  2.5  km  standard  deviation.  The  stippled  band  through 
the  magnetic  and  bathymetric  profiles  indicates  the 
inferred  region  of  Bruhnes  age  crust.  FZ  indicates 
an  inferred  fracture  zone. 
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of  the  ridge  to  the  age-depth  relation  of  Parsons  and  Sclater 
(1977)  and  find  them  consistent  to  within  the  300  m  amplitude 
of  the  short  wavelength  relief.  In  brief,  the  SEIR  has  all 
the  morphologic  aspects  of  an  archetypal  mid-ocean  spreading 
center. 

The  topography  of  the  CIR  is  more  rugged  and  variable  and 
less  well-defined  than  that  of  the  SEIR.  Since  it  is  spreading 
more  slowly  than  the  SEIR  (25  mm/yr  as  compared  to  30  mm/yr) 
this  is  in  keeping  with  the  general  observation  that  rougher 
topography  is  associated  with  slower  spreading  ridges.  The 
median  valley  of  the  CIR,  delineated  by  the  3500  m  contour 
and  trending  approximately  N30°W,  is  clearly  visible  from  the 
triple  junction  to  24.5°S,  growing  generally  shallower  toward 
the  north  (Figure  2).  North  of  24.5°S  the  median  valley  may 
still  be  traced,  but  less  obviously.  Bathymetric  profiles 
across  the  ridge  (Figure  3b)  show  that  the  valley  exists  and 
is  200  m  to  300  m  deep,  but  that  its  width  and  its  absolute 
depths  are  quite  variable,  making  it  difficult  to  follow  on 
a  contour  chart.  Near  25.2°S,  70°E  the  CIR  median  valley 
is  offset  20  km  in  a  right-lateral  sense  by  a  small  transform 
fault  trending  059°  +  4°.  This  gives  a  reasonably  well 
constrained  measure  of  the  local  direction  of  relative  motion 
between  the  African  and  Indian  plates.  There  are  no  other 
clearly  defined  transform  faults  on  the  CIR  within  300  km 
of  the  triple  junction,  though  at  22.5°S,  just  to  the  north 


FIGURE  3b 


As  in  Figure  3a,  but  the  profiles  are  across  the  Central 
Indian  Ridge  and  are  projected  onto  a  heading  of  N60°E. 
The  models  are  for  a  spreading  rate  of  25  mm/yr.  The 
dashed  portions  of  profiles  E  are  over  crust  thought  to 
have  been  formed  at  the  Southeast  and  Southwest  Indian 
Ridges . 
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of  the  area  of  Figure  2,  lies  the  southernmost  of  the  major 
fracture  zones  mapped  by  Fisher  et  al.  (1971) .  Between  this 
and  the  25.2°S  transform  fault  are  two  offsets  of  the  median 
valley  at  24.2°S  and  24.6°S.  These  offset  the  axis  by  only 
10  km  and  have  little  topographic  expression,  though  their 
existence  is  confirmed  by  magnetic  data.  In  summary,  the 
morphology  of  the  CIR  is  similar  to  that  of  SEIR.  It  consists 
of  segments  of  spreading  centers  (with  median  valleys)  oriented 
at  right  angles  to  the  direction  of  relative  plate  motion. 

The  relatively  confused  topography,  however,  does  not  allow 
one  to  rule  out  a  slight  obliquity  of  spreading,  perhaps  by 
as  much  as  10°. 

The  morphology  of  the  SWIR  stands  in  marked  contrast  to 
that  of  the  CIR  and  the  SEIR.  It  appears  as  a  great  triangular 
chasm  whose  floor,  at  about  4500  m  depth,  is  flanked  north  and 
south  by  2000  m  high  linear  mountain  ridges  (Figure  2) .  The 
ridges  converge  eastward  toward  the  most  acute  apex  of  the 
triangle,  located  at  the  triple  junction.  The  extreme  relief 
across  this  spreading  center  is  seen  dramatically  in  profile 
(Figure  3c) .  The  active  spreading  center  and  its  recently 
created  sea  floor  are  up  to  2000  m  deeper  than  predicted  by 
the  Parsons  and  Sclater  (1977)  equation.  The  precise  positions 
and  trends  of  ridge  axis  segments  are  not  at  all  well-defined 
by  the  rough  topography  of  the  SWIR,  though  the  general 
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FIGURE  3c 


As  in  Figure  3a,  but  across  the  Southwest  Indian  Ridge 
and  projected  onto  a  heading  of  N0°E.  The  models  are 
based  on  a  spreading  rate  of  8  mm/yr,  and  the  magnetized 
layer  is  taken  to  lie  between  4000  and  4500  m  depth. 

The  dashed  portions  of  the  profiles  are  over  crust 
thought  to  have  been  formed  at  the  Central  and  Southeast 
Indian  Ridges. 
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topographic  grain  trends  about  N70°E.  Neither  are  there  any 
obvious  fracture  zone  traces  within  300  km  of  the  triple 
junction  with  which  the  local  direction  of  relative  motion 
between  the  African  and  Antarctic  plates  might  be  determined. 
There  are  a  few  vague  indications  of  a  north-south  grain,  the 
most  prominent  being  the  N3°W  trending  widening  of  the  central 
deep  at  68.8°E,  but  these  are  no  more  than  merely  suggestive. 
Farther  to  the  west,  however,  Sc later  et  al .  (in  preparation) 
have  mapped  the  Atlantis  II  and  Melville  fracture  zones  at 
57.5®  and  60.5®E  respectively.  These  are  long,  deep  fracture 
zone  traces  of  large  offset  and  trend  within  2°  of  north-south, 
and  this  direction  of  relative  plate  motion  may  be  taken  to  be 
close  to  that  nearer  the  triple  junction.  Assuming  this 
direction  of  spreading,  it  is  not  obvious  from  the  morphology 
whether  the  SWIR  is  dominated  by  en  echelon  ridge  segments 
spreading  normal  to  their  strike  or  by  segments  spreading 
obliquely.  Both  interpretations  can  be  supported  by  various 
aspects  of  the  rough  topography. 

The  confluence  of  these  three  spreading  centers,  each 
of  different  character,  marks  the  site  of  the  triple  junction 
near  25.7°S,  70°E.  Here  the  continuous  median  valley  of  the 
Central  and  Southeast  Indian  Ridges  (delineated  by  the  3500  m 
contour  in  Figure  2)  deepens  and  changes  trend  in  a  distinct 
bend.  Just  at  this  bend,  the  4500  m  depths  of  the  SWIR 
approach  the  CIR-SEIR  median  valley,  separated  from  it  by  a 
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sill  of  about  3300  m  depth  in  a  narrow  east-west  trending 
channel.  The  intersection  of  this  channel  and  the  CIR-SEIR 
median  valley  marks  the  site  of  the  triple  junction. 

MAGNETIC  ANOMALIES 

During  R/V  Atlantis  II  cruise  93  the  total  geomagnetic 
field  strength  was  observed  continuously.  From  this  we 
subtracted  the  International  Geomagnetic  Reference  Field 
(IGRF)  to  obtain  the  magnitude  of  the  residual  geomagnetic 
field,  profiles  of  which  across  the  three  ridges  are  shown  in 
Figures  3a,  3b,  and  3c.  For  each  ridge,  the  profiles  have 
been  projected  along  the  approximate  direction  of  spreading. 
Also  displayed  is  a  synthetic  magnetic  anomaly  profile  for 
each  ridge  derived  from  the  geomagnetic  reversal  time  scale  of 
LaBrecque  et  al .  (1977)  and  a  two-dimensional  flat  layer  model 
of  500  m  thickness  at  the  average  depth  of  each  ridge.  We  have 
modeled  the  "contamination"  of  the  magnetic  signal  caused  by 
the  finite  width  of  the  active  zone  of  intrusive  and  extrusive 
addition  of  crustal  material  (Atwater  and  Mudie,  1973)  by 
convolving  the  block  model  magnetization  with  a  gaussian 
filter  of  standard  deviation  similar  to  that  of  the  injection 
process  (Schouten  and  Denham,  in  preparation;  Tisseau,  1978) . 
The  effect  of  contamination  is,  of  course,  greater  at  lower 
spreading  rates  because  of  the  narrower  width  of  single¬ 
polarity  blocks.  One  also  expects  the  standard  deviation  of 
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the  gaussian  filter  to  increase  with  decreasing  spreading  rate 
as  the  plate  motions  become  slower  and  perhaps  therefore  more 
episodic,  increasing  the  temporal  contribution  to  the 
variability  of  the  magnetization  (Schouten  and  Denham,  in 
preparation) . 

The  magnetic  anomalies  along  the  SEIR  are  excellently 
formed  (Figure  3a)  in  that  one  may  easily  correlate  and 
identify  them  almost  everywhere.  The  only  notable  exception 
is  the  disturbance  of  the  central  portion  of  profile  M  which 
may  be  explained  by  the  probable  occurrence  of  a  fracture  zone. 
The  standard  deviation  of  the  gaussian  contamination  is  quite 
small,  being  2  km  or  less,  and  the  phase  shift  of  the  magnetic 
anomalies  is  close  to  the  50°  value  expected  for  blocks 
oriented  along  N45°W  at  this  location.  On  the  northeastern 
ends  of  profiles  0  and  P  one  sees  the  complete  anomaly  sequence 
since  Anomaly  5  (9  Ha).  Spreading  has  been  nearly  symmetric  at 
a  half-rate  close  to  30  mm/yr  during  this  period  and  particularly 
since  the  time  of  Anomaly  2'  (3  Ma).  During  the  latter  period 
the  only  evident  complications  are  the  previously  mentioned 
fracture  zone  in  profile  M  and  an  inferred  ridge  jump  in 
profile  P  at  or  about  the  time  of  Anomaly  2.  The  period 
previous  to  3  Ma  is  only  seen  well  on  profiles  O  and  P. 

There  were  slight  asymmetries  of  spreading  and  a  slightly 
higher  average  spreading  rate  than  more  recently,  between  31 
and  32  mm/yr.  Overall,  however,  the  profiles  show  that 


magnetically,  as  well  as  morphologically,  the  SEIR  is  a 
classic  example  of  a  simple  symmetric  spreading  center. 

Just  as  the  CIR  is  morphologically  more  complex  than 
the  SEIR,  so  also  are  its  magnetic  anomalies  more  difficult 
to  interpret  than  those  of  the  SEIR.  They  are  less  distinct 
and  are  of  more  variable  shape  {Figure  3b) ,  and  seem  consistent 
with  a  zone  of  emplacement  with  a  standard  deviation  of  between 
2  and  3  km.  Their  phase  shift  is  impossible  to  determine 
accurately,  but  appears  close  to  the  30°  expected  of  a  ridge 
trending  N30°W.  In  spite  of  these  difficulties,  many  anomalies 
may  be  correlated  and  identified,  and  as  a  whole  the  profiles 
indicate  that  the  CIR  has  been  spreading  with  little  or  no 
asymmetry  at  a  half  rate  of  about  25  mm/yr  since  the  time  of 
Anomaly  2'  (3  Ma)  or  perhaps  Anomaly  3  (4  Ma) .  Each  profile 
has  some  complications,  however,  and  we  have  little  information 
about  the  spreading  rate  before  4  Ma.  Profile  E  shows  higher 
spreading  rate  and  considerable  asymmetry,  but  it  may  cross  a 
transform  fault  near  its  center,  and  its  outer  portions  are 
over  crust  formed  at  the  Southwest  and  Southeast  Indian  spreading 
centers.  Profile  D  runs  obliquely  across  the  25.2°S  transform 
fault.  Lines  C,  B,  and  A  indicate  symmetric  spreading  at  about 
25  mm/yr  since  time  of  Anomaly  2'  (3  Ma)  or  perhaps  Anomaly  3 
(4  Ma) .  The  western  end  of  profile  A  extends  to  Anomaly  5 
(9  Ma) .  Beyond  Anomaly  3'  (5.3  Ma) ,  it  indicates  a  slower 
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spreading  rate  of  20  mm/yr,  in  excellent  agreement  with  the 
rate  from  21  Ma  to  9  Ma  that  Sclater  et  al.  (in  preparation) 
find  for  the  western  side  of  the  CIR  (based  in  part  on  the 
same  line).  It  is  questionable,  however,  whether  this 
reflects  a  slower  rate  of  opening  or  is  an  effect  of 
asymmetric  spreading.  The  data  of  Fisher  et  al .  (1971)  do 
not  indicate  any  significant  asymmetry,  but  neither  do  they 
indicate  any  significant  change  in  rate.  Over  the  last 
4  million  years,  in  any  case,  the  Central  Indian  Ridge  near 
the  triple  junction  has  been  spreading  symmetrically  at  near 
constant  rate. 

The  magnetic  anomalies  along  the  SWIR  are  very  subdued 
and  difficult  to  interpret  (Figure  3c) .  The  central  anomaly 
is  clearly  recognizable  only  on  profiles  U  and  V,  as  is 
(probably)  Anomaly  2'  both  to  the  north  and  the  south  of  the 
axis.  Indications  of  a  more  subdued  central  anomaly  are 
present  on  profiles  W  and  X.  These  lines  are  consistent  with 
observations  made  on  the  SWIR  between  57 °E  and  68°E  by  Sclater 
et  al.  (in  preparation) ,  which  demonstrate  that  the  SWIR  has 
been  spreading  in  a  north-south  direction  at  a  half-rate  of 
8  mm/yr  since  at  least  10  Ma.  We  may  identify  several  of  the 
reasons  for  the  poor  quality  of  the  magnetic  anomalies.  The 
slow  spreading  rate,  particularly  when  coupled  with  gaussian 
contamination  of  3  to  4  km  standard  deviation,  leads  to  low 
amplitude,  indistinct  magnetic  anomalies.  The  phase  shift  is 
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high,  consistent  with  the  75°  expected  for  a  normally  spreading 
ridge.  Either  closely  spaced  fracture  zones  must  exist, 
disturbing  the  two-dimensional  magnetic  patterns,  or  the 
ridge  must  spread  obliquely  near  the  triple  junction  (also 
with  the  possibility  of  fracture  zones) .  This  would  have 
the  effect  of  reducing  the  spreading  rate  normal  to  ridge 
strike  and  thus  decrease  the  resolution  of  the  anomalies. 

The  phase  shift  expected  for  this  configuration  would  be 
within  a  few  degrees  of  that  for  the  normal  spreading  case. 

In  summary,  the  profiles  across  the  SWIR,  together  with 
those  of  Sclater  et  al .  (in  preparation) ,  indicate  that  the 
African-Antarctic  relative  motion  has  been  progressing  at  a 
half  rate  of  close  to  8  mm/yr  near  the  triple  junction.  Since 
10  Ma,  this  rate  has  changed  very  little,  and  spreading  has 
been  nearly  symmetric,  but  the  data  do  not  allow  confident 
measurement  of  departures  from  these  norms.  Finally,  it  is 
interesting  to  note  that  only  the  central  portions  of  the 
profiles  of  Figure  3c  are  over  crust  created  at  the  Southwest 
Indian  spreading  center,  and  this  portion  narrows  as  one 
approaches  the  triple  junction.  Crust  under  the  southern  ends 
of  the  profiles  was  created  at  the  SEIR,  and  the  magnetic 
profiles  Z,  Y,  X  and  W  show  respectively  parts  of  anomalies 
2',  3,  3'  and  possibly  4  created  at  that  ridge. 
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TECTONICS 

The  magnetic  anomaly  profiles  of  the  AII-93  lines  and 
of  other  cruises  in  the  area  are  plotted  at  right  angles  to 
the  ship's  track  and  are  displayed  in  Figure  4  superimposed 
on  selected  bathymetric  contours.  Also  shown  are  identifications 
of  anomalies  1  through  5,  as  well  as  the  locations  of  fracture 
zones  inferred  from  the  bathymetry  and  magnetics. 

The  local  evolution  of  the  Southeast  Indian  Ridge  is 
well-defined  by  the  data  and  appears  to  have  been  predominantly 
one  of  symmetric  spreading  at  right  angles  to  ridge  strike  at 
nearly  constant  rate  and  direction.  Here  the  tectonic  history 
is  displayed  almost  in  its  entirety,  and  such  subtleties  as 
the  pre-anomaly  2  ridge  jump  at  72 °E  are  evident. 

The  CIR  evolution,  though  generally  well  laid  out,  is 
somewhat  less  clear.  The  line  placement  and  the  quality  of 
the  magnetics  are  such  that  possible  small  obliquities  and 
asymmetries  of  spreading  can  not  be  resolved,  particularly 
in  the  older  anomalies  and  very  near  the  triple  junction. 

This  question  has  profound  implications  for  the  evolution  of 
the  triple  junction. 

Except  for  the  region  of  AII-93  profile  V  near  68.4°E, 
the  tectonic  structure  of  the  SWIR  is  not  readily  resolved. 

One  can  form  some  idea  of  the  position  of  the  spreading  axis 
from  the  magnetics  and  bathymetry  and  map  its  approach  to  the 
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FIGURE  4 


Magnetic  anomaly  profiles  along  track  superimposed  on 
the  generalized  bathymetry  in  the  triple  junction  area. 
The  lightly  shaded  regions  are  shallower  than  3000  m, 
and  the  stippled  regions  are  deeper  than  3500  m  in  the 
Central  and  Southeast  rift  valleys,  and  deeper  than 
4500m  in  the  Southwest  rift  valley.  Anomaly  identifi¬ 
cations  are  indicated  by  heavy  lines  and  numbers. 
Fracture  zones  are  also  indicated.  The  heavy  dashed 
lines  represent  the  approximate  extent  of  crust  formed 
at  the  Southwest  Indian  spreading  center.  The  earth¬ 
quake  epicenters  shown  in  Figure  1  appear  as  circles. 
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CIR-SEIR  axis.  Some  fracture  zones  are  indicated,  but  the 
actual  trend  of  the  spreading  axis  is  not  obvious.  The  extent 
of  the  SWIR  regime  is  constrained  by  the  presence  of  CIR  and 
SEIR  generated  crust  to  the  north  and  south  of  the  southwest 
axis.  The  boundaries  of  this  triangular  region  on  the 
African  and  Antarctic  plates  are  indicated  by  the  rough-smooth 
transition  in  bathymetry  and  the  approximately  congruent 
smooth-rough  transition  in  magnetic  texture. 

The  tectonic  diagram  of  Figure  4  serves  to  locate  the 
triple  junction  within  a  few  kilometers  of  25.66°S,  70.06°E. 

Here  the  SHIR  deep  joins  by  means  of  an  east-west  trending 
channel  with  the  continuous  CIR-SEIR  median  valley  at  the 
site  of  the  bend  in  the  latter.  Congruent  with  the 
morphologic  bend  is  the  angle  in  the  trend  of  the  CIR-SEIR 
central  anomaly.  The  topographic  and  magnetic  data  together 
show  that  the  three  ridges  meet  in  a  RRR  triple  junction  at  that 
point,  and  the  magnetic  data  indicate  a  stable  evolution  of 
the  junction  over  the  past  few  million  years. 

It  is  significant  that  the  tectonic  information  (Figure  4) 
gives  no  indication  whatsoever  of  a  local  system  of  microplates 
such  as  postulated  by  Barberi  and  Varet  (1977)  for  the 
incipient  triple  junction  at  the  head  of  the  East  African 
rift  system;  the  Indian  Ocean  triple  junction  is  indeed  a 
junction  of  three  plates.  It  is  peculiar,  however,  that  the 
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bathymetric  expression  of  the  Southwest  Indian  spreading 
center  is  not  more  clearly  continuous  with  those  of  the 
Central  and  Southeast  branches.  Perhaps  the  more  slowly 
spreading  Southwest  axis  does  not  exist  continuously  as  an 
identifiable  feature  in  the  young,  thin  crust  near  the 
triple  junction,  but  manifests  itself  only  during  especially 
active  pulses  of  motion  between  the  African  and  Antarctic 
plates. 

INSTANTANEOUS  RELATIVE  PLATE  MOTIONS 

Local  Relative  Plate  Motions  and  the  Velocity  Diagram 

Table  1  lists  the  directions  and  rates  of  the  relative 
motions  among  the  African,  Indian,  and  Antarctic  plates 
derived  from  bathymetric  and  magnetic  data  near  the  triple 
junction.  The  large  uncertainties  ascribed  to  the  Africa/ 
Antarctica  values  arise  from  the  fact  that  the  best  information 
along  this  boundary  comes  not  from  the  immediate  neighborhood 
of  the  triple  junction,  but  from  500  to  1500  km  to  the  west. 
Preferred  estimates  of  these  errors  are  about  half  the  given 
magnitude,  but  since  most  of  the  data  are  not  of  local  origin 
one  must  allow  for  greater  possible  discrepancy.  The 
information  in  Table  1  defines  the  velocity  triangle  of  the 
triple  junction  (Figure  5a) .  The  triangle  is  closed  well 
within  the  uncertainties  of  the  measurements,  indicating 
that  the  local  relative  plate  velocities  have  been 


TABLE  1.  Observed  Relative  Motions  at  the  Indian  Ocean  Triple  Junction 


Plate  Boundary 

Source  Data 

Direction 

deg. 

Africa/India 

25.2°S  fracture  zone 
AII-93  profiles  A,  B,  C 

059°  ±  4“ 

India/ 

Antarctica 

71 °E  fracture  zone 

AII-93  profiles 

J,  K,  L,  N,  0 

047°  ±  2° 

Africa/ 

Antarctica 

Melville  and  Atlantis  II 
fracture  zones 

AII-93  profiles  U  and  V 
and  Sclater  et  al. 

(in  preparation) 

000°  ±  10° 

Full  Rate 
mm/yr 


49  ±  3 


59  ±  2 


16  ±  5 


Ill 


FIGURE  5 

a.  A  vector  diagram  of  the  relative  plate  velocities 

at  the  triple  junction.  Stipple  areas  represent  the 
limits  of  the  possible  errors  in  the  vectors.  The 
error  at  the  African  end  of  the  Africa/Antarctica 
vector  is  the  compound  error  of  that  vector  and  the 
India/Antarctica  vector.  The  lines  representing  the 
best  estimates  of  the  vectors  have  been  emphasized. 
Directions  are  given  in  degrees  east  of  north,  and 
full  rates  are  given  in  mm/yr. 

b.  The  heavy  lines  represent  the  vector  diagram  of  the 
relative  plate  velocities  at  the  triple  junction 
calculated  from  the  instantaneous  relative  angular 
velocity  vectors  computed  in  this  paper.  These  are 
superimposed  on  the  best  estimate  vectors 
ciated  regions  of  error  of  Figure  5a. 


and  asso- 


VECTOR  DIAGRAM 


(a)  OBSERVED  VELOCITIES 

INDIA 


ANTARCTICA 


(b)  CALCULATED  VELOCITIES 

INDIA 


ANTARCTICA 


correctly  deduced.  In  viewing  the  velocity  triangle,  one 
sees  clearly  the  striking  similarity  of  the  Africa/India 
and  Antarctica/India  motions  and  the  markedly  contrasting 
slowness  of  the  north-south  Af rica-Antarctica  motion,  even 
allowing  for  errors  in  the  velocities. 

The  Global  Solution  and  the  Instantaneous  Poles 

Under  the  rigid  plate  hypothesis,  one  can  reduce  the 
uncertainties  in  an  observation  of  relative  plate  motion  by 
using  a  larger  set  of  observations  to  determine  the 
instantaneous  relative  angular  velocity  vector  which  describes 
the  relative  rotation  of  the  plates.  Using  the  method  of 
Forsyth  (Forsyth,  1975;  Sclater  et  al. ,  1976),  we  determined 
the  internally  consistent  set  of  relative  rotation  poles  which 
best  fit  published  observations  of  the  direction  and  rate  of 
the  present  relative  motion  across  the  common  boundaries  of 
the  African,  Indian,  Antarctic,  and  South  American  plates 
(Table  2).  These  poles  are  given  in  Table  3,  and  the  values 
one  calculates  for  the  data  using  these  poles  are  included  in 
Table  2. 

In  general,  the  poles  determined  here  fall  quite  close 
to  other  published  poles  (Table  3) ,  the  only  notable  exception 
being  the  Africa/Antarctica  pole,  for  which  there  is  considerable 
dispersion  among  various  determinations.  Our  pole  at  5.6°N, 
39.4°W  falls  35°  from  that  calculated  by  Minster  et  al. 


TABLK  2. _ Computation  of  Instantaneous  Poles _ 

_ Observed _  Calculated 

Plate  Position  Direction  Rate  Un-  Direction  Rate 

Boundary  Source  N.Lat.  E.Long. _ 0 _ tnm/yr  certainty _ 0 _ mm/yr 


Africa/  2  11.0 

India  3  -9.0 

3  -14.0 

3  -16.0 

3  -17.5 

3  -20.0 

3  -20.0 

1  -25.2 

4  7.0 

4  4.5 

4  4.0 

3  -10.25 

3  -12.0 

3  -13.0 

3  -20.0 

3  -20.75 

3  -22.25 

1  -23.8 

3  -24.0 

1  -24.5 

1  -24.75 

3  -25.0 

India/  1  -26.2 

Antarctica  4  -38.0 

4  -39.0 

4  -40.0 

4  -42.0 

4  -45.0 

5  -49.0 

6  -62.0 

1  -25.9 

1  -26.2 

1  -26.4 

4  -27.0 

4  -31.0 

4  -38.0 

4  -44.0 

7  -50.0 

5  -49.0 

7  -50.5 

8  -62.0 


57.5 

30 

67.5 

52 

66.0 

58 

66.5 

60 

66.0 

63 

67.0 

60 

68.0 

63 

70.0 

59 

60.0 

26 

62.25 

28 

63.5 

29 

66.5 

38 

66.0 

38 

66.5 

38 

66.0 

41 

68.0 

45 

69.0 

46 

69.6 

50 

70.0 

51 

69.8 

50 

69.8 

50 

70.0 

49 

71.0 

47 

78.0 

48 

78.0 

42 

82.0 

39 

85.0 

34 

96.0 

29 

125.0 

7 

156.0 

-22 

70.35 

58 

71.5 

59 

71.9 

59 

73.0 

64 

76.0 

60 

78.0 

70 

94.0 

74 

114.0 

74 

125.0 

74 

134.0 

74 

158.0 

68 

5 

28 

12 

4 

52 

36 

3 

59 

39 

5 

60 

41 

3 

62 

42 

10 

62 

45 

10 

60 

45 

5 

61 

50 

4 

38 

17 

4 

40 

21 

4 

39 

22 

4 

55 

36 

4 

58 

38 

4 

58 

39 

4 

63 

44 

5 

61 

46 

3 

60 

48 

2 

60 

49 

4 

60 

49 

2 

60 

50 

2 

61 

50 

4 

61 

50 

5 

46 

60 

15 

45 

68 

15 

45 

68 

15 

42 

70 

15 

40 

71 

15 

32 

73 

5 

11 

72 

15 

-21 

64 

2 

47 

60 

2 

46 

61 

2 

45 

61 

6 

45 

62 

6 

44 

65 

4 

45 

68 

5 

33 

72 

4 

18 

73 

4 

11 

72 

6 

3 

70 

4 

-23 

64 

TABLE  2  (continued) 
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Observed  Calculated 


Plate 

Position 

Direction 

Rate 

Un- 

Direction 

Rate 

Boundary  Source  N.Lat. 

E.Long. 

0 

mm/yr 

certainty 

0 

mm/yr 

Africa/ 

9 

-54.25 

2.0 

45 

2 

45 

‘16 

Antarctica 

10 

-54.25 

6.0 

40 

2 

41 

16 

10 

-53.  5 

9.0 

39 

2 

38 

16 

11 

-46.0 
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(1974)  and  20°  from  that  calculated  by  Forsyth  (reported  in 
Sclater  et  al . ,  1976) .  It  is  remarkably  close  to  that 
derived  by  Minster  and  Jordan  (1978)  and  Sclater  et  al . 

(in  preparation),  being  5°  and  4°  distant  from  these, 
respectively.  The  calculation  in  the  latter  paper  is  based 
upon  much  the  same  SWIR  data  as  the  present  calculation, 
though  the  method  of  computation  is  different.  As  these 
three  calculations  are  to  date  the  only  ones  based  upon  all 
the  currently  available  data  along  the  Africa/Antarctica  plate 
boundary,  and  in  view  of  the  fact  that  they  arrive  at  virtually 
the  same  Africa/Antarctica  pole  in  spite  of  differing  methods 
and  constraints,  one  must  accept  that  pole  in  preference  to 
earlier  determinations. 

Finally,  we  note  that  a  better  fit  to  the  data  could 
not  be  obtained  by  allowing  Africa  to  behave  as  independent 
Nubian  and  Somalian  plates  with  slight  extension  across  the 
East  African  rift  system.  This  is  in  agreement  with  the 
findings  of  Sclater  et  al .  (in  preparation) .  Evidently,  any 
incipient  plate  separation  in  East  Africa  has  not  yet  had 
an  impact  on  the  global  tectonic  system  measureable  at  the 
scales  investigated  here.  We  did  not  consider  the  possibility 
that  the  Indian  plate  might  in  fact  be  two  plates  in  slight 
relative  motion  across  the  Ninetyeast  Ridge  as  proposed  by 
Stein  and  Okal  (1978) .  Because  most  of  the  fracture  zones 


along  the  India/Antarctica  boundary  are  poorly  surveyed, 
the  data  used  here  are  inadequate  to  resolve  this  question. 

The  Calculated  Velocity  Diagram 

Table  3  lists  the  relative  plate  velocities  at  the  site 
of  the  Indian  Ocean  triple  junction  calculated  from  the  poles 
derived  here.  When  one  plots  these  as  a  vector  diagram 
superimposed  on  the  observed  velocities  and  uncertainties 
(Figure  5b) ,  one  sees  that  the  predicted  relative  motions 
agree  excellently  with  observations  in  the  neighborhood  of 
the  triple  junction.  The  agreement,  well  within  the 
uncertainties  of  the  observations,  demonstrates  the  consistency 
of  the  triple  junction  data  with  other  measurements  of  the 
relative  motions  among  these  plates  and  validates  the 
assumption  of  rigid  plate  behavior.  We  feel  confident  that 
the  errors  in  the  local  relative  plate  velocities  calculated 
here  are  small,  perhaps  2°  in  direction  and  2  mm/yr  in 
magnitude.  Therefore,  in  the  following  discussion  of  the 
stability  and  evolution  of  the  triple  junction,  we  take  the 
velocity  triangle  of  Figure  5b  to  represent  exactly  the 
instantaneous  relative  plate  velocities  at  the  triple  junction. 
Any  inaccuracy  introduced  by  this  simplifying  assumption  will 


be  small. 


STABILITY  AND  CONFIGURATION 


The  topographic  and  magnetic  data  presented  here  (Figure 
4)  indicate  that  the  Indian  Ocean  triple  junction  is  a  RRR 
type.  Furthermore,  it  appears  to  be  stable  in  the  sense  of 
McKenzie  and  Morgan  (1969) ,  since  its  geometrical  arrangement 
has  persisted  for  a  finite  period  of  geologic  time.  From 
Figure  4,  one  sees  that  the  configuration  of  the  junction 
has  remained  largely  unchanged  for  five  to  ten  million  years, 
and  Sclater  et  al .  (in  preparation)  present  less  detailed 
data  indicating  that  the  triple  junction  has  evolved  in  much 
the  same  manner  since  40  Ma. 

Stability  -  Normal  RRR 

Where  adequate  to  resolve  the  question,  the  sea-floor 
spreading  data  show  that  all  three  ridge  axes  are  spreading 
symmetrically  at  right  angles  to  their  strikes.  This  is 
the  situation  encountered  along  the  major  portion  of  the 
world's  mid-ocean  rift  system.  Assuming  this  to  hold 
rigorously,  one  can  construct  a  stability  diagram  of  the 
triple  junction  (Figure  6a)  after  McKenzie  and  Morgan  (1969) , 
using  the  velocity  triangle  of  Figure  5b. 

Figure  6a  shows  the  junction  to  be  stable,  of  course, 
as  it  must  be  under  the  assumptions.  The  speed  of  each  of 
the  three  plates  relative  to  the  junction  (Jn)  is  34  mm/yr, 
and  the  traces  of  the  junction  on  the  Indian,  Antarctic,  and 
African  plates  have  trends  of  N20°E,  N74°E,  and  N102°E, 
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FIGURE  6 

a.  The  stability  diagram  for  the  triple  junction  based 
on  the  vector  diagram  of  Figure  5b,  assuming  symmetric 
spreading  in  all  three  ridge  axes  in  directions  nor¬ 
mal  to  their  trends.  The  lighter  lines  represent 

the  frames  of  reference  which  are  stationary  with 
respect  to  each  plate  boundary.  They  intersect  at 
the  frame  of  reference  Jn  in  which  the  triple  junc¬ 
tion  is  stationary.  Dashed  lines  represent  the 
motion  of  each  of  the  three  plates  with  respect  to  the 
triple  junction 

b.  The  stability  diagram  as  in  Figure  6a,  but  allowing 
oblique  spreading  on  the  Central  and  Southwest  Indian 
Ridges  to  match  the  observed  trends  of  the  triple 
junction  traces  on  the  African  and  Antarctic  plates. 
The  triple  junction  is  stationary  in  frame  of 
reference  JQ. 


STABILITY  DIAGRAM 

(a)  NORMAL  RRR 

INDIA 


(b)  OBLIQUE  RRR 

INDIA 


ANTARCTICA 
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respectively.  These  traces  are  the  sutures  between  crust 
formed  at  different  plate  boundaries.  One  notes  that  as  the 
triple  junction  evolves,  the  Central  and  Southwest  Indian 
Ridges  grow  longer  by  22  mm/yr  and  33  mm/yr  respectively, 
while  the  SEIR  shortens  at  a  rate  of  15  mm/yr.  The  rapid 
lengthening  of  the  SWIR  compared  to  its  spreading  rate  causes 
the  acute  triangular  morphology  on  that  branch,  with  the 
rough  SWIR-created  crust  enclosed  by  the  converging  triple 
junction  traces  on  the  African  and  Antarctic  plates. 

Formally,  Figure  6a  represents  the  instantaneous  relative 
velocities  precisely  at  the  site  of  the  triple  junction. 

This  is  geophysically  inappropriate,  however,  and  it  is 
more  proper  to  consider  the  Figure  to  represent  an  average 
over  a  relatively  short  period  of  time  and  a  relatively 
small  area,  say  about  one  million  years  and  a  few  tens  of 
kilometers  square.  When  viewed  in  this  manner,  the  stability 
diagram  of  Figure  6a  cannot  alone  account  for  the  data  in  the 
area  of  the  triple  junction.  The  topographic  and  magnetic 
features  interpreted  as  marking  the  boundaries  between  sea 
floor  created  at  the  SWIR  and  sea  floor  created  at  the  CIR 
and  SEIR  trend  roughly  N90°E  and  N60°E  respectively, 

(Figure  4) ,  not  the  N102°E  and  N74°E  shown  on  the  stability 
diagram.  Furthermore,  the  magnetic  anomaly  data  (Figure  4) 
clearly  demonstrate  that  the  length  of  the  SEIR  has  changed 
relatively  little  over  the  last  five  to  ten  million  years, 


A 


certainly  not  by  the  75  to  150  km  decrease  predicted  by  a 
simplistic  application  of  Figure  6a.  One  may  account  for  these 
observations  in  a  number  of  ways,  but  they  first  lead  one 
to  consider  another  quasi-instantaneous  stability  diagram. 
Stability  -  Oblique  RRR 

Based  upon  the  crustal  boundary  trends  on  the  African  and 
Antarctic  plates  and  the  nearly  constant  length  of  the  SEIR, 
one  can  redraw  the  stability  diagram  of  the  Indian  Ocean 
triple  junction  (Figure  6b) ,  imposing  a  stable  condition  on 
the  junction  (Jq)  by  allowing  the  SWIR  and  the  CIR  to  trend 
obliquely  to  their  spreading  directions.  The  geometry  of  sea 
floor  evolution  thus  implied  is  in  excellent  accord  with  the 
general  magnetic  and  topographic  patterns  observed, 
particularly  along  the  Southwest  branch.  The  small  obliquity 
required,  less  than  15°  on  the  SWIR  and  only  10°  on  the  CIR, 
may  well  be  consistent  with  observations  very  close  to  the 
triple  junction. 

One  may  achieve  a  similar  result  by  introducing  asymmetry 
of  spreading  rather  than  obliquity,  but  the  degree  of  asymmetry 
required  is  excessive,  especially  in  the  SWIR,  and  may  be 
ruled  out  by  the  available  magnetic  data.  A  lesser  degree  of 
asymmetry  might  still  be  present  along  with  the  obliquity, 
however.  In  any  case,  the  weight  of  the  evidence  strongly 
suggests  that  the  SEIR  is  spreading  normally  and  symmetrically 


right  up  to  the  triple  junction. 
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Other  Configurations 

Bathymetric  and  magnetic  observations  indicate  that  the 
Indian  Ocean  triple  junction  is  a  RRR  junction,  and  this 
may  be  further  confirmed  by  considerations  of  stability. 

One  cannot  obtain  stable  ridge-fracture-fracture  (RFF)  or 
ridge-ridge-fracture  (RRF)  triple  junction  structures 
consistent  with  the  data. 

In  order  that  a  RFF  triple  junction  be  stable,  the  ridge 
axis  trend  must  pass  through  the  opposite  vertex  of  the 
velocity  triangle.  Reference  to  Figure  6  shows  that  this 
can  be  accomplished  most  easily  by  forcing  the  SWIR  trend  to 
pass  through  the  Indian  vertex,  but  even  this  demands  more 
asymmetry  and/or  obliquity  than  the  data  will  permit. 
Moreover,  there  is  absolutely  no  sign  of  a  fracture  zone  at 
the  triple  junction  on  the  India/Antarctica  plate  boundary, 
which  is  the  best  defined  of  the  three.  Neither  is  there 
one  on  the  Africa/India  boundary.  This  rules  out  the 
possibility  of  a  RFF  junction. 

Neither  can  the  triple  junction  have  a  stable  RRF 
configuration,  since  stability  in  this  case  requires  that  the 
two  ridge  axis  trends  intersect  on  the  third  side  of  the 
velocity  triangle.  Figure  6b  shows  that  the  SWIR  and  CIR 
trends  might  possibly  intersect  on  the  India/Antarctica  side 
of  the  velocity  triangle,  but  in  an  RRF  configuration,  that 
would  imply  a  fracture  zone  on  the  India/Antarctica  plate 


boundary  at  the  triple  junction.  Once  again,  there  is  no 
sign  of  this  on  the  SEIR.  This  rules  out  a  RRF  structure  and 
is  further  evidence  that  the  Indian  Ocean  triple  junction  is 
in  fact  a  stable  ridge-ridge-ridge  junction. 

STRUCTURE  NEAR  THE  TRIPLE  JUNCTION 

It  is  important  to  remember  that  the  stability  diagrams 
of  Figure  6  represent  quasi-instantaneous  relationships 
actually  at  the  triple  junction.  The  pattern  of  sea-floor 
evolution  with  which  we  are  comparing  them,  however,  was 
laid  down  over  a  comparatively  long  period  of  time,  about 
10  million  years,  and  in  part  was  created  quite  far  from 
the  coeval  triple  junction.  During  most  of  this  time  the 
normal  spreading  situation  of  Figure  6a  might  have 
predominated,  provided  that  periodically  the  triple  junction 
jumped  to  the  northwest  along  the  trend  of  the  Southeast 
Indian  Ridge,  creating  fracture  zones  in  the  Central  and 
Southwest  branches.  If  the  frequency  and  magnitude  of  these 
jumps  were  properly  scaled,  the  time-averaged  relationships 
would  resemble  those  of  Figure  6b  and  would  account  equally 
well  for  the  general  pattern  of  triple  junction  evolution. 

In  fact.  Figures  6a  and  6b  represent  end  members  of  a 
continuum  of  configurations,  any  of  which  might  dominate  at 
various  times.  In  order  to  determine  the  true  configuration 
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and  stability  conditions  prevailing  at  the  Indian  Ocean  triple 
junction,  one  must  examine  in  detail  the  structure  of  each  of 
the  three  ridges  as  they  approach  the  junction.  The 
conclusion  is  critically  dependent  on  the  structures  of  the 
Central  and  Southwest  branches  particularly,  and,  unfortunately, 
the  available  data  are  insufficient  to  allow  one  to  determine 
these  structures  to  the  desired  precision. 

The  Southwest  Indian  Ridge 

A  closer  examination  of  the  magnetic  and  topographic 
information  along  the  SWIR  immediately  reveals  that  it  is 
unable  to  answer  the  question  of  spreading  obliquity.  One 
may  easily  interpret  the  data  to  be  consistent  with  either 
extreme  (Figure  7a  and  b) .  In  general,  the  bathymetric 
contours  suggest  oblique  spreading  (Figure  7b)  except  in  a 
section  between  68 °E  and  69 °E  (across  which  runs  the  clearest 
magnetic  profile) .  Also,  there  is  an  east-west  (non-oblique) 
trending  grain  within  20  km  of  the  triple  junction.  The 
resolving  power  of  wide  beam  echo-sounding  records  is  too 
low,  however,  to  define  the  structure  in  adequate  detail. 

Sclater  et  al .  (in  preparation)  show  that  farther  to  the 
west  the  SWIR  trends  perpendicular  to  its  direction  of  spreading. 
This  has  no  certain  implications  nearer  the  triple  junction, 
however.  It  is  possible  that  a  given  section  of  spreading 
center  originates  as  an  oblique  spreading  center  at  the  triple 
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FIGURE  7 


An  interpretation  of  the  Southwest  Indian  spreading 
center  as  a  system  of  en  echelon  segments  trending 
normal  to  their  spreading  direction.  Shaded  areas 
are  shallower  than  3000  m.  Deeper  contours  within 
the  spreading  center  area  are  shown  at  500m  intervals 
Postulated  fracture  zones  are  indicated  by  heavy 
lines,  and  the  postulated  area  of  Bruhnes  age  crust 
is  stippled.  Heavy  dashed  lines  represent  the  boun¬ 
daries  of  crust  created  at  the  Southwest  Indian  Ridge 
An  interpretation  of  the  Southwest  Indian  spreading 
center  as  a  system  of  spreading  segments  trending 
oblique  to  their  direction  of  spreading.  Features 
are  as  in  Figure  7a. 


junction,  and,  as  the  evolving  junction  moves  away  toward  the 
east,  the  section  gradually  reorients  itself  to  lie  perpendicular 
to  its  spreading  direction. 

The  Central  Indian  Ridge 

Although  the  magnetic  anomaly  data  are  better  over  the 
CIR  than  over  the  SWIR,  it  is  again  impossible  to  distinguish 
between  normal  and  slightly  oblique  spreading  near  the  triple 
junction.  Both  bathymetric  contours  and  magnetic  correlations 
suggest  an  oblique  ridge  orientation  in  the  expected  sense  in 
the  20  km  north  of  the  triple  junction  (Figure  8) .  However, 
the  bathymetry  also  suggests  the  presence  of  a  small  transform 
fault  at  25.5°S  (Figure  8)  whose  existence  would  invalidate 
the  magnetic  correlation.  In  the  absence  of  a  more  detailed 
survey,  the  spreading  center  strike  can  be  determined  only 
to  within  a  few  degrees. 

The  Southeast  Indian  Ridge 

Figure  8  demonstrates  that  as  the  SEIR  approaches  the 
triple  junction,  it  is  indeed  spreading  symmetrically  in  a 
direction  normal  to  its  strike.  Moreover,  the  young  magnetic 
anomalies  on  the  Indian  plate  associated  with  the  SEIR 
(anomalies  2  and  2')  extend  farther  to  the  north  than  is 
consistent  with  the  present  length  of  the  SEIR.  This  indicates 
that  the  steady  shortening  of  the  SEIR  inferred  from  Figure  6a 
has  been  proceeding  since  3  Ma.  Magnetic  anomalies  on  the 
Antarctic  plate  (Figure  4),  however,  show  that  this  has  not 
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FIGURE  8 

A  detail  of  the  magnetic  anomalies  along  track  in  the 
vicinity  of  the  triple  junction,  superimposed  on  the 
bathymetric  contours.  Areas  deeper  than  3500  m  in 
the  axial  valleys  are  shaded.  The  magnetic  correlations 
are  discussed  in  the  text. 
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continued  uninterrupted  since  the  time  of  Anomaly  3'  (5  Ma) . 

Some  combination  of  ridge  jumps  and  oblique  spreading  must 
have  intervened. 

EVOLUTION  OF  THE  TRIPLE  JUNCTION  SINCE  10  Ma 

In  spite  of  not  being  able  to  determine  the  precise 
configuration  of  the  triple  junction,  one  can  still  map  out 
the  general  evolution  of  the  junction  since  10  Ma.  One  may 
construct  the  regional  pattern  of  crustal  isochrons  in  the 
following  manner. 

First,  one  assumes  that  the  poles  and  rates  of  the 
instantaneous  relative  rotations  among  the  Indian  Ocean  plates 
determined  here  (Table  3)  have  remained  constant  since  10  Ma. 
This  study  shows  the  assumption  is  justified.  Furthermore, 
Sclater  et  al .  (in  preparation)  have  shown  that  the  poles  in 
question  have  changed  very  little,  if  at  all,  since  40  Ma. 
Though  it  is  geometrically  impossible  that  all  three  poles 
remain  fixed  with  respect  to  their  associated  plates  for  a 
finite  period  of  time,  the  rotations  since  10  Ma  are  small 
enough  that  the  errors  introduced  by  this  simple  scheme  are 
insignificant.  Next,  each  spreading  axis,  as  delineated  by 
the  topographic  and  magnetic  data,  is  rotated  about  its  pole 
of  opening  to  define  the  position  of  the  plate  boundary  at 
various  times  in  the  past.  These  loci  represent  isochrons 
on  the  oceanic  crust. 


133 

Figure  9  displays  the  pattern  of  one  million  year 
isochrons  since  10  Ma  superimposed  on  the  3000  m  bathymetric 
contour.  Symmetric  spreading  at  constant  rate  has  been 
assumed,  except  where  ridge  jumps  occur  on  the  SEIR. 

Comparison  with  Figure  4  shows  this  simple  scheme  to  be  an 
excellent  representation  of  the  tectonic  history  of  the  triple 
junction  over  the  last  ten  million  years.  The  figure  is 
based  upon  a  configuration  of  normally  oriented  ridge  axes. 

Not  shown  is  the  pattern  of  isochrons  produced  by  spreading 
obliquities  near  the  triple  junction.  It  is  evident  from 
Figure  7  that  the  two  patterns  are  nearly  identical  except 
for  details  within  the  domain  of  the  Southwest  branch. 

The  traces  of  the  triple  junction  separate  crust  formed 
at  different  ridges  (Figure  9) .  Those  on  the  African  and 
Antarctic  plates  converge  in  an  acute  angle  at  the  triple 
junction  and  enclose  the  triangle  of  crust  formed  at  the 
SWIR.  The  geometry  is  a  consequence  of  the  rapid  extension 
of  the  SWIR  as  the  evolving  triple  junction  moves  swiftly 
eastward  with  respect  to  Africa  (Figure  10) .  This  is 
responsible  for  the  principal  topographic  features  along 
the  Africa/Antarctica  plate  boundary  between  6"’°E  and  the 
triple  junction.  The  same  process  acting  since  40  Ma  has 
created  the  main  topographic  grain  of  the  SWIR  from  as  far 
westward  as  the  Melville  Fracture  Zone  at  61°E  (Sclater  et 
al . ,  in  preparation) . 
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FIGURE  9 


The  isochron  pattern  produced  by  the  evolving  triple 
junction  of  normally  spreading  ridges  superimposed  on 
the  generalized  bathymetry.  The  isochrons  are  at 
1  m.y.  intervals,  and  their  ages  are  given  by  the  num¬ 
bers.  Crust  formed  since  1  Ma  is  stippled.  Heavy  lines 
represent  fracture  zones.  Heavy  dashed  lines  represent 
traces  of  the  triple  junction.  The  construction  of  the 
Figure  is  discussed  in  the  text. 
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FIGURE  10 

A  schematic  representation  of  the  evolving  Indian  Ocean 

triple  junction  at  times  8  Ma,  6  Ma,  4  Ma,  and  2  Ma. 

The  African  plate  has  been  held  fixed  in  its  present 

position.  The  isochrons  are  at  1  m.y.  intervals,  and 

in  each  reconstruction  crust  of  less  than  1  m.y. 

* 

age  is  stippled. 


I 
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The  triple  junction  traces  on  the  African  and  Antarctic 
plates  are  associated  with  broken  yet  lineated  topographic 
highs,  and  so  to  a  lesser  degree  is  the  trace  on  the  Indian 
plate.  There  may  be  a  causal  relationship  in  this  association, 
in  that  the  highs  may  be  an  effect  of  the  contrast  between 
crust  created  at  different  spreading  rates.  Another  possible 
cause  for  lineated  highs  along  the  triple  junction  traces  is 
that  they  were  built  by  excess  magmatism  at  or  near  the  triple 
junction  in  a  process  similar  to  the  construction  of  aseismic 
ridges.  Phillips  and  Tapscott  (in  preparation)  postulate  a 
similar  origin  for  the  Yermak  Plateau  and  the  Morris  Jessup 
Rise  in  the  Arctic  Ocean. 

SUMMARY 

We  have  located  the  Indian  Ocean  triple  junction  to  within 
5  km  and  have  shown  it  to  be  a  ridge-ridge-ridge  junction.  We 
have  succeeded  in  determining  the  instantaneous  relative 
angular  velocity  vectors  of  the  Indian  Ocean  plates  and  have 
calculated  the  velocity  triangle  of  the  triple  junction. 
Unfortunately,  our  inability  to  determine  the  detailed 
structure  of  the  Southwest  Indian  Ridge  prevents  us  from 
determining  the  true  stability  conditions  of  the  junction. 

Both  postulated  configurations  have  probably  prevailed  at 


different  times. 


Regardless  of  the  precise  structure  of  the  triple  junction 
we  have  been  able  to  map  its  evolution  and  have  demonstrated 
that  the  morphology  of  the  Africa/Antarctica  plate  boundary 
in  the  region  is  a  direct  consequence  of  that  evolution. 
Finally,  we  have  pointed  out  the  intriguing  possibility  that 
the  triple  junction  traces  on  the  three  plates  may  have 
topographic  expression. 
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CHAPTER  4 

THE  FINITE  ROTATIONS  OF  A  THREE  PLATE  SYSTEM 


145 


THE  FINITE  ROTATIONS  OF  A  THREE  PLATE  SYSTEM 
WITH  APPLICATIONS  TO  THE  LABRADOR  SEA 


C.  Tapscott 

Department  of  Earth  and  Planetary  Sciences 
Massachusetts  Institute  of  Technology 
Cambridge,  MA  02139  and 
Woods  Hole  Oceanographic  Institution 
Woods  Hole,  MA  02543 


146 


ABSTRACT 

A  method  is  presented  for  determining  the  best  recon¬ 
struction  of  the  past  relative  positions  of  three  tectonic  plates 
around  a  triple  junction  according  to  a  fitting  criterion 
applied  to  sea-floor  spreading  data.  The  method  enables  one 
to  set  confidence  limits  on  the  optimal  set  of  internally  con¬ 
sistent  finite  rotations  based  upon  the  probable  errors  in  the 
data.  The  technique  is  illustrated  by  reconstructions  of  the 
plates  about  the  now  extinct  Labrador  Sea  triple  junction  at 
the  times  of  anomalies  21  (50  Ma)  and  24  (56  Ma) . 


INTRODUCTION 


Over  the  last  several  million  years  the  relative  motions 
of  the  lithospheric  plates  have  been  the  principal  mechanism 
shaping  the  face  of  the  earth.  For  this  reason,  one  of  the  major 
thrusts  of  geophysical  research  since  the  advent  of  the  new 
global  tectonics  has  been  toward  determining  as  accurately 
as  possible  the  relative  positions  of  the  plates  at  various  times 
in  the  past.  Early  in  this  endeavor,  Bullard  et  al.  (1965) 
recognized  that,  by  a  theorem  due  to  Euler,  such  plate  recon¬ 
structions  could  be  described  in  terms  of  finite  rotations  about 
axes  passing  through  the  center  of  the  earth.  The  problem  of 
reconstructing  the  past  relative  positions  of  two  plates  is 
then  the  problem  of  discovering  the  finite  rotation  that  best 
brings  into  coincidence  two  patterns,  one  on  each  plate,  that 
represent  the  adjacent  edges  of  the  plates  at  the  time  for  which 
the  reconstruction  is  valid.  Many  such  reconstructions  have 
been  made  fitting  together  the  margins  of  the  Gondwana  continents 
(e.g.,  Bullard  et  al . ,  1965;  Smith  and  Hallam,  1970)  and  magnetic 
anomaly  isochron  lineations  and  fracture  zone  traces  in  the 
ocean  basins  (e.g.,  McKenzie  and  Sclater,  1971;  Pitman  and 
Talwani,  1972) . 

No  two  plates  exist  in  isolation.  All  pairs  are  part  of 
the  global  system,  and  all  the  pairwise  relative  rotations  form 
an  internally  consistent  set  if  the  plates  are  indeed  rigid 
bodies.  All  closed  cyclic  products  of  relative  rotations  must 
result  in  unity.  The  interdependence  of  the  finite  rotations 
provides  additional  constraints  which  can  be  an  advantage  in 
seeking  the  best  reconstruction,  in  that  a  paucity  of  data 
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between  two  plates  can  sometimes  be  overcome  by  enlarging  the 
system  to  include  three  or  more  plates  around  one  or  more 
triple  junctions.  With  adequate  data,  one  may  test  the  validity 
of  the  rigid  plate  hypothesis  for  the  plates  involved  in  such 
a  multiple  plate  system,  questioning  it  if  no  acceptable  recon¬ 
struction  can  be  achieved. 

This  paper  extends  the  fitting  technique  of  Hellinger 
(1979)  to  deal  with  three  plates  around  a  triple  junction.  The 
method  presented  here  enables  one  to  find,  with  a  minimum  of 
effort,  the  internally  consistent  set  of  finite  rotations 
describing  the  best  possible  reconstruction  of  a  three  plate 
system  according  to  the  fitting  criterion.  Perhaps  most  impor¬ 
tant,  Hellinger* s  (1979)  technique  allows  one  to  set  statisti¬ 
cally  meaningful  limits  on  the  region  of  valid  rotation  tensors. 
The  three  plate  method  is  illustrated  by  reconstructions  of  the 
plates  around  the  Labrador  Sea  triple  junction  at  the  times  of 
Anomaly  24  (56  Ma)  and  Anomaly  21  (50  Ma) . 

THE  METHOD 

The  best  fit  of  patterns  depends  on  the  statistical 
criterion  defining  the  quality  of  the  fit.  Criteria  commonly 
used  include  the  simple  but  effective  method  of  "eyeballing" 
the  patterns,  variations  on  minimizing  the  sum  of  the  squares 
of  the  distances  of  each  point  from  the  other  lineation  (Bullard 
et  al . ,  1965;  Pilger,  1978),  and  minimizing  the  overlap  and 
underlap  area  between  the  patterns  (McKenzie  and  Sclater,  1971) . 
These  techniques  have  both  advantages  and  disadvantages  relative 


to  one  another.  All  have  the  common  disadvantage  that  they 
fail  to  take  account  of  the  known  uncertainties  in  the  data. 
Hellinger  (1979)  has  recently  introduced  a  method  correcting 
this  shortcoming,  and,  moreover,  Hellinger' s  technique  uses 
the  data  and  their  uncertainties  to  set  statistically  meaning¬ 
ful  limits  on  the  region  of  valid  finite  rotations. 

The  technique  is  particularly  well-suited  to  the  type 
of  data  usually  encountered  in  sea-floor  spreading  studies.  The 
data  must  satisfy  the  following  conditions  or  assumptions 
(Hellinger,  1979) : 

1.  The  patterns  or  lineations  to  be  matched  on  the  two 

plates  are  known  only  at  a  finite  number  of  discrete 

points  (shiptrack  crossings  of  magnetic  anomalies  or 
fracture  zones) . 

2.  The  location  of  each  datum  is  uncertain  (to  the  limit 
of  resolution  in  interpretation)  and  the  distribution 
of  errors  may  be  adequately  described  by  a  bivariate- 
normal  probability  density  function.  These  errors  are 
uncorrelated. 

3.  The  data  on  each  plate  may  be  divided  into  a  number  of 

sections  (e.g.,  a  fracture  zone  trace  or  a  magnetic 

anomaly  lineation  between  adjacent  fracture  zones) , 
each  of  which  may  be  well-represented  by  a  segment  of 
a  great  circle. 

4.  One  may  identify  a  one-to-one  correspondence  between 
the  sections  of  data  on  one  plate  and  those  on  the 
other  plate. 
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The  simplest  and  most  readily  calculated  criterion  of 
fit  applicable  to  the  method  (criterion  4>0  in  Hellinger,  1979) 
is  easily  described.  After  a  trial  rotation,  the  best  fit 
great  circle  is  calculated  for  each  section  of  data  (now  the 
united  set  of  data  from  both  plates,  see  Figure  1) ,  weighting 
each  point  according  to  its  uncertainty.  The  proper  great 
circle  is  that  which  minimizes  the  sum  of  the  squares  of  the 
weighted  distances  of  the  points  from  the  plane  of  the  great 
circle.  The  measure  of  the  misfit  of  each  section  is  just 
this  sum  of  squared,  weighted  distances,  and  the  measure  of 
misfit  of  the  reconstruction  is  the  sum  of  the  measures  of  the 
individual  sections.  That  is,  if  there  are  n^j  points  in 
section  j  on  plate  i,  then  the  measure  of  misfit  of  that  section 
for  a  rotation  r  is 

4oj  -  3i(  ersijk  )2^jk  +  3i<Pj'r!i2jk|2o23*  '  (1) 


where  is  the  unit  vector  perpendicular  to  the  plane  of  the 

best  fit  great  circle  of  section  j,  the  §. ,,  are  the  unit  vectors 

13k 

through  the  data  points  of  section  j  on  plate  i,  and  the  a 

ljk 

are  the  standard  deviations  of  their  error  distributions,  scaled 
to  a  unit  sphere.  If  there  are  a  total  of  J  sections,  then  the 
measure  of  misfit  for  the  reconstruction  defined  by  the  rotation 
r  is 

*0  "  .  =  *oj  •  (2) 

D=1 


The  best  reconstruction  is  the  one  which  minimizes  this  measure 
of  misfit.  The  method  makes  minimal  assumptions  about  the 
actual  shape  of  the  ancient  plate  boundary  and  is  well-suited 


FIGURE  1 


A  schematic  diagram  illustrating  the  method  of  pattern  fitting. 
The  data  on  one  plate  are  shown  as  circles  and  on  the  other 
plate  as  squares,  and  corresponding  sections  of  data  are 
indicated  before  rotation.  After  a  trial  rotation,  great 
circles  are  fit  to  each  section  of  plate  boundary,  using  the 
data  from  both  plates.  After  Hellinger  (1979) . 
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to  the  sparse  data  sets  often  encountered  in  investigations  of 
sea-floor  spreading. 

The  measure  of  misfit  of  a  reconstruction  of  several 

plates  is  merely  the  sum  of  the  measures  of  misfit  of  each 

plate  pair.  This  paper  presents  a  computer  program  (Appendix) 

which  calculates  the  measure  of  misfit  for  a  reconstruction 

of  a  three  plate  system  and  searches  for  the  interdependent  set 

of  finite  relative  rotations  describing  the  best  reconstruction 

according  to  criterion  4>  (Hellinger,  1979)  .  This  technique 

o 

provides  a  valuable  tool  for  reconstructing  the  past  relative 
positions  of  several  tectonic  plates  and  for  defining  the 
confidence  limits  we  may  set  on  our  knowledge  of  their  relative 
motions . 

THE  LABRADOR  SEA  -  ANOMALY  24 

Srivastava  (1978)  has  presented  the  results  of  a  geo¬ 
physical  survey  in  the  Labrador  Sea  and  has  identified  several 
marine  magnetic  anomaly  lineations, giving  a  set  of  finite 
rotations  which  restore  the  Greenland,  North  American,  and 
Rockall  plates  to  their  relative  positions  at  the  time  of 
Anomaly  24  (56  Ma) .  Phillips  and  Tapscott  (in  preparation) 
also  present  a  set  of  finite  rotations  for  the  same  purpose, 
basing  their  study  solely  on  published  data  from  the  North 
Atlantic  Ocean  (Laughton,  1971;  Johnson  et  al^. ,  1972;  Johnson 
and  Egloff,  1973;  Vogt  and  Avery,  1974;  Kristof fersen  and 
Talwani,  1977). 

After  compiling  a  data  set  of  locations  of  Anomaly  24 


mm 
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from  these  sources  (Table  1,  Figure  2),  one  may  apply  the 
criterion  of  fit  described  above  to  both  sets  of  rotations 
(Table  2) .  Overall,  the  rotations  of  Srivastava  (1978)  yield 
a  much  smaller  measure  of  misfit,  though  those  of  Phillips  and 
Tapscott  (in  preparation)  produce  a  slightly  better  fit  for  the 
Greenland-Rockall  and  North  America-Rockall  boundaries.  This 
is  obvious  when  one  plots  the  reconstructions  (Figures  3  and  4) . 
The  Phillips  and  Tapscott  (in  preparation)  rotations  take  North 
America  much  too  close  to  Greenland,  while  the  Srivastava  (1978) 
rotations  move  Rockall  slightly  too  far  toward  the  west.  By 
arbitrarily  fixing  the  Greenland-North  America  and  Greenland- 
Rockall  poles  and  making  slight  adjustments  in  the  angles  of 
rotation,  one  may  significantly  improve  both  reconstructions, 
though  that  derived  from  Srivastava  (1978)  is  clearly  superior. 

Starting  from  this  point  and  searching  for  the  best  set 
of  finite  rotations  (down  to  a  resolution  of  0.25°  in  pole 
location  and  0.01°  in  rotation  angle),  one  finds  a  minimum  in 
the  measure  of  misfit  for  a  set  including  a  Greenland-North 
America  pole  close  to  that  of  Srivastava  (1978) .  The  other  two 
rotation  poles,  however,  have  moved  toward  the  northwest, 
giving  a  much  better  fit  of  Rockall  to  the  Greenland-North 
America  combination  (Table  2,  Figure  5).  The  set  of  finite 
rotations  derived  here  clearly  yields  an  improved  reconstruction 
about  the  Labrador  Sea  triple  junction  for  Anomaly  24  time. 

Having  discovered  the  optimal  set  of  finite  rotations  for 
the  data,  one  may  determine  the  confidence  limits  on  the  finite 
rotation  poles.  One  method  of  doing  this  is  to  apply  the 


TABLE  1.  Anomaly  24  Locations  Used  in  this  Study 
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♦From  Phillips  and  Tapscott  (in  preparation).  All  other  data  are  from  Srivastava  (1978). 
All  data  are  estimated  to  be  accurate  to  within  10  km. 
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FIGURE  2 


Locations  of  Anomaly  24  along  the  margins  of  Greenland  (circles) , 
North  America  (squares) ,  and  Rockall  Bank  (stars)  used  in  this 
study.  Corresponding  sections  of  data  are  indicated  by  joined 
balloons.  The  positions  of  the  continental  masses  are  indicated 
by  shading  the  regions  within  the  1000  fathom  bathymetric 
contour.  The  data  and  sources  are  listed  in  Table  1. 
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TABLE  2.  Comparison  of  the  Quality  of  Reconstructions  for  Anomaly  24. 


Source 

Plate 

Pair* 

Finite  Rotation 

N.  Lat.  E.  Long.  Angle 

Relative 

Measure 

Total 

Relative 

Measure 

Srivastava 

N/G 

26.0 

1.5 

-2.80 

5.23 

(1978) 

R/G 

50.2 

130.9 

-10.97 

18.9 

91.2 

R/N 

40.0 

145.0 

-11.40 

67.0 

Phillips  and 

N/G 

-9.1 

25.6 

-3.22 

542 

Taps co tt 

R/G 

35.6 

127.9 

-9.86 

13.0 

584 

(in  prep.) 

R/N 

35.4 

147.8 

-11.15 

29.1 

Revised  after 

N/G 

26.0 

1.5 

-2.80 

5.23 

Srivastava 

R/G 

50.2 

130.9 

-10.77 

13.9 

46.7 

(1978) 

R/N 

39.8 

145.2 

-11.20 

27.6 

Revised  after 

N/G 

-9.1 

25.6 

-2.37 

100 

Phillips  and 

R/G 

35.6 

127.9 

-10.11 

39.0 

180 

Tapscott 

R/N 

35.8 

142.7 

-10.97 

41.3 

(in  prep.) 

This  paper 

N/G 

22.0 

4.0 

-2.62 

5.30 

R/G 

63.6 

120.6 

-12.01 

13.7 

35.2 

R/N 

56.2 

141.0 

-11.90 

16.2 

This  paper 

N/G 

23.5 

2.5 

-2.69 

5.00 

— 

(two  plates) 

*N/G  signifies  North  American  relative  to  Greenland; 
R/G,  Rockall  relative  to  Greenland; 

R/N,  Rockall  relative  to  North  America. 


FIGURE  3 


A  reconstruction  of  Anomaly  24  about  the  Labrador  Sea 
junction  according  to  Srivastava's  (1978)  rotations, 
are  as  in  Figure  2. 


triple 

Symbols 
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FIGURE  4 


A  reconstruction  of  Anomaly  24  about  the  Labrador  Sea  triple 
junction  according  to  rotations  of  Phillips  and  Tapscott  (in 
preparation).  Symbols  are  as  in  Figure  2. 
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FIGURE  5 


A  reconstruction  of  Anomaly  24  about  the  Labrador  Sea  triple 
junction  using  the  rotations  derived  in  this  paper.  Symbols 
are  as  in  Figure  2. 


rigorous  and  exhaustive  technique  of  Hellinger  (1979)  to  each 
plate  boundary  in  turn.  A  somewhat  more  qualitative  but  no 
less  useful  determination  of  the  regions  of  confidence  may 
easily  be  made  by  exploring  a  few  rotations  in  the  neighborhood 
of  the  best-fit  set.  In  the  present  study,  for  example,  the 
testing  of  a  few  poles  in  the  neighborhood  of  the  optimal 
Greenland-North  America  pole  shows  that  when  the  best  measure 
of  misfit  obtainable  for  a  pole  is  greater  than  about  1.25 
times  the  measure  for  the  optimal  set  of  rotations,  the  fit  has 
degraded  to  the  limit  of  acceptability.  Such  degradation  is 
generally  concentrated  in  the  fit  along  one  plate  boundary  (in 
this  case  the  Greenland-North  America  plate  boundary)  and  is 
usually  due  to  poor  fits  of  one  or  two  sections  of  data. 

Figure  6  displays  the  contour  of  the  measure  at  the 
critical  value  about  the  preferred  Greenland-North  America  pole. 
The  region  of  confidence  of  this  pole  is  a  long,  narrow  area 
roughly  joining  the  poles  of  Srivastava  (1978)  and  Phillips  and 
Tapscott  (in  preparation) .  This  is  not  unexpected,  as  poles  of 
rotation  are  usually  well-constrained  in  a  direction  perpendi¬ 
cular  to  the  general  trend  of  the  plate  boundary,  and  poorly 
constrained  along  that  trend.  In  this  case,  because  the  data 
are  confined  to  a  small  area,  one  is  unable  to  rule  out 
absolutely  any  Greenland-North  America  finite  poles  along  a 
50°  arc.  If  it  is  considered  in  combination  with  Greenland- 
Rockall  and  North  America-Rockall  rotations  close  to  those  of 
the  optimal  set,  the  Greenland-North  America  pole  of  Srivastava 
(1978)  is  not  significantly  different  from  that  derived  here. 
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FIGURE  6 


Regions  of  confidence  for  the  Greenland-North  America  finite 
rotation  poles  for  Anomaly  24.  The  continuous  contour  denotes 
the  region  of  confidence  for  the  pole  derived  in  this  paper  for 
the  three  plate  system,  and  the  dashed  contour  is  that  for  the 
two  plates  alone.  The  Greenland-North  America  finite  rotation 
poles  for  Anomaly  24  discussed  in  the  text  are  also  shown. 
Continental  coastlines  are  shown  for  scale  and  orientation. 
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It  is  of  interest  to  compare  these  results  to  those 
following  from  the  consideration  of  a  system  of  only  two  plates. 
Application  of  Hellinger's  method  shows  that  a  rotation  of  -2.69° 
about  a  pole  at  23.5°N,  2.5°E  best  reconstructs  the  data  along 
the  Greenland-North  American  plate  boundary  without  regard  to 
the  fit  of  the  Rockall  plate  (Table  2) .  This  pole  lies  between 
that  of  Srivastava  (1978)  and  the  pole  derived  in  the  present 
study  of  the  three  plate  system.  Contouring  the  region  of 
confidence  in  a  similar  manner,  one  sees  that  it  closely  resembles 
that  for  the  three  plate  system  (Figure  6) .  The  two  regions 
trend  somewhat  differently,  however,  and  the  confidence  region 
for  the  three  plate  system  is  more  restricted  in  latitudinal 
extent.  These  differences  result  from  the  additional  constraints 
imposed  by  the  other  two  plate  boundaries.  In  this  example, 
these  constraints  are  relatively  weak,  yet  their  effects  are 
significant  in  limiting  the  uncertainty  in  the  Greenland-North 
America  pole  of  rotation. 

THE  LABRADOR  SEA  -  ANOMALY  21 

Locations  of  Anomaly  21  in  the  Labrador  Sea  and  adjoining 
regions  may  be  compiled  from  the  same  sources  as  above  (Table  3, 
Figure  7).  Testing  Srivastava's  (1978)  reconstruction  for 
Anomaly  21  time  (Table  4),  one  finds  that  once  again  these 
rotations  move  Rockall  Bank  too  far  westward  (Figure  8)  and 
that  a  slight  adjustment  of  the  Greenland-Rockall  rotation  angle 
significantly  improves  the  fit  (Table  4) .  Simultaneous  consi¬ 
deration  of  the  three  plate  boundaries  leads  to  an  optimal  set 


TABLE  3.  Anomaly  21  Locations  Used  in  this  Study. 
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FIGURE  7 


Locations  of  Anomaly  21  about  the  Labrador  Sea  triple  junction. 
Symbols  are  as  in  Figure  2. 


TABLE  4.  Comparison  of  the  Quality  of  Reconstructions  for  Anomaly  21 


Source 

Plate 

Pair* 

Finite  Rotation 

N.  Lat.  E.  Long.  Angle 

Relative 

Measure 

Total 

Relative 

Measure 

Srivastava 

N/G 

13.0 

2.0 

-1.10 

3.13 

(1978) 

R/G 

55.0 

135.0 

-9.59 

15.2 

35.5 

R/N 

50.8 

142.8 

-9.87 

17.1 

Revised  after 

N/G 

13.0 

2.0 

-1.10 

3.13 

Srivastava 

R/G 

55.0 

135.0 

-9.44 

7.02 

10.5 

(1978) 

R/N 

50.7 

142.9 

-9.72 

0.341 

This  paper 

N/G 

16.8 

1.3 

-1.15 

2.75 

R/G 

65.5 

127.6 

-10.41 

2.10 

5.42 

R/N 

61.5 

138.6 

-10.44 

0.574 

This  paper 

N/G 

-18.6 

20.0 

-0.95 

2.24 

(two  plates) 

*  N/G  signifies  North  America  relative  to  Greenland; 
R/G,  Rockall  relative  to  Greenland; 

R/N,  Rockall  relative  to  North  America. 
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FIGURE  8 

A  reconstruction  of  Anomaly  21  about  the  Labrador  Sea  triple 
junction  using  rotations  of  Srivastava  (1978)  .  Symbols  are 
as  in  Figure  2. 
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of  finite  rotations  with  a  Greenland-North  America  pole  close 
to  that  of  Srivastava  (Table  4) .  The  other  rotation  poles  have 
moved  toward  the  northwest,  as  was  the  case  for  Anomaly  24. 

This  yields  a  better  reconstruction  (Figure  9) . 

Figure  10  displays  the  uncertainty  region  of  the  Greenland- 
North  America  pole  for  the  Anomaly  21  reconstruction.  It  is 
similar  to  that  for  Anomaly  24  (Figure  6) ,  but  shows  two 
interesting  facts.  First,  though  the  pole  of  Srivastava  (1978) 
is  within  the  acceptable  region,  it  is  close  to  its  limits 
rather  than  along  the  axis  of  the  ellipse.  The  difference 
between  the  two  poles  is  more  significant  than  is  the  case  for 
the  Anomaly  24  reconstructions.  Second,  and  more  striking,  is 
the  vast  difference  between  the  optimal  rotations  for  the  three 
plate  system  and  the  two  plate  system.  The  finite  rotation 
which  best  fits  Anomaly  21  in  the  Labrador  Sea  without  regard 
to  any  data  in  the  North  Atlantic  lies  40°  away  from  the  best 
pole  for  the  three  plate  system  (Table  4,  Figure  10).  Moreover, 
it  is  far  more  poorly  constrained.  This  dramatically  demonstrates 
the  importance  of  using  systems  of  several  plates  whenever 
possible  to  determine  finite  rotations. 

DISCUSSION  AND  CAUTIONS 

The  examples  discussed  here  demonstrate  that  the  method 
presented  in  this  paper  provides  a  useful  and  powerful  tool 
for  determining  the  finite  rotations  necessary  to  best  recon¬ 
struct  multiple  plate  systems  in  plate  tectonics.  However,  in 
using  the  computer  program  (Appendix) ,  it  is  well  to  be  wary  of 


FIGURE  9 


A  reconstruction  of  Anomaly  21  about  the  Labrador  Sea  triple 
junction  using  the  rotations  derived  in  this  paper.  Symbols 
are  as  in  Figure  2. 
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FIGURE  10 


Regions  of  confidence  for  the  Greenland-North  America  finite 
rotation  poles  for  Anomaly  21.  The  continuous  contour  denotes 
the  region  of  confidence  for  the  pole  derived  in  this  paper 
for  the  three  plate  system,  and  the  dashed  contour  is  that  for 
the  two  plates  alone.  The  Greenland-North  America  finite 
rotation  poles  for  Anomaly  21  discussed  in  the  text  are  also 
shown.  Continental  coastlines  are  shown  for  scale  and 
orientation. 


several  potential  pitfalls. 

In  general,  the  measure  of  misfit  varies  smoothly  when 
considered  as  a  function  of  its  six  independent  variables. 

Under  some  configurations  of  data,  however,  the  measure  may 
vary  quite  rapidly  as  a  function  of  one  or  more  of  the  variables, 
especially  the  rotation  angles.  When  this  is  the  case,  the 
measure  must  be  examined  at  suitably  fine  intervals  of  the 
variables.  Moreover,  the  measure  of  misfit  may  have  several 
local  minima  for  some  configurations  of  data.  A  survey  of 
several  poles  can  expose  the  problem  if  it  exists.  Finally,  one 
must  beware  of  a  situation  in  which  one  of  the  finite  rotations 
is  small  compared  to  the  other  two, in  the  sense  that  it  is  the 
difference  between  two  rotations  by  nearly  the  same  angle  about 
nearly  the  same  pole.  As  implied  in  the  discussion  of  Hellinger 
(1979) ,  this  situation  is  unstable  because  a  small  perturbation 
in  either  of  the  two  similar  rotations  can  lead  to  a  large 
change  in  the  third  rotation.  Such  a  condition  can  also  be 
an  advantage,  however,  if  the  smaller  rotation  can  be  well- 
determined  in  its  own  right.  This  in  turn  can  constrain  the 
two  larger  rotations. 

All  the  pitfalls  may  be  safely  avoided  most  easily  if  one 
first  reconstructs  each  individual  plate  boundary  and  under¬ 
stands  the  problems  associated  with  the  fits  of  the  separate 
plate  pairs.  The  method  presented  in  this  paper,  together 
with  the  two-plate  finite  rotation  method  of  Hellinger  (1979) , 
enables  one  to  discover  reliable  finite  rotations  to  reconstruct 
the  ancient  positions  of  several  plates,  and  to  place  limits  on 


their  reliability.  The  method  is  applicable  to  a  wide  class 
of  sea-floor  spreading  data  and  should  prove  generally  useful 
in  refining  our  knowledge  of  the  tectonic  history  of  the  earth. 
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PROGRAM  PLATE3 

****************************************************************** 

PROGRAM  PLATE3  JANUARY  1979  C.TAPSCOTT 

THE  PROGRAM  ALLOWS  ONE  TO  CALCULATE  THE  MEASURE  OF  FIT  OF  SEA¬ 
FLOOR  SPREADING  DATA  FOR  A  RECONSTRUCTION  OF  THREE  PLATES  AROUND  A 
TRIPLE  JUNCTION,  AND  TO  SEEK  THE  SET  OF  FINITE  ROTATIONS  DEFINING 
THE  BEST-FIT  RECONSTRUCTION. 

THE  MEASURE  OF  FIT  OF  EACH  PLATE  PAIR  IS  CALCULATED  ACCORDING  TO 
THE  FITTING  CRITERION  PHI-0  OF  HELL INGER  11979).  THE  TOTAL  MEASURE 
FOR  THE  RECONSTRUCTION  IS  THE  SUM  OF  THE  MEASURES  OF  THE  PLATE 
PAIRS. 

ONE  SPECIFIES  TWO  OF  THE  INITIAL  FINITE  RELATIVE  ROTATIONS.  THE 
THIRD  IS  CALCULATED  BY  REQUIRING  THE  CYCLIC  PRODUCT  TO  BE  UNITY. 

DATA  ON  EACH  PLATE  ALONG  EACH  PLATE  BOUNDARY  ARE  SUBDIVIDED  INTO 
CORRESPONDING  SECTIONS,  EACH  OF  WHICH  MAY  BE  WELL-APPROXIMATED  BY 
A  GREAT  CIRCLE  1  HELL  I NGE R,  1979). 

EACH  DATUM  HAS  ASSOCIATED  WITH  IT  AN  UNCERTAINTY  IN  LOCATION, 
DEFINED  BY  A  STANDARD  DEVIATION  EXPRESSED  IN  KILOMETERS.  IN  THE 
CALCULATION  OF  THE  MEASURE*  EACH  DATUM  IS  WEIGHTED  BY  THE  INVERSE 
OF  THE  SQUARE  OF  ITS  STANOARD  DEVIATION  OF  ERROR. 

REFERENCE: 

HELLINGER,  S.,  THE  STATISTICS  OF  FINITE  ROTATIONS  IN  PLATE 

TECTONICS,  PH.  D.  THESIS,  MASS.  INST.  OF  TECHNOL.,  1979. 

****************************************************************** 


OPERAT ION: 

ALL  INPUT  ON  DEVICE  1. 

ALL  OUTPUT  ON  DEVICE  108  CLINE  PRINTER). 

ALL  INPUT  IS  IN  FREE-FIELD  FORMAT  CG  FORMAT). 


INPUT: 

READ  INSECT (I), 1-1,3) *WF ACTR , 1  PUT  IN  15G)  FORMAT. 

N SECT ( I )  THE  NUMBER  OF  SECTIONS  FOR  PLATE  BOUNDARY  I. 

AT  PRESENT,  THE  PROGRAM  IS  LIMITED  TO  NO  MORE 
THAN  TEN  SECTIONS  PER  BOUNDARY. 


o  n  n 
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C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


WFACTR  ARBITRARY  SCALING  FACTOR  TO  KEEP  THE  MAGNITUDE 

OF  THE  VARIANCES  NEAR  1.  DEFAULT  VALUE  =  100. 

I  PUT  IF  I  PUT  =  0,  INPUT  DATA  POINTS  ARE  LISTED. 

OTHERWISE,  POINTS  ARE  NOT  LISTED.  DEFAULT  =  0. 

READ  DATA  POINTS  ( AL AT  III , ALON (I » , AERR £  I ) , 1= 1 , 5 >  IN  ( 15G  I  FORMAT. 

LATITUDE  AND  LONGITUOE  IN  DECIMAL  DEGREES  {POSITIVE  NORTH  AND 
EAST),  ANO  ERROR  IN  KILOMETERS. 


THE  FIRST  RECORD  IS  POINTS 

CF 

BOUNDARY 

I, 

PLATE 

1. 

SECTION 

I. 

THE  SECOND  IS 
ETC. 

POINTS  OF 

BOUNDARY 

1  , 

PLATE 

1, 

SECTION 

2. 

THEN 

BOUNDARY 

1. 

PLATE 

2, 

SECTION 

1. 

BOUNDARY 

1. 

PLATE 

2, 

SECTION 

2. 

ETC. 

THEN  THE  SAME 

FOR  BOUNDARY 

2, 

AND  THEN 

BOUNDARY 

3. 

AT  PRESENT,  THE  PROGRAM  IS  LIMITED  TO  NO  MORE  THAN  FIVE 
POINTS  PER  SECTION  FOR  EACH  PLATE  {AND  EACH  BOUNDARY ) . 

IF  A  GIVEN  SECTION  FOR  A  GIVEN  PLATE  AND  BOUNDARY  CONTAINS 
FEWER  THAN  FIVE  POINTS,  THE  STRING  MUST  BE  TERMINATED  BY  A 
DUMMY  POINT  WITH  A  LATITUDE  GREATER  THAN  90.  OR  A  ZERO  ERROR. 

READ  THE  OPTION  NUMBER  IN  ( 1 G >  FORMAT. 

SEE  BELOW  FOR  DESCRIPTIONS  OF  THE  OPTIONS. 

READ  THE  INITIAL  ROTATIONS  FOR  BOUNDARIES  1  AND  2  IN  16G)  FORMAT. 

LATITUDE,  LONGITUOE,  ANGLE,  LATITUDE,  LONGITUDE,  ANGLE. 

THE  PROGRAM  WILL  COMPUTE  THE  RESULTANT  THIRD  ROTATION. 


OTHER  INPUT  FOLLOWS  FOR  SOME  OF  THE  OPTIONS. 

OPTION  1. 

FOR  THE  GIVEN  SET  CF  ROTATIONS,  COMPUTE  AND  OUTPUT  THE  MEASURE  FOR 
EACH  OF  THE  BOUNDARIES,  AND  THE  TOTAL  MEASURE. 


OPT  I CN  -1. 

THE  SAME  AS  OPTION  l,  EXCEPT  THAT  A  DECOMPOSITION  OF  THE  MEASURE 


m 
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IS  ALSU  OUTPUT. 

THE  DECOMPOSITION  CONSISTS  OF  THE  EIGENVECTOR  AND  MEASURE  OF  EACH 
SECTION,  AND  THE  DEVIATION  OF  EACH  POINT  [RELATIVE  TO  ITS  ERROR) 
FROM  THE  BEST-FIT  GREAT  CIRCLE  FCR  ITS  SECTION. 


OPTION  2. 

KEEPING  THE  FIRST  ROTATION  CONSTANT  AND  FIXING  THE  POLE  OF  THE 
SECOND  ROTATION,  FIND  THE  ANGLE  FOR  THE  SECOND  ROTATION  THAT 
MINIMIZES  THE  MEASURE. 

INPUT: 

READ  THE  LARGER  AND  SMALLER  INCREMENTS  OF  THE  ANGLE 
SEARCH  AND  THE  HALF-RANGE  Or  THE  ANGLE  SEARCH  IN  (3G) 
FORMAT. 


OPTION  -2. 

THE  SAME  AS  OPTION  2,  EXCEPT  THAT  A  DECOMPOSITION  OF  THE  MEASURE 
IS  GIVEN  FOR  THE  RESULT. 


OPTION  3. 

FIXING  THE  POLES  OF  THE  FIRST  AND  SECOND  ROTATIONS,  SEARCH  FOR  THE 
COMBINATION  OF  ROTATION  ANGLES  THAT  MINIMIZES  THE  MEASURE. 

INPUT: 


AS  IN  OPTION  2. 


OPTION  -3. 

THE  SAME  AS  OPTION  3,  EXCEPT  THAT  A  DECOMPOSITION  OF  THE  MEASURE 
IS  GIVEN  FOR  THE  RESULT. 


OPTION  4. 

FIXING  THE  POLE  OF  THE  FIRST  ROTATION,  SEARCH  FOR  THE  COMBINATION 
OF  FIRST  ROTATION  ANGLE  AND  SECOND  ROTATION  POLE  AND  ANGLE  THAT 
MINIMIZES  THE  MEASURE. 

THE  PROGRAM  SEARCHES  A  HEXAGONAL  GRID  OF  SPECIFIED  RADIUS  ABOUT 
THE  STARTING  POLE,  RECENTERS  THE  GRID  ON  THE  POLE  GIVING  THE  BEST 
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c 
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c 

c 

c 

c 

c 

c 


MEASURE »  AND  REPEATS.  IF  THE  BEST  POLE  IS  IN  THE  CENTER  OF  THE 
GRID,  THE  GRID  IS  REDUCED  BY  THE  SPECIFIED  FACTOR.  THE  SEARCH 
PROCEEDS  UNTIL  THE  GRID  IS  SMALLER  THAN  THE  SPECIFIED  LIMIT. 

STEP  NUMBER  AND  BEST  RESULT  ARE  OUTPUT  FOR  EACH  GRID. 

I NPL T  : 


GRID  SEARCH  PARAMETERS  IN  I3G)  FORMAT. 

LARGEST  AND  SMALLEST  GRID  RADII  IN  DEGREES,  AND 
GRID  SIZE  DECREASE  FACTOR  (DEFAULT  VALUE  =  0.5>. 

ANGLE  SEARCH  PARAMETERS  IN  (3G)  FORMAT. 

AS  IN  OPTION  2. 


OPTION  -4. 

THE  SAME  AS  OPTION  4.  EXCEPT  THAT  A  DECOMPOSITION  OF  THE  MEASURE 
IS  GIVEN  FOR  THE  RESULT. 


OPTION  5. 

SIMILAR  TJ  OPTION  4,  BUT  THE  PROGRAM  SEARCHES  FOR  THE  BEST 
COMBINATION  OF  POLES  AND  ANGLES  FOR  BOTH  ROTATIONS. 

THE  RADII  OF  THE  SEARCH  GRIDS  ARE  REDUCED  WHEN  BOTH  POLES  OF  THE 
BEST-MEASURE  COMBINATION  LIE  IN  THE  CENTERS  OF  THEIR  RESPECTIVE 
GRIDS. 


INPUT  : 


AS  IN  OPTION  4. 


OPTION  -5. 

THE  SAME  AS  OPTION  5,  EXCEPT  THAT  A  DECOMPOSITION  CF  THE  MEASURE 
IS  GIVEN  FOR  THE  RESULT. 


OPTION  6. 


FIXING  TH:  SECOND  POLE,  THE  PRGG  RAM  EXPLORES  A  L ATI T UDE- LONG  I T UDE 
GRID  FOR  THE  FIRST  POLE,  FINDING  THE  BEST  COMBINATION  OF  ROTATION 
ANGLES  FOR  EACH  POLE  OF  THE  GRID  AND  LISTING  THE  RESULT. 
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C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 
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c 

c 

c 
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c 
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c 

c 

c 
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c 
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c 

c 

c 
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c 

c 

c 

c 

c 
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INPUT: 

READ  THE  GRID  PARAMETERS  IN  160)  FORMAT. 

LEAST  LATITUDE,  GREATEST  LATITUDE,  LATITUDE 
INCREMENT,  LEAST  LONGITUDE,  GREATEST  LONGITUDE, 
LONGITUDE  INCREMENT. 

READ  THE  ANGLE  SEARCH  PARAMETERS  IN  13G)  FORMAT. 

AS  IN  OPTION  2. 

OPTION  -6. 

SIMILAR  TO  OPTION  6,  BUT  FOR  EACH  POLE  OF  THE  GRID,  THE  PROGRAM 
SEARCHES  FOR  THE  BEST  SECOND  ROTATION  AS  IN  OPTION  4. 

INPUT: 


LATITUDE-LONGITUDE  GRID  PARAMETERS  IN  16G)  FORMAT. 

AS  IN  OPTION  6. 

GRID  SEARCH  PARAMETERS  IN  13G)  FORMAT. 

AS  IN  OPTION  4. 

ANGLE  SEARCH  PARAMETERS  IN  1  3G )  FORMAT. 

AS  IN  OPTION  2. 

*************************************************  ***  ************** 


AFTER  PERFORMING  ITS  TASK,  THE  PROGRAM  IS  READY  TO  READ  ANOTHER 
OPTION  NUMBER  AND  CONTINUE. 

AN  OPTION  GREATER  THAN  90  WILL  TERMINATE  THE  PROGRAM. 
****************************************************************** 
MESSAGES: 


*  BEST  ANGLE  AT  END  OF  RANGE  * 

IN  THE  COARSE  SEARCH  OF  ROTATION  ANGLES,  THE  BEST  MEASURE  WAS 
ACHIEVED  WITH  AN  ANGLE  AT  THE  LIMIT  OF  THE  HALF-RANGE. 


non  noon 


190 
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c 

C  *  DOUBLY  DEGENERATE  MATRIX  POSSIBLE  * 

C 

C  THE  BEST-FIT  GREAT  CIRCLE  F CR  A  SECTION  IS  POORLY  DETERMINED, 

C  AND  MAY  6E  IN  ERROR.  BOUNDARY  AND  SECTION  ARE  GIVEN. 

C 

C 

C  *  NOT  THREE  UNEQUAL  ROOTS  * 

C 

C  THE  ROTATION  PRODUCES  A  COVARIANCE  MATRIX  FOR  A  SECTION  THAT 

C  DOES  NOT  HAVE  THREE  UNEQUAL  REAL  EIGENVALUES.  IT  IS  NOT 

C  PJSSIBLE  TO  COMPUTE  THE  BEST-FIT  GREAT  CIRCLE.  BOUNDARY  AND 

C  SECI  ION  ARE  GIVEN. 

C 

C  *  ********  *****  ******  ************  **  ****  **********  ♦  ***  ♦  ♦**  ********** 

C 

C 

C  PROGRAM  PLATE3 

C 

C 

C  ********  *  *****  **************************************  ******  ******** 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

C 

COMMCN/XMAIN/  DATA ( 3 , 5 , 10 , 2, 3  I  , S IGMA I S , 1 0 , 2, 3) , WT I  5 , 10,2 ,3) , 

L  NPTSI 10, 2,3) , NSECTI 3) 

COMMON/SUtil/  P(3,10,3) 

COMMON/SUB2/  EPI10.3) 

C0MM0N/SUB3/  DMESURl 10,3 ) ,SMESUR 13 ) ,XFACTR 
COMMON/SU  84/  RM(3,3,3I 
DIMENSION  AROT I3,3),SR0T(3) 

DATA  RAD. RKM/1.T4532 9251 994330D-02, 1.56953 90-04/ 

DATA  1 1 N,  I  OUT/  l,  108/ 

C 

WRITE! IOUT, 1  ) 

1  FORMAT !/// LX, ’PROGRAM  PL ATE3 *, 10 X ,’ JANUARY  1979’, 10X, 

Z  ’C.TAP SCOTT’,/) 

READ  NUMBER  OF  SECTIONS  FOR  THREE  PLATE  BOUNDARIES,  ANO  READ 
WFACTR . 

READ! I  IN, 2)  INSECT! I ) ,1  =  1,3) .WFACTR, IPUT 

2  FORMAT !5G) 

IF  (WFACTR  .LT.  0.1)  WFACTR=100. 

XFACTR=RKM*RKM*toFACi  R 

READ  THE  DATA  INTO  THE  DATA  ARRAY. 

CALL  GETEM!WFACTR, IPUT) 

C 


onoooo  non  non  ooo  o  o  r» 
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* 

4  - 

READ  THE  CPTION  CODE  AND  READ  THE  STARTING  ROTATIONS. 

999  READII  IN, 3  )  NOP 

3  FORMAT (1GJ 

IF  INOP  .GT.  90)  GO  TO  1000 
NC0MP=0 

IF  ( NOP  .LT.  0)  NCOMP  =  i 
WRITE! IOUT ,4)  NOP 

4  FORMAT (//IX, ‘OPTION  =  *,I3) 

NOP= I ABS I  NOP ) 

READII IN, 5)  1 1 AROTII , J) , I  =  1 , 3 > , J  =  l , 2 ) 

5  FORMAT (6GI 

IF  I IDABS IAROT 1 1,1  ))  .GT.  90.)  .OR.  I DAB  SI AROTI 1,2)1  .GT.  90.)) 
2  GO  TO  1000 

COMPUTE  THE  THIRD  ROTATION  AND  THE  MEASURE,  AND  WRITE  THEM. 

CALL  RSUMIO.AROT) 

CALL  MEASURIO, SUM) 

CALL  PUTEM 1 0, AROT , SM ESUR  »SUM ) 

GO  TO  THE  PROPER  OPTION. 

GO  TO  < IOC, 200, 200, 300, 300.400) ,  NOP 

OPTION  1.  COMPUTE  THE  MEASURE  FOR  GIVEN  ROTATIONS. 

100  IF  INCOMP  .NE.  0)  CALL  DECOMP 
GO  TO  999 

OPTIONS  2  AND  3.  FIND  THE  BEST  ANGLES  FOR  1  OR  2  OF  THE  GIVEN 
POLES. 

READ  THE  ANGULAR  SEARCH  PARAMETERS. 

200  READII IN, 20)  AL ARG, ASMAL ,ARANG 

20  FORMAT  I 3G) 

ALARG2 DABS  I ALARG) 

ASMAL2 DABS (ASMAL ) 

ARANGsDABS (ARANG) 

IF  ( ALAKG  .LT.  1.0D-05I  GO  TO  1000 
IF  (ASMAL  .GE.  ALARG)  ASMAL=0. 

IF  (ARANG  .LT.  ALARG)  ARANG2 ALARG 
C 

GO  TO  1201,201,201,301,301,401),  NOP 
C 

201  WRITE! IOUT, 21)  ALARG, ASMAL, ARANG 

21  FORMAT  I //IX, ‘ANGLE  SEARCH:  LARGER  STEP  *,F8.3/15X, 

Z  ‘SMALLER  STEP  • ,F 8.3/ 15X, ‘RANGE ' , 10X, F8.3 ) 


OOOO  O  O  O  O  o  o 
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i 


CALL  SANGL  (  (3-NOP)  ,ARCT  , ALARG, AS MAL, ARANG, SUM,  1) 

C 

K=NOP— 2 

WR I T  E  4  I  OUT  *22)  K 

22  FORMAT  I// IX, 'BEST  ANGLES  FCR  BOUNDARIES  *,11,'  2  3') 

CALL  P UT  t;  M  { 0  *  A  RO T »  SM  E S UR  ,  S  UM  ) 

IF  INCOMP  .NE.  0)  CALL  DECOMP 
GO  TO  999 
C 

C  OPTIONS  4  AND  5.  SEARCH  ONE  OR  BOTH  PCLES. 

C 

300  READ! I  IN#  20)  XLARG , X SMAL ,XMULT 
XLARG=DABS  tXLARG) 

XSMAL*DABS (XSMAL ) 

XMULT=DAti  S IXMULT  ) 

IF  ( (XLARG  .GT.  90.)  .OR.  (XLARG  .LT.  1.0D-05))  GO  TO  1000 
IF  UXSMAL  .GT.  XLARG)  .OR.  (XSMAL  .LT.  1.0D-05))  XSMAL=XL ARG 
IF  ( XMULT  .GE.  1.)  X  MULT  =0 . 

GO  TO  209 

301  WRITE! lOJT  ,30)  XLARG , ALARG, X $MAL  ,ASMAL , XMULT , ARANG 
30  FORMAT!/ /IX, 'POLE  SEARCH: • ,/ 3X, • LARGEST  STEPS  • ,2! 2X,F8.3) , /3X, 

2  'SMALLEST  STEPS  •  ,2  ( 2X  ,  F8 .3  )  , /3X,  •  F  ACTOR  AND  RANGE  '  ,  2!  F8 . 3 , 2X  )  ) 

CALL  SPOL  I  (5 -NOP)  ,  AROT ,  XLARG , XS.MAL ,X MULT  , A LARG ,  ASMAL ,  ARANG,  SUM ,  I ) 

IF  ( NCOMP  .NE.  0)  CALL  OECOMP 
GO  TO  999 

OPTION  6.  EXPLORE  A  GRID  FOR  THE  FIRST  POLE.  EITHER  FIX  THE 
LOCATION  OF  THE  SECOND  POLE,  OR  SEARCH  FOR  THE  BEST  FIT. 

READ  THE  GRIO  PARAMETERS. 

400  READ! I  IN, 5  )  BOT , TOP, VST EP , XL  EFT, R IGHT ,HS TE P 
WRITE! IOJT,40)  BOT , T OP, VST EP ,XLE FT ,R IGHT ,HSTEP 

40  FORMAT  (//IX, 'GRID*.  •  ,20X,' LIMITS'  ,7X,'GRID  S I  2E  •  , /6X ,  '  LAT  I TUDE*  , 

2  4X , 3! 2X , F  8. 3) , /6X , ' LONG  I TUDE ' ,3 X, 31 2X ,F 8. 3 ) ) 

READ  THE  POLE  SEARCH  PARAMETERS  IF  NECESSARY,  AND  READ  THE  ANGULAR 
SEARCH  PARAMETERS. 

IF  (NCOMP  .EQ.  0)  GO  TO  200 
GO  TO  300 

401  IF  INCOMP  .EQ.  0)  WR  I  TE (ICUT , 2 1 )  ALARG, ASMAL, ARANG;  GO  TO  402 
WRITE! IOUT  ,30)  XLARG , ALARG ,XSMAL  , ASMAL , XMULT , ARANG 


4C2  NV*(TOP-BOT) /VSTEP 


OOO  OOOOOO 


NH=(RIGHT-XLEFT)/HSTEP 

C 

DO  600  I V  =  0»  NV 

AROT (1.1) =80T*VSTEP*FL0AT( IV) 

C 

00  500  IH=0»NH 

A ROT (2*1 )  =  XLEFT*HSTEP*FL0AT( IH) 

IF  (NC0MP  .NE.  0)  GO  TO  450 
C 

CALL  SANGL(0,AROT,ALARG,ASMAL,ARANG,SUM,0) 

CALL  PUTEM(IH,AROT,SMESUR,SUM) 

GO  TO  490 
C 

450  CALL  SPOL < 1 , AROT .XLARG, XSMAL , XMULT , ALARG ,A$MAL , ARANG .SUM , 0 ) 

C 

490  IF  (IH  .ME.  0)  GO  TO  500 
SANG=AROT (3,1) 

DO  491  1=1,3 

491  SROT (I )=AROT (1,2) 

C 

500  CONTINUE 
C 

AROT (3, 1) =SANG 
00  501  1=1,3 

501  AROTII ,2)=SR0T(II 
C 

600  CONTINUE 
GO  TO  999 
C 

1000  WRITE! I OUT, 1001) 

1001  FORMAT!//) 

STOP 
END 

******  ****** ********  ************************************ ********** 
SUBROUTINE  OECOMP 

SUBROUTINE  OECOMP  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  COMPUTES  AND  WRITES  THE  QUANTITIES  INVOLVED  IN 
THE  DECOMPOSITION  OF  THE  MEASURE. 

IMPLICIT  OOUBLE  PRECISION  (A-H.O-Z) 

C 

COMMCN/XMAIN/  DATA(3«5*10»2«3 ) ,S IGMA (5 » 10, 2, 3) ,WT(5, 10,2,3), 

Z  NPTSI 10,2,3), NSECTI 3 ) 

COMMON /SUB1/  PI  3,10, 3) 

C0MM0N/SU83/  DMESURI 10,3 ) , SMESUR (3 ) , XF ACTR 
COMMON/ SUB 4/  RM(3,3,3) 
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0  I  MENS  ION  DEV!  5)  ,T { 3) 

DATA  RAO. RKM.DKM/l. 74532 92519943300-02,1.5695390-04,  1.  1120*02/ 
DATA  1  OUT / 108/ 

C 

MR  1 TE ( I OUT *  1 ) 

1  FORMAT {//IX, ‘DECOMPOSITION  OF  ME ASURE* ,//5X, 

L ‘MEASURE  =  SUM  OF  S9  .  OEV.  OF  POINTS  IREL.  TO  ST.  DEV.  OF  ERROR)') 
WRITE! 10UT .2) 

2  FORMAT (/24X»* EIGENVECTOR* , /IX, 

l  ‘BOUNDARY  SECTION  LATITUDE  LONGITUDE  MEASURE') 

C 

DO  100  13=1,3 
WRI TE( IOJT.iO) 

10  FORMAT (•  •) 

DO  100  IS=1,NSECT( 13 ) 

C 

PL AT  =OAS I N (P ( 3 , IS, IB ) J/RAD 
PLON=DATAN2(P(2. IS . I B) ,P t 1 , I S , IB))/RAD 
C 

D=DMESURI I S, IB)/XFACTR 
WRITEUGUT.il)  IB,IS»PLAT,PLON,D 

11  FORMAT  C4X,I1,7X,I2,6X,2(F8.3,2X) ,D12.6) 

C 

100  CONTINJE 
C 

WRITE! IOUT ,3) 

3  FORMAT  I/5X. ‘DEV.  OF  POINTS  FROM  GREAT  CIRCLE  IREL.  TO  ST.  DEV.  OF 
ZERROR) •) 

WRITE l IOUT ,4  )  11,1=1,5) 

4  FORMAT I/IX.'BNDRY  SECT  PLATE ' , 5! 5X, ‘ PT . * , l 1, IX ) ) 

C 

DO  400  113  =  1,3 
WRITE! IOUT. 10) 

DO  400  1 S  =  1 ,NS  ECT ( 13 ) 

C 

NA=NPTS! I S , 1 , I B) 

DO  200  J=1,NA 

DEV(J)  =  DATA(l,J,IS,l,I3)*P{l,IS,IB)«-DATA(2,JtIS»l,lB)*P(2,IS,IB)-'- 
Z  DATA!3,J»IS»1»IB)*P!3»IS»IB) 

DEV! J)=DACOS !DEV { J )) /RAD-90. 

200  DEV!  J)  =DE  V I  J )*DKM*RKM/SI GMAl J,  IS, l, IB) 

C 

K=1 

IF  I  IB  .ED.  3)  K =2 

WRITE! 10UT.21)  I B , IS , K , { DEV l J ) , J= l , NA ) 

21  FORMAT {3X»I1,4X,I2«4X,II,5X,5(F8.3,2X) ) 

C 

NA=NPTS! IS, 2, IB) 

DO  300  J=1 ,NA 

CALL  ROT AT E!RM I 1,1, IB), DAT A! 1 , J, IS , 2 , I B) , T ) 


onnoooo  non 
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DEVI  JJ  =  H  1)*P{  1»IS»IB)*T(2)*P(2»IS»IB)*-T(3)*P(3,IS»IB) 
DEV  I J) -DA COS (D£V(J))/RAD-90. 

300  DEV(J)=DEV(J) *DKM*RKM/S I GMAI J, IS *2. IB) 


K=3 

IF  (IB  .EO.  1)  K =2 

WRITE(I0UT«21)  IB, I S ,K , ( DEVI J » , J =1 ,NA ) 

C 

400  CONTINUE 
C 

RETURN 
END 

****************************************************************** 


SUBROUTINE  SPOL(N»ARQT»XL»XS»XM,  AL , AS , AR , FSUM, I SUM) 

SUBROUTINE  SPOL  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  SEARCHES  FOR  THE  BEST  PAIR  OF  POLES  AND  ANGLES  IN 
HEXAGONAL  GRIDS  DECREASING  IN  SIZE  BY  FACT CR  XM  FROM  XL  TO  XS. 

IF  N  .NE.  0  ONLY  THE  SECOND  POLE  IS  VARIED. 

IMPLICIT  DOUBLE  PRECISION  (A-H,0>Z) 

C 

COMMON/ SUB 3/  DM£SUR( 10, 3  I *SMESUR C3 ) • XFAC TR 
JlMEAS I GN  AROT (3 , 3  I , R I  3, 3 1 »PLAT l 7, 2 ) , PLON ( 7, 2 ) 

DIMENSION  IAC0DEI7I,  IBC0DEI7) 

DATA  IUUT/108/ 

C 

M=N 

IF  (M  .NE.  0)  M- 1 
M=M*l 

IF  1 1 SUM  .NE.  01  GO  TO  1 
CALL  RSUM (0, AROT ) 

CALL  MEASUR(0*FSUM) 

C 

1  SIZE*XL 
NSTEP*0 
IA»1 
I  TA*1 
I  TB*1 

DO  LOO  1  =  1,3 
00  100  J* 1  ,2 
100  R ( I , J)  =  AROT ( I , J  ) 

WRITE! I OUT, 10) 

10  FORMAT  </ IX, 'BOUNDARY  LATITUDE  LONGITUDE  ANGLE *,6X, 
l  ‘MEASURE  TOTAL  MEASURE*) 

00  101  1*1,7 
I ACOOE 1 1 ) *0 


196 


101  I BC00L (11=0 
C 

GO  TO  1200,201 1,  M 
C 

?C0  CALL  PSfc  T (SIZE, ARGT (1,1) » AROT ( 2, 1), PL AT (1*1) , PLUN  (  l , 11) 

201  CALL  PSEKSIZE.AROTI 1,2) ,ARGT ( 2, 2) ,PLAT( 1,2) ,PL0N(1,2I) 

GO  TO  (300,301),  M 
C 

3C0  00  401  l  A*  1, 7 

R(1,1)=PLAT(1A,1) 

R(2,1)=PL0N(IA, 1) 

C 

301  00  400  10*1,7 

IF  ((IACOOE(IA)  .NE.  0)  .AND.  UBCODE(IB)  .NE.  0))  GO  TO  400 
K(  1,2)  =PL  AT  ( 10,2  ) 

R(2,2)=PLCN(IB,2) 

C 

CALL  SANGL (0,R,AL,AS ,AR,SUM,0) 

C 

IF  (FSUM  .LE.  SUM)  GO  TO  400 

I  TA*I  A 

1 T  B*  1  d 

FSUM=SUM 

00  302  1*1,3 

00  302  J= 1 »  2 

302  AROT l I , J) =R( I , J) 

C 

4C0  CONTINUE 

GO  TO  (401,402),  M 

401  CONTINUE 
C 

402  NSTEP=NST EP*1 
CALL  RSUM (0, AROT ) 

CALL  MEASUR I  0, FSUM) 

C 

I TC*ITA*I TB 

IF  ( IT C  .NE.  1)  GO  TO  500 
SIZE=SIZE*XM 

IF  (SIZE  .GE.  XS)  GO  TO  500 
WRITE! IOUT, 50) 

50  F0RMATI/1X,' SEARCH  STOPPED') 

500  WRITE! IOUT, 51)  NSTEP 

51  FORMAT (/IX, 'STEP* ,15) 

CALL  PUTEMd, AROT, SMESUR, FSUM) 

IF  (SIZE  .LT.  XS)  GO  TO  1000 
C 

GO  TO  (600,601),  M 
C 

600  IF  ((ITC  .EQ.  1)  .OR.  (ITA  .NE.  11)  CALL  PSET (SI ZE»AROT( l, l ) , 
Z  AROT (2,1 ),PLAT(1,1)  ,PL0N(1, 1) ) 


ooooooo  non 


6C1  CALL  ISETCITC»ITA*  I A  CODE  ) 

C 

IF  <  ( I  TC  .EQ.  1)  .OR.  (ITB  .NE.  1)>  CALL  PSE  T  C  SI  ZE,  AROT(  l,  2 1  , 
l  AROT <2»2)tPLAT(l»2)» PLON 11*2)1 
CALL  IStU  ITCt  ITB.IBCODE) 

C 

I  TA=1 
I  TB=  1 

GO  TO  ( 300*301 ) «  M 
C 

1000  RETURN 
ENO 

****************************************************************** 

SUBROUTINE  SANGL IN. AROT , AL#A S , AR »FSUM, I  SUM) 

SUBROUTINE  SANGL  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  SEARCHES  FOR  THE  BEST  COMBINATION  AF  ANGLES  FOR  THE 
GIVEN  POLES.  IF  N  .NE.  0  THEN  ONLY  THE  ANGLE  OF  THE  SECOND  POLE 
IS  VARIED. 

IMPLICIT  DOUBLE  PRECISION  CA-H,0~Z) 

C 

DIMENSION  AR0TC3,3),RC3,3},CC2) 

DATA  I0UT/108/ 

C 

M=N 

IF  CM  .NE.  01  M= I 
IF  (ISOM  .NE.  0)  GO  TO  1 
CALL  RSUMCO ♦ AROT ) 

CALL  MEASURt  0»  FSUM ) 

C 

1  ITA=0 
I  TB=0 
I  A=0 
DIF=AL 
NL»AR/AL 

IF  (NL  .EQ.  0)  NL= 1 

IF  CAS  .LT.  1.0D-05I  NS*0;  GO  TO  100 
NS=1.5*AL/A$ 

IF  CAS  .GE.  AL)  NS=0 
C 

100  LTI ME*1 

DO  101  Jsl»  2 
DO  101  I  =1  »3 

101  RCI » J)=AROTI I» J) 

C 


999  CCliaAR0TC3,l) 


o  r>  n  o  o  n  n 
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T 


C (2  ) =AROT (  3,2) 

IF  IM  .NF  .  0)  GO  TO  102 
C 

DO  201  I  A  =  -NL , NL 
R  ( 3 , 1)  =C  (  1 )  *-FLOAT  <  IA ) I F 
C 

102  00  200  IO--NUNL 
K(3»2)=C(2)*-FL0ATt  I B ) *D I F 

C 

MM=  l 

IF  l(M  .FO.  0)  .AND.  I  IB  .EQ.  -NL))  MM=0 
CALL  RSUM(MM,R) 

CALL  ME AS UR I  MM , SUM ) 

IF  (FSUM  .LE.  SUM)  GO  TO  200 
C 

AROT  (3,2)=  R  (3, 2) 

IF  (M  .EJ.  0)  AROT 13 ,1)=R(3, 1) 

DO  103  1=1,3 

103  AROT ( I ,3)=R( 1,3) 

FSUM=SUM 
ITA-IABSl IA) 

1TB=1ABS( IB) 

C 

200  CONTINUE 
IF  (M  .N£.  0)  GO  TO  202 

201  CONTINUE 
C 

202  GO  TO  1300,400)  LTIME 
C 

300  IF  KITA  .FQ.  NL)  .OR.  (  ITB  .EQ.  NL )  )  WR 1  TE  (  IOUT ,  30) 

30  FORMAT I/IX,'*  BEST  ANGLE  AT  END  CF  RANGE  *') 

IF  (NS  .EQ.  0)  GO  TO  400 
NL=NS 
OIF=AS 
LT I ME=2 
GO  TO  99 ) 

C 

4C0  CALL  RSUM  CM , AROT ) 

CALL  ME AS UR ( M, FSUM ) 

RETURN 
END 

****************************************************************** 
SUBROUTINE  MEASURIN, SUM) 

SUBROUTINE  MEASUR  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  CALCULATES  THE  MEASURE  CF  FIT  FOR  EACH  BOUNDARY, 


L 


■  ■ 


noon 
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i 

EXCEPTING  BOUNDARY  N,  FOR  THE  ROTATION  MATRICES  RM.  THESE  ARE 
SUMMED  TO  OBTAIN  THE  OVERALL  MEASURE,  ACCORDING  TO  S.HELLINGER* S 
CRITERION  PHI-O. 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

COMMGN/XMAIN/  OATAI3 ,5, 10,2, 3 I,S IGMA ( 5. 10. 2,  31  ,  WT ( 5, 10*2,3) , 

Z  NPTSI 10, 2*3), NSEC TI3) 

C0MM0N/SUB1/  P  (3 ,10, 3  ) 

COMMON/ SU £2/  EP(10,3) 

COMMON /SUB 3/  DMESURI 1 0, 3 ) , SME SUR 13 ) ,XFACTR 
C0MM0N/SUB4/  RM(3,3,3) 

DIMENSION  DMERGE (3,10),WMERGE(10) 

DATA  I  OUT/ 108/ 

DO  300  18=1,3 
IF  IN  .Ed.  IB)  GO  TO  300 

SMESUR ( 131=0. 

DO  200  K=1,NSECT (IB) 

NA=NPT$ IK, 1,18) 

DO  100  1=1, NA 
H MERGE (I)=WTtl,K«l*IB) 

DO  100  J= 1 , 3 

100  DMERGE ( J , I )=DATA(J,I ,K,1,IB) 

NB=NPT  S (K , 2, IB) 

DO  101  1=1 »NB 
W  MERGE  (  I  *NA )=WT  I  I,K,2,IB) 

101  CALL  R0TATE(RM(l,l,lB),DATA(l,I,K,2,IB),DMERGE(L*l+NA)) 

NAB=NA+N3 

CALL  E I GEN( DMERGE, WMERGE ,NAB ,P(1,K,IB),EP(K, IB),  NFL AG) 

GO  TO  (104,103,1021,  NFL  AG 
C 

102  WRITE! I OUT ,10)  IB,K 

10  FORMAT (IX,**  NOT  THREE  UNEQUAL  ROOTS  **, 12X, * BOUNDARY,  SECTION  =  • 
Z  ,11, IX, 12) 

GO  TO  104 

103  WRITE! IOJT, 11)  IB.K 

11  FORMAT ( IX, '*  DOUBLY  DEGENERATE  MATRIX  POSSIBLE  *  BOUNDARY,  SECTIO 
ZN  =  *, 11, IX, 12) 

C 

104  DMESURI K, !B)=EP(K, IB) 

C 

200  SME SUR ( IB) = SME SUR ( IB ) ♦DMESURI K, I B) 

C 

300  CONTINUE 
C 


on  ooooooooooooooooooooono  non 


5 
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SUdROUTINfc  EIGENIX.WT.N.P, EP,NFLAG> 


SUBROUTINE  EIGEN 


NOVEMBER  1978 


S.HELL I NGER 


P,0, R  ARE  EIGENVECTORS  ASSOCIATED  WITH  EIGENVALUES  EP » EQ  » ER 
X  =  L  INEAR  ARRAY  CONTAINING  DATA  POINTS  IN  CARTESIAN  COORDINATES 
WT  =  LINEAR  ARRAY  WHICH  CONTAINS  WEIGHTS  FOR  DATA  POINTS.  WTIK) 
PROPORTIONAL  TO  SQUARE  OF  STANDARD  DEVIATION  OF  ERROR  FUR  K  TH 
POINT.  1 WT  =  INTEGER  VARA  I ABLE  TO  INDICATE  WHETHER  DATA  IS  TO  BE 
WEIGHTED  ACCORDING  TO  WT  ARRAY.  IWT=0  INDICATES  WEIGHTING  10  BE 
USED.  IF  IWT.NE.O  THEN  NO  WEIGHTING  IN  FINDING  EIGENVECTORS 
AND  EIGENVALUES. 

N=NUMBER  OF  DATA  POINTS 

V EC  =  3X 3  ARRAY  CONTAINING  EIGENVECTORS,  STORED  COLUMN-WISE 
0-L I NE  AR  ARRAY  CONTAINING  EIGENVALUES.  ELEMENT  L  CONTAINS 
EIGENVALUE  FOR  EIGENVECTOR  IN  COLUMN  OF  ARRAY  VEC. 

A=3X3  COVARIANCE  MATRIX.  SIX  INDEPENDENT  ELEMENTS  STORED  IN 
LINEAR  ARRAY  S16). 

S II >=A(1,  1  ) 

S(2)sA(l*2)-A(2tL) 

S  I3)=A(2,2I 
S(4)=A(  l,3)  =  A(3,  1) 

S<5) -A (2, 3)=A(3, 21 
SI6»=A13,3» 

IMPLICIT  REAL*8  (A-E.P-Z) 

DIMENSION  VEC13.3) ,A(3«3),P(3)VQ(3),R(3) ,S(6),D(3),CD(3) 
DIMENSION  WTIlOI.Xm 

DATA  C/0.  74074074074 074  1D-0 1/, Cl 20/ 0.  20943951023931 9Dm/ 

DATA  C24J/0. 41 88 7902 04786390*01/ .XZERO/1 .00-05/ 

NFLAG= 1 

CREATE  SYMMETRIC  COVARIANCE  MATRIX  FROM  DATA. 


IS 


**.N0  WEIGHTING 

FOR  N 

=  2** 

DO  5  I=l» 

6 

S (11=0. 

IFIN.E0.2 

)  GO 

TO 

8 

WEIGHTED 

SOLUT 

ION 

DO  6  M=1 * 

N 

MM=3*M 

sm=sm 

♦X  { MM 

-2 

)*X( 

MM- 

2 

J/WTIM) 

S I  2 )  =S  1 2) 

♦•XI  MM 

-2 

)*X( 

MM- 

1 

I/WTIM) 

SI3I=S  C3J 

♦  X(MM 

-I 

)*X< 

MM- 

1 

)/ WT (M) 

SI4MSC4) 

♦  XI  MM 

-2 

l*X( 

MM) 

/ 

WTIM) 

S  C  5  >  =S  C  51 

♦  XI  MM 

-1 

)*XI 

MM) 

/ 

WTIMJ 
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6  S(6)=S(6) +X(MM)*X(MM)/WT ( M) 

GO  TO  9 

C  UNWEIGHTED  SOLUTION 

8  DO  LO  M=i,N 
MM=3*M 

SU)=S(1)  *X(  MM-2)  *X(  MM-2) 

S (2 )=S ( 2)  *-X  ( MM-2  )  *X(  MM-1 ) 

S13)=S(3) *X( MM-1 )*X( MM-1 ) 

S ( 4) =S l 4)  *-X(  MM-2  )  *X(  MM) 

SIS)  — S  (31  *-X  (MM-1  )  *X(  MM  ) 

10  S(6)=S(6)  *-X(MM)*X(  MM) 

9  CONTINUE 

C  SOLVE  CUBIC  OE TERM  IN ANTAL  EQUATION  FOR  EIGENVALUES 

PP=-(S(1)  *S(3)fS(6>) 

GJ=S(1  )*SI6USI1)*SI  3)*S(3)*S(6)-SI5)*SC5)-S(2)*SI2)-S(4)*S(4) 
RR=S ( l ) *S I  5) *S15)+S( 2)*S(2)*S(6) fS(4)*SI 4)  *S  ( 3)- S 1 1 )  *S  ( 3  ) * S(  6) 

2  -2.*S(2)*S(4)*S(5 ) 

C  SOLVE  REOUCEO  EQUATION 

AA=Q0-PP*PP/3. 

BB=C*PP*PP*PP-PP*0G/3.  *■  RR 
A A=AA/ 3  « 

BB=8B/2. 

C  CHECK  FOR  THREE  REAL  UNEQUAL  ROOTS 

aaa«aa*aa*aa 

SP=BB*  BB* AAA 
IFISP.GE.O.)  GO  TO  50 

C  HAVE  THREE  REAL  UNEQUAL  ROOTS.  FIND  THEM 

PHI -DARCQS I  -BB/DSQRT I-AAA)  ) 

PHIsPHI/3. 

AAA=2.*0SGRT(-AA) 

0(1 )=AAA*  DCOS 1  PH  I ) 

D(2)=AAA*DC0Sl PHI +C 120) 

0(3) =AAA$DCOS ( PHI ♦C240J 

C  D(1)»D(2)»D(3)  ARE  ROOTS  OF  REDUCED  EQUATION.  FIND  ROOTS 

C  OF  ORIGINAL  EQUATION. 

PP=PP/3. 

0(1  l=D( 1 1-PP 
D(2)=D(2)-PP 
D 1 3 )-D ( 3 ) -PP 

C  FIND  MINIMUM  EIGENVALUE  AND  PLACE  IT  IN  DU) 

IF(  (D(  1)  .LE.D  (2  ) )  .AND.  (0(2.)  . LE.O(  3)  )  )  GO  TO  30 
MI  N=3 

I F(  0(2) »LE.D( 3)  )  M I N=2 
C  EXCHANGE  EIGENVALUES 

C0(  1)  =  D(  1  ) 

D(i)=D(MIN) 

D ( MI N) =CD (I I 

C  COMPUTE  EIGENVECTORS 

C  HAVE  ROOTS  OF  ORIGINAL  EQUAT I CN=E IGENVALU ES  OF  COVARIANCE  MATRIX 

C  FIND  EIGENVECTORS  ASSOCIATED  WIOH  EIGENVALUES.  DONE  BY  FINDING 
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C  LINEARLY  INDEPENDENT  ROWS  OF  MATRIX  A-D(J)*I,  J=l,3  ,  1= IDENTITY 

C  MATRIX.  EIGENVECTOR  IS  THEN  IN  DIRECTION  OF  CROS S-PRUDUC 1  OF 

C  LINEARLY  INDEPENDENT  ROV-  VEC  TORS  . 

30  A ( 1 » 2 ) =S ( 2  ) 

A I  2  ,  L )  -  S  l  2  ) 

A ( I »  3) =S ( 4 ) 

A! 3,i)=S<4) 

A  1  2, 3)  =S(  b  > 

A  l  3  » 2 ) =S I b  ) 

J  =  1 

A ( 1  *  1 ) =$ (  1J-DI J) 

A  ( 2 , 2 ) =S ! 3 I-0( J ) 

A  I 3»  3) =S I 6)-Dt  J) 

C  TEST  FOR  LINEAR  INDEPENDENCE  OF  ROW-VECTORS  (=COLUMN-VEC 10RS ) 

DC  440  M - 2 ♦ 3 
DO  44  L=  L  »  3 

44  R!L)=A(1,L)/A(M,L) 

MROW=M 

IF!  DABS  I  R I  1 >  — R 121  l.GT.XZERO  >  GO  TO  441 
If (  DABS  I  Rill — R 131  l.GT.XZERO  )  GO  TO  441 

440  CONTINUE 

C  MAY  HAVE  DOUBLY  DEGENERATE  MATRIX.  PRINT  WARNING  MESSAGE.  COMPUTE 

C  VECTOR  PRODUCT  OF  PAIRS  CF  ROWS.  CHOOSE  VECTOR  WITH  GREATEST 

C  MAGNITUDE  TO  BE  DIRECTION  OF  EIGENVECTOR.  NORMALIZE  LENGTH  TO  1. 

NFLAG=2 

CALL  CROSS!  A! 1, 1)  ,A < 1,2)  ,P  ) 

^ ALL  CROSS!  A! L.l)  ,A { 1,3  1,0  ) 

CALL  CROSS!  A ! L , 2 )  ,  A ( 1 , 3  ) , R  ) 

C0ll)=PU)*P!  1>*P( 2)*P(2)*PI3)*P(3) 
CD(2)=U(l)*Q(l)*Q(2)*C(2)+Q(3)*Q(3) 

C0(3)=R(1)*K(1)»R(2) *R(2)*R( 3)*R(3) 

MM=  I 

IF!  CD(2)  .GT.CDI  I)  )  MM=2 
IF!  CD(MM).LT.CD(3>  )  MM=3 
SP=DS0RT I  CD(MM)  ) 

GO  TO  (4402,4403,4404),  MM 

4402  DO  44021  1=1,3 
44021  VEC( I»J)=P(I )/SP 

GO  TO  40 

4403  DO  44031  1=1,3 
44031  V  EC  (  I,  J )  =  w(  I  )/SP 

GO  TO  40 

44C4  DO  44041  1=1,3 
44041  V  EC  (  I,  J)=R(I  )/SP 
GO  TO  40 

C  RCW  (MROW)  IS  LINEARLY  INDEPENDENT  OF  ROW  1.  COMPUTE  EIGENVECTOR. 

441  CALL  CROSSt  A  (  1*1)  ,  A  (  1 ,  M  RO  W )  ,  V  EC  (  1 »  J  )  ) 

C  VEC(ltJ)  NCW  CONTAINS  FIRST  ELEMENT  OF  UN-NORMAL  I ZEO  EIGENVECTOR. 

C  MAKE  LENGTH  =1. 

AAA  =  VEC1 l » J ) vVEC (1 , J ) *VEC (2, J) *VEC (2, J )*VEC( 3, J ) *VEC I  3 » J ) 


o  o  o  o  o  o  ooo 
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A 

V 


AAA=DS JRT ( AAA) 

DO  4413  1*1,3 

4415  VECII,J)=VEC(I , J ) / AA A 
C 

40  00  41  L=l,3 

41  P(L)=VEC(L,1) 

EP=DI 1 ) 

RETURN 

C  TEST  FOR  THREE  UNEQUAL  ROOTS  FAILED.  PRINT  MESSAGE. 

50  NFL AG=  3 

RETURN 
•  END 

*********  ********************************************************* 
SUBROUTINE  PSETI S, XL AT ,X LON, PL AT ,PLONI 

SUBROUTINE  PSET  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  RETURNS  A  HEXAGONAL  GRID  OF  LATITUDE-LONGITUDE 
COORDINATES  OF  RADIUS  S  CENTERED  AT  XLAT ,  XLON. 

IMPLICIT  DOUBLE  PRECISION  CA-H.O-Z) 

C 

DIMENSION  PLAT  (7  I , PLONI7 I ,R< 3, 3)  ,PI3>,Q<  31 
DATA  RAD/ 1. 7 45 32 92 5 l 99433 0D- 02/ 

C 

PLAT (1 ) =XLAT 
PLONI1 ) =XLON 
C 

P(l)=0. 

P(2I=XL0N*90. 

P(3)=90.-XLAT 
CALL  RMAT (P,R! 

C 

T=0. 

DT=60.*RAD 
P  1 3 )=DCOS I S*RAD) 

DP*DSINIS*RAD) 

C 

00  100  1*2,7 
PC1)=DP*DC0SIT I 
P12I*DP*DSIN(T ) 

CALL  ROT ATEIRtPtQ) 

PLAT (I)=90.-(DACOSCQI3) J/RAD) 

PL0N1I ) *0ATAN2 ( Q12 )* Q ( 1) I/RAD 
ICO  T*TFDT 


(  C 


RETURN 

END 


204 


C 

C  * ■?/**** ***  ******##***++$******$*******  *****#**********  +  * 

C 

SUBROUTINE  RSUM { N»  AROT ) 

C 

C  SUBROUTINE  RSUM  JANUARY  1979  C.TAPSCOTT 

C 

C  THE  SUBROUTINE  RETURNS  THE  RESULTANT  POLE  AND  ANGLE  AND  ROTATION 

C  MATRICES  FCR  ROTATIONS  L  AND  2,  LEAVING  THE  ROTATION  MATRIX  FOR 

C  ROTATION  N  =  0 » 1 »2  UNCHANGED. 

C 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-ZI 
C 

COMMON /SUB4/  RM(3,3,3) 

DIMENSION  AKOT (3,3), 3(3, 3) 

C 

IF  IN  .EJ.  1)  GO  TO  100 
CALL  RMAT (AROT (1 » 1 ) »  RM( 1,1,11) 

C 

100  IF  (N  .EU.  2)  GO  TO  200 

CALL  RMAT ( AROT (1,2), RM (1,1, 2)) 

C 

200  DO  201  1-1,3 
DO  201  J*  1 , 3 
2C1  0(I,J)*RM(J,I,1) 

C 

CALI.  MM33  1 0,  RM  (  1, 1 , 2  )  ,RM  1 1 , 1 , 3 II 
CALL  R I NVERT (KM (1*1,3), AROT (1,3) ) 

C 

RETURN 

END 

C 

C  ****************************************************  ********  ****** 

C 

SUBROUTINE  ISET(L,M,NJ 
C 

C  SUBROUTINE  ISET  JANUARY  1979  C.TAPSCOTT 

C 

C  THE  SUBROUTINE  SETS  THE  ELEMENTS  OF  THE  ARRAY  N  AS  CODES  TO  AVOID 
C  UNECESSARY  RECALCULATION  IN  SUBROUTINE  SPOL. 

C 

DIMENSION  N l 71 
C 

N  ( 1 1  — 1 


DO 

100 

I  -2  i 

,7 

Nil 

1 1*0 

IF 

IL  , 

.ED. 

1) 

RETURN 

IF 

IM  , 

,GT. 

1) 

GO  TO  1C2 

DO 

101 

1*2, 

*7 

oooono  o  o  o 
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101  N  1 1 J  =1 
RETURN 

102  DO  103  1  =  M«-2,M*4 
J=I 

IF  CJ  .GE.  8)  J=J-6 
1C3  NIJ)=l 

RETURN 

END 

************ ****** ***^** ********************************  ********** 


SUBROUTINE  GET EM ( W , I  PUT ) 

SUBROUTINE  GETEM  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  REAOS  THE  DATA  SECTIONS*  CONVERTS  THE  POINTS  TO 
CARTESIAN  COORDINATES,  AND  PLACES  THEM  IN  THE  DATA  ARRAY. 


( 


IMPLICIT  DOUBLE  PRECISION  <A-H,0-Z> 

C 

COMMCN/XM AIN/ DAT  A(3*5*L0t2»3),SI GMA (5* 10 *2*3), WT (5*10*2*3) , 
Z  NPTSI 10, 2,3)* NS EC T{ 31 
DIMENSION  ALATt5),AL0N(5),AERR(5) 

DATA  RAD, RKM/1. 7453292519943300-02, 1.569539D-04/ 

DATA  I  IN  *  I  OUT/ 1*108/ 

C 

WRITE! I OUT  » 1 1 

1  FORMAT!/, 'BOUNDARY  PLATE  SECTION  POINTS',/) 

C 

DO  200  I B-  1 , 3 
DO  200  IP-1.2 
K=IP 

IF  ( 1 1 B  .EO.  2)  .AND.  (IP  .EQ.  2))  K=3 
IF  (IB  .EC.  31  K=IPH 
00  200  I S=1 * NSECT( IB ) 

C 

READ! I  IN. 10)  ! ALAT ( J ), ALON! J),AERRIJ),J= 1,5) 

10  FORMAT (15G) 

IF  ( IPUT  .NE.  0)  GO  TO  99 

WRITE! 10UT  *11)  IB.K, 1 S, ! ALAT I J ) * ALON ( J) , AERR ! J ) ,J=1,5) 

11  FORMAT (1X*I1*1X«I1*1X*I2*3X*5(F6.2*1X*F7.2,1X*F4.1*5X)) 

C 

99  DO  100  J-1,5 

IF  ( (ALAT ( J)  .GT.  90.)  .OR.  (AERR! J)  .LT.  0.01))  GO  TO  101 
C 

CLAT*DCQS ! RAD* AL AT ( J ) ) 

SLAT«OSI N( RAO* ALAT (J ) ) 

CLGN-DCOS ( RAD* ALON IJ ) ) 


ooooo  o  o  o  o  r  o  o  o  noooo  ono 
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SLCN=DSl N<  RAD*AL0N  !J  )  ) 

C 

DA  T  A ( 1 f  J ,  IS, IP,I B)=CLAT*CLON 
DAT  A 12, J,  IS, IP, IB) =C  LAT*  SLON 
DAT  A I  3 , J , IS, IP, IB) =  $LAT 
SIGMAl J , IS  ,IP, IB)  =  AERR( J)*RKM 
100  WT(J,lS,IP,l B)=AERRl J)*AERR< J)/W 
C 

1C1  NPTSI IS, 1 P, IB) =J-1 
C 

200  CONTINUE 
C 

RETURN 
END 

****  *****  *******************  ********************:****  ************** 
SUBROUTINE  PUTEM II T, A ,S, T) 

SUBROUTINE  PUTEM  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  WRITES  THE  ROTATION  POLES  AND  ASSOCIATED  MEASURES. 
DOUBLE  PRECISION  A (3 ,3) , SI  3) ,T 
DATA  I0UT/108/ 

IF  I  IT  .NE.  0)  GO  TO  100 
WRITE! I OUT, l  ) 

l  FORMAT (//IX, 'BOUNDARY  LATITUDE  LONGITUDE  ANGLE*, 6X, 

Z  *  MEASURE  TOTAL  MEASURE*) 

ICO  WRITE! IOUT  ,10)  ( J , ( A  II , J ) , 1  =  1 , 3 )  ,  S  l  J ) , J= 1 • 3) ,T 
10  FORMAT I3l/4X,I1,5X,3IF8.3,2X),D12.6),/54X,D12.6) 

RETURN 
END 

********************************************************  ********** 
SUBROUTINE  CROSS(A,B,C) 

SUBROUTINE  CROSS  NOVEMBER  1978  S.HELLINGER 

THE  SUBROUTINE  COMPUTES  THE  CROSS  PRODUCT  CF  TWO  VECTORS  AS  AXB=C. 

OOUBLE  PRECISION  A { 3 ) , 3 l 3 ) , C ( 3 ) 

C 

CI1)=AI2) *B!3)-A(3)*B!2) 
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C(2l=A(3l*Bl 1 >  —  A  1 1  )*B  ( 3 ) 

Cl3)=Atl)*R(2) -A (2)*Btl) 

RETURN 

END 

******£*********************************************************** 
SUBROUTINE  MM33(A,B,C) 

SUBROUTINE  MM33  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  MULTIPLIES  3*3  MATRICES  AS  A*B=C. 

DOUBLE  PRECISION  A (3 ,3) . B( 3,3) ,C (3,3) 

DO  10G  1-1,3 
DO  100  J  =  1 ,3 
C(I,J>=0. 

DO  100  K=1 *3 

100  C(I,J)=C(I,J!»A(I,KI*B(K,J) 

RETURN 

END 

****** ************************************************************ 
SUBROUTINE  RINVERT(R,P) 

SUBROUTINE  RINVERT  JANUARY  1979  C.TAPSCOTT 


THE  SUBROUTINE  CALCULATES  THE  ROTATION  AXIS  AND  ANGLE  OF  THE 
ROTATION  MATRIX  R. 

IMPLICIT  DOUBLE  PRECISION  |A-H,0-Z> 

DIMENSION  R(  3, 3 )  ,  P  ( 3  ) 

R90-0ASINII.) 

CA=(R(l»l)*-R(2»2)«-RI3»3)-l.)/2. 

IF  (OABS(CA)  .GT.  1.1  CA=DSIGN( 1 . ,CA ) 

IF  (Cl. -CA)  .GT.  (1. 50-10))  GO  TO  100 

THETA=R90 

PHI*0. 

ALPHA=0. 

GO  TO  10O0 


208 


100  CTMR13.  i  >-CA)/(l.-CA) 

If-  (CT  .LT.  0.)  cr=o. 

IF  (CT  .oT.  1.)  CT -I . 

CT  =  DS0KT 1 CT) 

THETA=0ACOS(CT  > 

C 

IF  (THETA  .CT.  1. 750-05)  GO  TO  200 

PHl=0. 

SA=(R(2tl)-R(lf2))/(2.*CT) 

GO  TO  9iJ 
C 

2CJ  IF  (CT  .CT.  1.750-05)  GO  TO  300 

IF  ( (CA  +  1 •  )  .LT.  1. 50-10)  GO  TO  250 

C 

SP=R(1.3 J-R(3, 1) 

CP-R  (  3»2)-R(2#3) 

IF  (SP  .EU.  0.1  PH  1  =  R90-QS IGN ( R9 0 »CP ) ;  GO  TO  201 
IF  (CP  .eQ.  0.1  PH  1=  OS  I GN ( R9  0#  SP  )  ;  GO  TO  202 
PHl=DATAiM2(SPtCP) 

IF  (DABS ( CP)  .LT.  DABS (S  P) 1  GO  TO  202 

201  SA=CP/(2.*DSIN(PHI  >*DSIN (THETA) ) 

GO  TO  903 

202  SAsSP/(2. *DCOS (PHI )*0S IN (THETA) ) 

GO  TO  900 

C  , 

250  THET A=R93 

ALPHA=2.*R90 

CP=(R(i,im.)/2. 

IF  (CP  .LT.  0.1  CP=0. 

IF  (CP  .GT.  1.1  CP=1. 

CP*DSQRT(CP) 

PMI-0SI GN  (  OACOS  (  CP  )  »  (R(  1  *2  )  «-R  (  2 1  1))) 

GO  TO  1000 
C 

3 00  CP*(«Ul#3)fR(3.1)  )/(2.*(  l.-CA)) 
SP*(R(2»3)»R(3t2))/(2.*I l.-CA) ) 

IF  (SP  .t C.  0.1  PH I*R90-DSIGN(R90»CP ) ;  GO  TO  301 
IF  (CP  .cC.  0.)  PH  1=  OS IGN ( R90 »  SP ) ;  GO  TO  301 
PHI=DATAN2($P,CP) 

C 

301  A»K(2. 1)-R(1.2) 

B=Ri It  3)-R(3,l) 

CSR ( 3t  2 )-R  (2 » 3  ) 

IF  (DABS(A)  .LT.  DABS! B ) )  GO  TO  302 
IF  ( DABS ( A )  .LT.  OABS(C))  GO  TO  303 
SA=A/(2.*CT) 

GO  TO  903 

3C2  IF  (OABS(B)  .LT.  OABS(CI)  GO  TO  303 
SA«8/(2.*0SIN(PHII *OSIN( THETA) > 

GO  TO  900 


nnoooo  o  no 
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3C3  S A=C/( 2 • *DC0S ( PH  I  )  *DS I N(  THETA)  ) 

C 

9 CO  IF  (SA  .EQ.  0.)  ALPH A=R90~DS IGN { R9Q, CA ) ;  GO  TO  1000 
IF  (CA  .EO.  0.)  ALPHA=0SIGN(R90, SA);  GO  TO  1000 
ALPHA3 JAT  AN2 ( SA »  CA ) 

C 

1000  P(1)=90.-(THETA*90./R90) 

P(2)=PHI *90. /R90 
P(3)-ALPHA*90./R90 
C 

RETURN 
ENO 

******-***  **  *  ********  ************  ********  ****  ********  ****  ********** 

SUBROUTINE  RMAT (Q»R) 

SUBROUTINE  RMAT  JANUARY  1979  C. TAP  SCOTT 

THE  SUBROUTINE  COMPUTES  THE  ROTATION  MATRIX  R  ASSOCIATED  WITH  THE 
ROTATION  Q  =  ILAT.  LONG.  ANGLE). 

IMPLICIT  DOUBLE  PRECISION  (A-H.O-Z) 

C 

DIMENSION  0(3) ,R (3,3) 

DATA  RAD/ 1.745329251994330-02/ 

C 

T= ( 90. -0 ( 1 ) ) *R  AD 
PSQ 1  2) *RAD 
A=Q(3)*RA0 
C 

CT=DCOS(T) 

ST=DSIN(T ) 

CP=DCOS(P ) 

SP®DSI NIP ) 

CA»OCOS(A) 

SA=DSI M( A ) 

W  A= 1 ,-C A 

RU,ll3CAMWA*CP*CP*ST*ST) 

R(2»2)=CAMWA*SP*SP*ST«'ST) 

R(3»3)=CAMWA*CT*CT) 

C 

U=WA*CP*SP*ST*ST 

V»CT*SA 

R(l,2)-U-V 

R(2.1)3U*V 

C 

U-WA*CP*CT*ST 

V«SP*ST*SA 


ooooo  ooo 


mmmm 


ill 


fa- 
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R<l,3J=UfV 
R l 3  » 1 ) =U- V 
C 

U=WA*$P*C I *ST 
v=cp*st+sa 

R  <2,3)=U-V 
R  (  3*2)  =U»- V 
C 

RETURN 
END 

**4**»*«.^  **«***$**«»  **$*****♦****$***$**********  **********«*#*$*** 

SUBROUTINE  ROT AT E < R , V , VR ) 

SUBROUTINE  ROTATE  JANUARY  1979  C.TAPSCOTT 

THE  SUBROUTINE  RATATES  THE  VECTOR  V  TO  VR  BY  THE  MATRIX  R. 

DOUBLE  PRECISION  R 1 3 , 3 J , V( 3) ,VR< 3 1 
C 

DU  1GU  1=1,3 
VRl I)=0. 

DO  ICO  J  = 1  *  3 

ICO  VRm  =  VRU)HlII,  J)*V(J) 

C 

RETURN 

END 
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